THE NONLINEAR STABILITY OF THE IRROTATIONAL EULER-EINSTEIN SYSTEM 
WITH A POSITIVE COSMOLOGICAL CONSTANT 



IGOR RODNIANSKI AND JARED SPECK 



Abstract. In this article, we study small perturbations of the family of Friedmann-Lemaitre-Robertson-Walker 
cosmological background solutions to the Euler-Einstein system with a positive cosmological constant in 1 + 3 
dimensions. The background solutions describe an initially uniform quiet fluid of positive energy density evolving 
in a spacetime undergoing accelerated expansion. We show that under the equation of state p — c 2 s p, < c 2 s < 1/3, 
the background solutions are globally future asymptotically stable under small irrotational perturbations. In 
particular, we prove that the perturbed spacetimes, which have the topological structure [0, oo) x T 3 , are future 
causally geodesically complete. 
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1. Introduction 

The irrotational Euler-Einstein system models the evolution of a dynamic spacetimd3 (Ai,g) containing a 
perfect fluid with vanishing vorticity. In this article, we endow this system with a positive cosmological constant 
A and consider the equation of state p = c 2 p, where p is the fluid pressure, p is the proper energy density, and 
the non-negative constant c s is the speed of sound. As is fully discussed in Section [3l under these assumptions, 
the irrotational Euler-Einstein system comprises the equations!! 

(1.0.1a) R^-^Rg^ + Ag^=Tj^ alar \ Qu,*/ = 0,1,2,3), 

(1.0.1b) D a (a s D a ®) = 0, 

where is the spacetime metric, D is the covariant derivative operator corresponding to g^ u , R^ u is the Ricci 
curvature tensor, R = g a @ R a p is the scalar curvature, <I> is the fluid potential j^ ca ' ar ) = 2a s {d^){d v ^) + 
g^v{ s + l) _1 o" s+1 is the energy-momentum tensor of an irrotational fluid, a = —g a P(d a &)(di3&) is the enthalpy 
per particle, and s = (1 — c 2 )/(2c 2 ). The fundamental unknowns are (A4,g,d&), while the pressure and proper 
energy density can be expressed as p = j?p[0~ s+1 , p = 2 g^ cr s+1 . In this article, we will mainly restrict our 

attention to the case s > 1, which is equivalent to < c s < a/1/3. Although we limit our discussion to the 
physically relevant case of 1 + 3 dimensions, we expect that our work can be easily generalized to apply to the 
case of 1 + n dimensions, n > 3. 

Remark 1.0.1. Due to topological restrictions arising in the application of Poincare's lemma (see Section 
I3.1|) . the function may only be defined locally. However, only its spacetime partial derivatives d& enter into 
equations (jl.O.lap - (jl.O.lbp and into our estimates. Since the quantities that can be locally expressed as d& 
do not suffer from the same topological obstructions as $, we abuse notation and use the single symbol "<9<I>" 
to denote them globally. In fact, we could in principle avoid introducing d<&, choosing instead to work with the 
1-form see equation (|3.1.5|) . The advantage of working with d& is that the wave equation (jl.O.lbp allows us 
to quickly and compactly obtain energy estimates using a standard integration by parts argument (see Section 

The choice of an equation of state is necessary in order to close the relativistic Euler equations. Our choice 
of p = c 2 p is often made in the mathematics and cosmology literature. As is explained in Section [H under such 
equations of state, there exists a family of Friedmann-Lemaitre- Robertson- Walker (FLRW) solutions to (|1.0.1a|) 
- (jl.O.lbp that are frequently used to model a fluid-filled universe undergoing accelerated expansion; these are 
the solutions that we investigate in detail in this article. The cases p = and p = (l/3)p, which are known as 
the 'dust" and "radiative" equations of state, are of special significance in the astrophysical literature. The latter 
is often used as a simple model for a "radiation-dominated" universe, while the former for a "matter-dominated" 
universe. Unfortunately, as we will see, these two equations of state lie just outside of the scope of our main 
theorem. Our results can be summarized as follows. We state them roughly here; they are are stated more 
precisely as theorems in Sections [TT] and [12J 

Main Results. If < c s < ^JT/3 (i.e. s > 1), then the FLRW background solution ([0, oo) x T 3 ,g, 5$) 
to (jl.O.lap - (jl.O.lbp . which describes an initially uniform quiet fluid of constant positive proper en- 
ergy density, is globally future-stable under small perturbations. In particular, small perturbations 
of the initial data corresponding to the background solution have maximal globally hyperbolic devel- 
opments that are futurdj causally geodesically complete. Here, g = —dt 2 + e 2 ^ 1 ^ Yli=i(d x3 ) 2 > an d 
<9$ = f (d t $,d 1 <!t>,d2<S>,d 3 <S>) = (^e-^W, 0,0,0), where * is a positive constant, x = 3/ (2s + 1) = 3c 2 , 
and Q(t) defined in (I4.1.9P - Furthermore, in the wave coordinate system introduced in 



By spacetime, we mean a 4 —dimensional Lorentzian manifold M together with a metric g M „ on M of signature ( — , +, +, +). 
2 Throughout the article, we work in units with 8nG = c = 1, where c is the speed of light and G is Newton's universal 
gravitational constant. 

"^In [10] . Christodoulou refers to $ as the wave function. 
^Throughout this article, dt is future-directed. 
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Section [5.1[ the components g^ u of the perturbed metric, its inverse g^ u , the fluid potential spacetime 
gradient d&, and various coordinate derivatives of these quantities converge as t — > oo. 

Remark 1.0.2. We made the irrotationality assumption only for simplicity. In a future article, we will remove 
this restriction. 

Remark 1.0.3. The kind of stability result stated above, in which the solution is shown to converge as t — >■ oo, 
is known as asymptotic stability. 

Remark 1.0.4. Note that our results only address perturbations of fluids featuring a strictly positive proper 
energy density. We have thus avoided certain technical difficulties, such as dealing with a free boundary, that 
arise when p vanishes. 

Remark 1.0.5. We have only shown stability in the "expanding" direction (t —> oo). 

Remark 1.0.6. We define T d = f [— ir, n] d with the ends identified. 

We would now like to make a few remarks about the cosmological constant. We do not attempt to give 
a detailed account of the rich history of A, but instead defer to the discussion in [5]; we offer only a brief 
introduction. While the cosmological constant was originally introduced by Einstein |14) to allow for static 
solutions to the Einstein equations in the presence of matter, the present day motivation for introducing A > 
is entirely different. In 1929, Hubble discovered of the expansion of the universe |20j. In brief, Hubble's 
"law," which was formulated based on measurements of redshift, states that the velocities at which distant 
galaxies are receding from Earth are proportional to their distance from Earth. Furthermore, the present day 
explanation is that the cause of these velocity shifts is the expansion of spacetime itself. For example, a metric 
of the form g = —dt 2 + a 2 (t) Ylf=i(dx 1 ) 2 , with a > 0, creates a redshift effect. In the 1990's, experimental 
evidence derived from sources such as type la supernovae and the Cosmic Microwave Background suggest that 
the universe is in fact undergoing accelerated expansion. Our main motivation for introducing the positive 
cosmological constant is that it allows for spacetime solutions of (|1.0.1a|) - (jl.O.lbp that feature this effect. A 
simple example of a solution to the Einstein-vacuum equations that features such accelerated expansion is the 
metric g = —dt 2 + e 2Ht Y^=i(dx 1 ) 2 on the manifold (— oo, oo) x T 3 , where H = a/A/3. 

The introduction of a positive cosmological constant is not the only known mechanism for generating so- 
lutions to Einstein's equations with accelerated expansion. In particular, Ringstrom's work |28| . which is the 
main precursor to this article, shows that the Einstein-scalar field system, with an appropriate nonlinear- 
itj@ V($), has an open family of solutions with accelerated expansion. More specifically, the system stud- 
ied by Ringstrom can be obtained by replacing (jl.0,lb|) with g a @ D a Dp§ = V'(<&) and setting T^alar) _ 
(<9 /t <I> )(<%<!>) — [^g a/3 (d a ^)(d^^) + V(&)]g^ in equation (ll.O.laj) . where D denotes the covariant derivative 
induced by g. V is required to satisfy V(0) > 0, V(0) = 0, V" (0) > 0, so that in effect, the influence of the 
cosmological constant is emulated by V(&). His main result, which is analogous to our main result, is a proof 
of the future global stability of a large class of spacetimes featuring accelerated expansion. 

The Main Results stated above allude to both the existence of an initial value problem formulation of the 
Einstein equations, and the existence of a "maximal" solution. These notions are fleshed out in Section 13.21 
but we offer a brief description here. One of the principal difficulties in analyzing the Einstein equations is 
the lack of a canonical coordinate system. Intimately connected to this difficulty is the fact that due to the 
diffeomorphism invariance of the equations, their hyperbolic nature is not readily apparent until one makes some 
kind of gauge choice. One way of resolving these difficulties is to work in a special coordinate system known as 
wave coordinates (also known as harmonic gauge or de Bonder gauge), in which the Einstein equations become 
a system of quasilinear wave equations. One advantage of such a formulation is that local-in-time existence for 
a system of wave equations is immediate, because a standard hyperbolic theory based on energy estimates has 
been developed (consult e.g. \19\ Ch. VI], [34|, Ch. 16], |30| . |31|). Although the use of wave coordinates is often 
attributed solely to Choquet-Bruhat, it should be emphasized that use of wave coordinates in the context of 
the Einstein equations goes back to at least 1921, where it is featured in the work of de Donder |13| . However, 



Ringstrom refers to V($) as a "potential." However, we call V($) a "nonlinearity" in order to avoid confusing it with our 
"potential" <E>, and in order to avoid confusing it with a term of the form V(x)&, which is often referred to as a "potential" term in 
the literature. 
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the completion of the initial value formulation of the Einstein equations is in fact due to Choquet-Bruhat [8], 
her main contribution being a proof that the wave coordinate condition is preserved during the evolution if it 
is initially satisfied. 

The initial data for the irrotational Euler-Einstein system consist of a 3— dimensional Riemannian manifold 
S, a Riemannian metric g on S, a symmetric two-tensor K on S, and initial data <9<i>, \& for the tangential and 
normal derivatives of <I> along S. A solution consists of a 4— dimensional manifold Ai, a Lorentzian metric g 
and the spacetime derivatives d& (see Remark 11.0. lj) of a scalar field on A4 satisfying (jl.O.lap - (ll.O.lbJ . and 
an embedding S C A4 such that g is the first fundamental formal of E, K is the second fundamental forrrQ of E, 
and the restriction of 9$ and to E are 9<l and ^ respectively. Here denotes the derivative of in 

the direction of the future-directed unit normal N to E. It is important to note that the initial value problem 
is overdetermined, and that the data are subject to the Gauss and Codazzi constraints: 

(1.0.2a) R - K ah K ab + (g ab K ab ) 2 = Ao s & - -^-a s+ \ 

s + 1 

(1.0.2b) D a K aj - g ab DjK ab = 2a s 4fd j $, (j = 1, 2, 3), 

where R is the scalar curvature of g, D is the covariant derivative induced by g, and a = <t|vj. 

Remark 1.0.7. In this article, we do not address the issue of solving the constraint equations for the system 
(fTTTTal) - (ll.O.lbp . 

17 years after the initial value problem formulation was understood, Choquet-Bruhat and Geroch showed [7] 
that every initial data set satisfying the constraints launches a unique maximal globally hyperbolic development. 
Roughly speaking, this is the largest spacetime solution to the Einstein equations that is uniquely determined 
by the data. This result is still a local existence result in the sense that it allows for the possibility that the 
spacetime might contain singularities. In particular, future-directed, causal geodesies may terminate, which 
in physics terminology means that an observer (light ray in the case of null geodesies) may run into the end 
of spacetime in finite affine parameter. For spacetimes launched by initial data near that of the background 
solution, our main result rules out the possibility of these singularities for observers (light rays) traveling in the 
"future direction." 

We offer a few additional remarks concerning wave coordinates. The classic wave coordinate condition is 
the algebraic relation = 0, where the are the contracted Christoffel symbols of the spacetime metric. 
In this article, we use a version of the wave coordinate condition that is closer in spirit to the one used by 
Ringstrom in |28| . which was itself inspired by the ideas in |17| . Namely, we set = T^, where is the 
contracted Christoffel symbol of the background solution metric. Simple computations imply that = 3uj5q, 
where to = u(t) ~ yA/3 is defined in (|4.1.3p . It follows that in our wave coordinate system, the (geometric) 
wave equation g a P D a Dpv = for the function v is equivalent to the modified (also known as the "reduced") 
wave equation g a ^d a dpv = 3ui(dtv) 2 , which features the dissipative source term u(dtv) 2 . We provide a more 
detailed discussion of this modified scalar equation in Section 11.21 Furthermore, in Section [5J we modify the 
irrotational Euler-Einstein system in an analogous fashion, arriving at an equivalent hyperbolic system featuring 
dissipative terms. More precisely, the modified system is equivalent to the Einstein equations if the data satisfy 
the Einstein constraint equations and the wave coordinate condition. 

1.1. Comparison with previous work. First, it should be emphasized that the behavior of solutions the 
fluid equation (ll.O.lbp on exponentially expanding backgrounds is quite different than it is in flat spacetime. In 
particular, if one fixes a background metric on [0, oo) x T 3 near g, then our proof shows that the fluid equation 
(jl.O.lbp on this background with < c s < -y/l/3 has global solutions arising from data that are close to that of 
an initial uniform quiet fluid state, which is represented by d&. This is arguably the most interesting aspect of 

^Recall that the first fundamental form of E is defined to be the metric g on E such that g\ x (X, Y) — g\ x (X, Y) for all vectors 
X,Y e T X E. 

^Recall that K is defined at the point x by K(X,Y) ^ g(DxN,Y) for all X, Y £ T^E, where N is the future-directed unit 
normal to S at i. 
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our main result. In contrast, Christodoulou's monograph |10| shows that on the Minkowski space background, 
shock singularities can form in solutions to the irrotational fluid equation arising from data that are arbitrarily 
close to that of a uniform quiet fluid state. Our original intuition for this article was that rapid spacetime 
expansion should smooth out the fluid and discourage the formation of shocks. 

Second, we would like to compare our nonlinear stability result to the well-known linear instability results of 
Sachs and Wolfe [29]. In this work, they consider the Euler-Einstein system with A = under the equations 
of state p = and p = ^p, which correspond to borderline cases of the equations of state that we consider 
in this article. Sachs- Wolfe then consider a family of background solutions to this system on the manifold 
(— oo, oo) x M 3 . We remark that these well-known background solutions are of FLRW type, and can be obtained 
as special casea^ of solutions we have presented in Section HJ The first main result of |29] is that the linearized^] 
system is unstable; i.e., it features solutions that grow like t , where C < 1 depends on the equation of state. 
Because our theorem does not cover the borderline cases p = 0,p = -^p, it is not yet entirely clear whether or 
not the Sachs- Wolfe instabilities are caused exclusively by the absence of a cosmological constant, or whether 
the borderline equations of state also play a role; see Section 11.21 for some additional remarks concerning the 
borderline cases. As a side remark, we mention the most well-known result of the Sachs- Wolfe paper: they 
showed that the growing density perturbations couple back into the metric. The resulting variations in the 
metric lead to anisotropies in the Cosmic Microwave Background (CMB). In particular, the amount by which 
photons are gravitationally shifted varies with direction in the sky. The theoretical predictions of this effect, 
which is known was the Sachs- Wolfe effect, are consistent with the variations in the CMB detected by the 
Mather-Smoot team's COBE satellite in 1992 [3]. 

Next, we note that Brauer, Rendall, and Reula have shown [3] a Newtonian analogue of our main result. 
More specifically, they studied Newtonian cosmological models^ with a positive cosmological constant and with 
perfect fluid sources under the equation of state p = Cp]^ ewt , where pNewt > is the Newtonian mass density, 
and 7 > 1. They showed that small perturbations of a uniform quiet fluid state of constant positive density 
lead to a global solution. It is of particular interest to note that they do not require the fluid to be irrotational. 
This suggests that our main result can be extended to allow for (small) non-vanishing vorticity. As discussed 
in Remark 11.0.21 we will address this issue in an upcoming article. 

We also note a curious anti-correlation between our results and some well-known stability arguments for the 
Euler-Poisson system (a non-relativistic system with vanishing cosmological constant) which may be found e.g. 
in Chapter XIII of Chandrasekhar's book [6]. Chandrasekhar considers a simple model for an isolated body in 
equilibrium, namely a static compactly supported solution to the Euler-Poisson equations under the equation 
of stated p = Cn 7 , where n denotes the fluid element number density, C > is a constant, and 7 > 0. He uses 
virial identity arguments to suggest that such a configuration is stable if 7 > 4/3. However, since (13.0. 7ft implies 
that the equation of state p = c 2 s p (here p denotes the proper energy density, a relativistic quantity) is equivalent 
to p = Cn 1+Cs , our main results show that our background solution is stable under irrotational perturbations if 
1 < 1 + (? s < 4/3; i.e., our results seem to anti-correlate with the aforementioned non-relativistic onj^l. 

In addition to the previously mentioned work of Ringstrom, we would also like to mention other contributions 
related to the question of nonlinear stability for the Einstein equations with a positive cosmological constant. 
The first author to obtain stability results in this direction was Helmut Friedrich, first in vacuum spacetime |15] 
in 1 + 3 dimensions, and then later for the Einstein-Maxwell and Einstein- Yang-Mills equations |16] , Anderson 
then extended the vacuum result to cover the case of 1 + n dimensions, n odd [I]. Their work was based on the 
conformal method, which reduces the question of global stability for the Einstein equations to the much simpler 



^The Sachs- Wolfe solutions have spatial slices diffeomorphic to R 3 , while our solutions have spatial slices diffeomorphic to T 3 . 

9 That is, the system formed by linearizing the full Euler-Einstein equations around one of the FLRW background solutions. 
Their models were based on Newton-Cartan theory, which is a slight generalization of ordinary Newtonian gravitational theory 
that can be endowed with a highly geometric interpretation. 

^More correctly, his equation of state is p = Cp 1 Newt . However, in the case of the Euler-Poisson equations, pNewt is proportional 
to the number density n; i.e., p = Cn 7 is equivalent to p = CpJ iewt . 

12 We temper this statement by noting that our problem differs in several key ways from that of Chandrasekhar; e.g., Chan- 
drasekhar studied compactly supported data for a non-relativistic system on a flat background, while here we study relativistic 
fluids of everywhere positive energy density on an expanding background. 
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question of local-in-time stabilit\T 3 l for the conformal field equations, which were developed by Friedrich. In 
contrast, Ringstrom has stated that one of his main motivations for his wave coordinate approach in |28| is that 
the conformal method cannot necessarily be easily adapted to handle matter models other than Maxwell and 
Yang-Mills fields. Our work can be viewed as an example of the robustness of his methods. 

Finally, we compare our work here to the body of work on the stability of Minkowski space, which is the 
most well-known solution to the Einstein-vacuum equations in the case A = 0. This groundbreaking work, 
which was initiated by Christodoulou and Klainerman |12j , covered the case of the Einstein- vacuum equations 
in 1 + 3 dimensions. Their proof, which is manifestly covariant, relied upon several geometric foliations of 
spacetime, including maximal t = const slices and also a family outgoing null cones. In particular, it was 
believed that wave coordinates were unstable and therefore unsuitable for studying the global stability for 
Minkowski space. However, Lindblad and Rodnianski have recently devised yet another proof for the Einstein- 
vacuum and Einstein-scalar field systems |23| . |24| . which is much shorter but less precise, and which shows 
stability in the wave coordinate gauge = 0. In particular, Lindblad and Rodnianski were the first authors 
to show that a wave coordinate system can be used to prove global stability results for the Einstein equations. 
As we will explain in the next section, our result was technically simpler to achieve than either of these results. 
More specifically, in 1 + 3 dimensions with A = 0, the Einstein-vacuum equations involve nonlinear terms that 
are on the boardei0 of what can be expected to allow for global existence. As we will see, the addition of 
A > to the Einstein equations, together with our previously mentioned wave coordinate choice, will lead to 

the presence of energy dissipation terms. Consequently in the parameter range < c s < \J\, we do not have to 
contend with the difficult "borderline integrals" that appear in the proofs of the stability of Minkowski space. 
A more thorough comparison of the proofs of the stability of Minkowski space to the proofs of the stability of 
exponentially expanding solutions can be found in the introduction of |28| . Moreover, we remark that readers 
interested in results related to those of Christodoulou-Klainerman and Lindblad- Rodnianski can consult [2], 
[22], [25], [36]. 

1.2. Comments on the analysis. The basic idea behind our global existence argument can be seen by ana- 
lyzing the inhomogeneous wave equation g a ^D a D^v = F for the model metric g = —dt 2 + e 2t Yl^=i(^ xl ) 2 on 
the manifold-with-boundary Ai = [0, oo) x T 3 . Here, we are using standard local coordinates on T 3 . An omitted 
computation implies that relative to this coordinate system, the wave equation can be expressed as follows: 



(1.2.1) -8 2 v + e- 2t S ab d a d b v = 3(d t v) 2 + F. 

To estimate solutions to (jl.2.ip . one can define the "usual" energy E 2 (t) = ^ fj 3 (dtv) 2 + e~ 2t S ab (d a v)(dbv) d 3 x, 
and a standard integration by parts argument leads to the inequality 

(1.2.2) ± E <-E+\\F\\ L2 . 

It is clear from (jl.2.2p that sufficient estimates of ||-F 11^,2 in terms of E will lead to energy deca)£3- 

Our estimates for the modified irrotational Euler-Einstein equations are in the spirit of the above argument. 
The bulk of the work lies in estimating the inhomogeneous terms and in ensuring that the perturbed solution 
remains close to that of the background solution. We remark that the main tools used for estimating the 
inhomogeneous terms are Sobolev-Moser type estimates, which we have placed in the Appendix for convenience. 
Our analysis of the spacetime metric components closely parallels the work [28] of Ringstrom. In particular, 
based on Ringstrom's work, we provide (in Section [9]) a Counting Principle based on the number of spatial 
indices upstairs and downstairs in a product of metric components, Christoffel symbols, and first derivatives 
of these quantities. This heuristic mechanism can be used to quickly (and roughly, but good enough to prove 
small-data local existence) determine the rates of decay of products of such terms. Despite the availability of 

13 The conformal field equations are symmetric hyperbolic, and for such systems, local-in-time stability is a standard result. 

^They contain nonlinear terms that, on the basis of their order alone, might be expected to produce blow-up. However, these 
terms satisfy a version of the null condition, which means that they have a special algebraic structure that allows for small-data 
global existence. 

l^In our work below, we work with rescaled energies that are approximately constant in time. 
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this tool, we carefully derive many of the estimates in complete detail. We remark that the Counting Principle 
is not precise enough to detect the improved decay estimates derived in Section [12j 

Although Ringstrom's framework is useful for analyzing the metric components, our analysis of the fluid 
variables d^Q, {n = 0,1,2,3), involves additional complications beyond those encountered in his analysis. 
We would now like to briefly discuss these complications, and also to indicate why we make the assumption 

< c s < We believe that the breakdown of our proof in the case c s = is merely an artifact of our 
methods; we cannot address the equation of state p = here because we use the pressure as an unknown in 
our equations. In a future article, we will investigate this case using alternate methods. On the other hand, 
for a technical example of where the proof breaks down in the case c s > 1/3, see Section [11.31 We expect that 
stability can be shown for the equation of state p = (l/3)p using the Friedrich's conformal method; we plan to 
investigate this issue in a future article. On the other hand, the question of what to expect in the case c s > 1/3 
is open. 

The principal difficulty we encounter in our analysis is that our background solution fluid variable d& = 
(^e - ^, 0, 0, 0) and the associated quantity a = — g a P (d a Q)(d/3&) both exponentially decay to as t — > oo. 
Therefore, the quantity a corresponding to a slightly perturbed solution will also decay. Furthermore, by 
examining the fluid equation (ll.O.lbjl . we see that a = corresponds to a degeneracy. In order to deal with 
these difficulties, and in order to close our bootstrap argument for global existence, we have to ensure that a 
doesn't decay too quickly, and that it never becomes in finite time. In fact, for the global solutions that we 
construct, we show that the perturbed a will decay at the same rate as a. Since this fact plays a key role in 
our estimates for the fluid equation, we now heuristically outline our approach its proof. Assuming that the 
perturbed g^ u is near that of g^, (jj,, v = 0, 1, 2, 3), it follows that a « -\d t <S>) 2 + e~ 2n{ ^5 ab {d a ^){d b ^). Thus, 

to conclude that a ~ <j, we i) show that dt<& ~ dt&, and ii) introduce the quantities zj = e~ n ^dj&/dt&, 
(j = 1,2,3), which we assume are small initially, and then show that they remain small. In fact, to close our 
bootstrap argument, we show that they exponentially decay to in L°°. Combining i) and ii), we deduce that 

The above mathematical conditions have a physical interpretation. To elaborate, we first note that the 
four- velocity u of the fluid is connected to the fluid potential via (13,1, 3ft and (|3.1.5p . which imply that = 
—a~ l l 2 d[ jL <&. From this fact and the previous remarks, it thus follows that the smallness/decay condition on Zj 
ensures that the fluid decays towards the "low velocity" regime, in which the time component u° is dominant 
over the spatial components u 3 . 

Finally, to understand the assumption < c s < y/l/3, we have to explain some of the details of our estimates. 
If one tries estimate Zj in L°° using only energy estimates plus Sobolev embedding, then one can only conclude 
that ||j£j||loo is of size Ce during the time of existence when the rescaled energy is of size e. Unfortunately, this 
estimate is not strong enough to close our bootstrap argument. However, if we use the bootstrap estimated 
< Ce" n , \\dtdj$>\\i,°° < Cee~ xn , and the initial condition ||<9j<1>(0)||loc < Ce, then we can estimate 
by first estimating ||9j-$||x,oo through time integration using the bound on ||c?i<9j-<&||£oo. This leads to the 
estimate H^H^oo < Cee^ -1 '^*' , which decays in time if x < 1. The condition < x < 1 holds exactly when 

1 < s < oo, which is equivalent to < (? s < 1/3; this parameter range of stability is precisely the one mentioned 
in the Main Results statement above. 

1.3. Applications to spatial topologies other than T 3 . The model metric g = —dt 2 + e 2t Y^,=i(dx 1 ) 2 has 
another feature that is of crucial relevance for possible extensions of our work. To illustrate our point, let us 
consider g to be a metric on [0, oo) x M 3 , a Lorentzian manifold with boundary that has the Cauchy hypersurface 

S = f {t = 0|. Simple computations imply that the causal past of the causal future of a point intersects S in a 
precompaci 17 ] set. This is in stark contrast to the situation encountered in Minkowski space, where the causal 
future of a point is the forward null cone emanating from that point, and the causal past of this null cone 
includes the entire Cauchy hypersurface {t = 0}. One consequence of this fact is that the study of solutions to 
wave equations on exponentially expanding spacetimes is a "very" local problem; i.e., if we makes assumptions 

lf ^These estimates will follow easily from the assumption that our rescaled energy is of size e. 

-^For example, in this model spacetime, the causal past of the causal future of the origin is contained in the set {(t, x 1 , x 2 , x 3 ) | t > 

o,E- =1 (^) 2 <4}. 
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about the data in large enough ball in E, then we can control the solution in a non-compact region of spacetime 
that includes a cylinder of the form [0, oo) x B, where B is a spatial-coordinate ball. 

Using these observations, Ringstrom was able to prove the future-stability of various solutions to the Einstein- 
scalar field system for many spatial topologies in addition to T 3 . The main idea is to choose local coordinate 
patches on the spatial slices on which the problem is quantitatively close to the study of case of S = T 3 , and to 
piece together the future development of these patches into a global spacetime. The most difficult part of his 
argument is the global existence theorem on T 3 . However, his patching argument requires that one use cut-off 
functions, which introduces regions in which the Einstein constraint equations are not satisfied. To deal with 
this difficulty, he constructs his modified system of equations in such a manner that one can still conclude global 
existence, even if the constraint equations are not satisfied in the cut-off regions. Finally, after patching, these 
artificially-introduced regions are of course "discarded" and are not part of the spacetime. 

The modified system (I5.4.2ap - (15. 4, 3ft that we study is similar to Ringstrom's modified equations in that 
our global existence argument depends only on a smallness condition on the data, and not on whether or not 
the constraint equations are satisfied. As noted above, this is the main step in Ringstrom's work. For these 
reasons, it is very likely that his patching arguments can be used to extend our result to other spatial topologies. 
However, for sake of brevity, we do not explore this issue in this article. 

1.4. Outline of the structure of the paper. 

• In Section[2j we describe our conventions for indices and introduce some notation for differential operators 
and Sobolev norms of tensor-fields. 

• In Section [U we introduce the irrotational Euler-Einstein system. 

• In Section |4j we use a standard ODE ansatz to derive a well-known family of background Friedmann- 
Lemaitre-Robertson- Walker (FLRW) solutions to the irrotational Euler-Einstein system. 

• In Section [5] we introduce wave coordinates and use algebraic identities valid in such a coordinate 
system to construct a modified version of the irrotational Euler-Einstein system. We then discuss how 
to construct data for the modified system from data for the un-modified system. Finally, we discuss 
classical local existence for the modified system, including the continuation principle that is used in 
Section [TTJ 

• In Section [6j we introduce the relevant norms and the related energies for the modified system that we 
use in our global existence argument. We also provide a preliminary analysis of the derivatives of the 
energies, but the inhomogeneous terms are not estimated until Section [9j 

• In Section we introduce some bootstrap assumptions on the spacetime metric g^. We then use these 
assumptions to prove some linear-algebraic lemmas that are useful for analyzing g^ u and the inverse 
metric g^ v . 

• In Section [HJ we introduce our main bootstrap assumption, which is a smallness condition on Stv, a norm 
of difference between the perturbed solution and the background solution. We also define the positive 
constants q and rj m i n , which play a fundamental role in the technical estimates of the following sections. 

• Section [9] contains most of the technical estimates. We assume the bootstrap assumptions from the 
previous sections and use them to deduce estimates for g^ v , g^ v , m) iv , and for the nonlinearities appearing 
in the modified equations. Here, m fil ' denotes the reciprocal acoustical metric, which is the effective metric 
for the irrotational fluid equation. 

• Section [10] is a very short section in which we show that the Sobolev norms we have defined are equivalent 
to the energy norms. 

• In Section \TT\ we use the estimates from the previous sections to prove our main theorem, which is 
a small-data global (where "small" means close to the background solution) existence result for the 
modified equations. We then discuss the breakdown of our proof in the case c s > 1/3. Finally, we use 
the global existence theorem to prove a related theorem, which states that initial data satisfying the 
irrotational Euler-Einstein constraints, the wave coordinate condition, and a smallness condition lead to 
a future-geodesically complete solution of the irrotational Euler-Einstein system. 

• In Section PT2| we prove that the global solution from the main theorem converges as t — > oo. The main 
idea is that once we have a global small solution to the modified system, we can revisit the modified 
equations and upgrade some of the estimates proved in Section [TTJ 
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2. Notation 



In this section, we briefly introduce some notation that we use in this article. 

2.1. Index conventions. Greek indices a, (3, ■ ■ ■ take on the values 0, 1, 2, 3, while Latin indices a, b, ■ ■ ■ (which 
we sometimes call "spatial indices") take on the values 1, 2, 3. Pairs of repeated indices, with one raised and one 
lowered, are summed (from to 3 if they are Greek, and from 1 to 3 if they are Latin). We raise and lower 
indices with the spacetime metric g^ and its inverse g^ v . The only exceptions are equations (|1.0.2ap - (jl.0.2b[) 
and (|3.2.2ap - ()3.2.2bj) . in which we use the 3— metric gj^ and its inverse gi k to raise and lower indices, and in 
Section [12j in which all indices are raised and lowered with g^ u and g^ v except for the 3— metric g^^ , which 
has gj^\ as its corresponding inverse metric. 

2.2. Coordinate systems and differential operators. Throughout this article, we work in a standard local 
coordinate system x^x^x 3 on T 3 . Although strictly speaking this coordinate system is not globally well- 
defined, the vectorfields dj = f Tjfj are globally well-defined. This coordinate system extends to a coordinate 
system x°, x^x 2 , x 3 on manifolds with boundary of the form A4 = [0, T) x T 3 , and we often write t instead of 
x°. We write to denote the coordinate derivative ^J^, and we often write dt instead of 8q. Throughout the 
article, we will perform all of our computations with respect to the frame {<9^} _ 1 2 3 - 

If a = (n±, n2, n 3) is a triplet of non- negative integers, then we define the spatial multi-index coordinate 

differential operator d& by d& = f 9™ 1 d% 2 d% 3 . We denote the order of a by \a\, where |a| = n\ + ri2 + n^. 

We write D T Ul '" Ur — r) T Vv " Vr -4- V r V Ua r T Vl '" Va ~ lolVa+ll ' r — V s V a T Ul '" Ur tn denote the 

components of the covariant derivative of a tensorfield on Ai with components ■ 

We write ^T^:*; to denote the array containing of all of the N th order spacetime coordinate derivatives 
(including time derivatives) of the component T^'."^- Similarly, we write 0WT£:::£ to denote the array 
containing of all N th order spatial coordinate derivatives of the component T^'..,^" s ■ We omit the superscript 
when N = 1. 

2.3. Norms. All of the Sobolev norms we use are relative to the coordinate system x 1 , x 2 , x 3 introduced above. 
We remark that they are not coordinate invariant quantities. If T is a tensorfield on Ai with components T^'.'.'.^ r s 
relative to this coordinate system, then we define the following norm of T^'. '. '.u r s on t ne t = cord slices of Ai : 

(2.3.1) P%~z\\ s » iv) = ( E EE/ idzT^zox-^dxwy 2 . 

V Hl,-H s =0 vi,-v r =0 \a\<N Jx ^ J 

Note that the above norm is in general a function t. If the set U is omitted, then it is understood that 
U = T 3 . This is a slight abuse of notation since the coordinate system x , x 2 ,x 3 is not globally well-defined 
on T 3 . Furthermore, when writing integral expressions similar to the one on the right-hand side of (|2.3.ip . we 
suppress the coordinate map (x 1 , x 2 , x 3 ) : M ->• T 3 . This should cause no confusion since we are working relative 
to the frame {d^j _ Q 123 - 

Using the above notation, we can write the N th order homogeneous Sobolev norm of the component T^Zu r s 

as 



(2.3.2) \\d {N) TZ'Z\\v = E P* 1 ? 



dcf 

"^^-fJ.l—fJ.sWL 2 ' 

\a\=N 



If 8. C W 1 or ^CT™, then Cj^(&) denotes the set of A— times continuously differentiable functions (either 
scalar or array- valued, depending on context) on Int(^) with bounded derivatives up to order A that extend 
continuously to the closure of .ft. The norm of a function F G C]^(&) is defined by 

(2.3.3) |*W= E sup|fi r F(.)|, 

\I\<N 
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where dj- is a mult i- indexed operator representing repeated partial differentiation with respect to the arguments 
• of F, which may be either spacetime coordinates or state-space variables depending on context. When N = 0, 
we also use the notation 



(2.3.4) |F|*^sup|F(.)|. 

•ei? 

Furthermore, we use the notation 

(2.3.5) |i^U d = f 

\I\=N 

The quantity |-F^^ is a measure of the size of the N th order derivatives of F on the set 

If / C M is an interval and X is a normed function space, then we use the notation C N (I,X) to denote the 
set of Al-times continuously differentiable maps from I into X. 

2.4. Running constants. We use C to denote a running constant that is free to vary from line to line. In 
general, it can depend on N (see (|8.1.2|) ) and A, but can be chosen to be independent of all functions (g^ u , d^&), 
(/i, v = 0, 1, 2, 3), that are sufficiently close to the background solution (guv, d^<&) of Section^ We sometimes use 
notation such as C(N) to indicate the dependence of C on quantities that are peripheral to the main argument. 
Occasionally, we use c, C#, etc., to denote a constant that plays a distinguished role in the analysis. We remark 
that many of the constants blow-up as A — > + . 

2.5. A warning on the sign of Although we often choose notation that agrees with the notation used 
by Ringstrom in |28| . our reduced wave operator n g = g a ^d a dp has the opposite sign of the one in |28| . 

3. The Irrotational Euler-Einstein System 

The Einstein equations connect the Einstein tensor — ^g^R, which depends on the geometry of the 
spacetime manifold A4, to the energy-momentum-stress-density tensor (energy-momentum tensor for short) 
Tavi which models the matter content of spacetime. In 1 + 3 dimensions, they can be expressed as 

(3.0.1) - -Rg^v + Ag^ = T uv , (//, v = 0, 1, 2, 3), 

where R av is the Ricci curvature tensor, R is the scalar curvature, and A is the cosmological constant. We 
remark that the stability results proved in this article heavily depend upon the assumption A > 0. 

Recall that the Ricci curvature tensor and scalar curvature are defined in terms of Ruau , the Riemann curvature 
tenson^\. which can be expressed in terms of r" , the Christoffel symbols of the metric. In a local coordinate 
system, these quantities can be expressed as follows: 



(3.0.2a) R m J> d = f d a Tf v - d^£ v + T^Y^ - rf x T a \, (a, 0, fi, v = 0, 1, 2, 3), 

(3.0.2b) R av = R^ av a = d a T a a u — d^T^ v + T^ X T^ U — T^ X T^ V , (//,, u = 0,1, 2, 3), 

(3.0.2c) R= { g^R aP , 

(3.0.2d) r^/= \g aX (d^g Xu + d y g^-d x g^), (a, v = 0, 1, 2, 3). 

The Bianchi identities (see |35j ) imply that the left-hand side of (|3.0.ip is divergence free, which leads to the 
following equations satisfied by T^ u : 

(3.0.3) D v T> xu = 0, (/i = 0,1,2,3). 

By contracting each side of (|3.0,ip with g^ v , we deduce that R = 4A — T, where T == g a ^T a ^ is the trace of 
Tav- From this fact, it easily follows that f|3.0. lj) is equivalent to 

18 Under our sign convention, D^D v X a — D v D^X a — R^aXp. 
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(3.0.4) R„u = Ag^ + T ljlu --Tg lxl/ , (ji, v = 0, 1, 2, 3). 

It is not our aim to give a complete discussion of the notion of a perfect fluid. A thorough introduction to 
the subject, including a discussion of its history, can be found in Christodoulou's survey article [9]. We confine 
ourselves here to a brief introduction. In general, the energy-momentum tensor for a perfect fluid is 

(3.0.5) T^ luid) = {p + p)u»u» + PS T, (/x, u = 0,l, 2, 3), 

where p > is the proper energy density, p > is the pressure, and u is the four velocity, a unit-length (i.e., 
u a u a = — 1) future-directed vectorfield on M. The Euler equations, which are the laws of motion for a perfect 
relativistic fluid, are the four equations ([3.0.3P together with a conservation law ()3.0.6b|) for the number of fluid 
elements. In a local coordinate system, they can be expressed as follows: 



(0 = 0,1, 2, 3), 



(3.0.6a) D v T^ uid) = 0, 

(3.0.6b) D u {nu v ) = 0, 

where n is the proper number density of the fluid elements. We also introduce the thermodynamic variable r\, 
the proper entropy density, which we will discuss below. 

Unfortunately, even in a prescribed spacetime (M,g), the equations (|3.0.6ap - ()3.0.6b|) do not form a closed 
system. The standard means of closing the equations is to appeal the laws of thermodynamics, which imply the 
following relationships between the state-space variables (see e.g. |18j): 

(1) p > is a function of n > and rj > 0. 

(2) p > is defined by 

(3.0.7) p = „| 



where the notation |_ indicates partial differentiation with • held constant. 
(3) A perfect fluid satisfies 



(3.0.8) 



(3.0.9) 



dp 
dn 



>0, 



dp 
dn 



>0, 



dp 
On 



> with " =" iff rj = 0. 



As a consequence, we have that £, the speed of sound 19 ! in the fluid, is always real for n > : 

dp/dn^ 



£2 dg dp 
dp 



> 0. 



v dp/dn\ v 

(4) We also demand that the speed of sound is positive and less than the speed of light whenever n > and 
7] > 0: 

(3.0.10) n > and r? > =>- < C < 1. 

Postulates (1) - (3) express the laws of thermodynamics and fundamental thermodynamic assumptions, while 
Postulate (4) ensures that at each i£M, vectors that are causal with respect to the sound conCj in T x Ai are 
necessarily causal with respect to the gravitational null concF 1 ! in T x A4 . See |32| for a more detailed analysis of 
the geometry of the sound cone and the gravitational null cone. 

We note that the assumptions p > 0, p > together imply that the energy-momentum tensor (I3.0.5P satisfies 
both the weak energy condition (T^fl w X a X@ > holds whenever X is timelike and future-directed with 

respect to the gravitational null cone) and the strong energy condition ({Tj£, luid ^ 
holds whenever X is timelike and future-directed with respect to the gravitational null cone). Furthermore, if 



l/2g^T { J' uld) g^}X^X^ > 



19 In general, £ is not constant. However, for the equations of state we study in this article, £ is equal to the constant c s . 

2 ^The sound cone is defined to be the subset of tangent vectors X £ T X M such that m a pX a X 13 = 0, where m M „ is the acoustical 
metric. The matrix m p „ is the inverse matrix of the reciprocal acoustical metric m' 1 ", which is introduced in Section [5.41 i.e., m M „ 
is not obtained by lowering the indices of with g M „. 



21 



By gravitational null cone at x, we mean the subset of tangent vectors X £ T X M such that g ol /3X a X ,z 



0. 
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we assume that the equation of state is such that p = when p = 0, then f|3.0.9j) and (|3.0.10p guarantee that 
p < p. It is then easy to check that < p < p implies the dominant energy condition (—g^Tat^ 1 ^ X v is causal 
and future-directed whenever X is causal and future-directed with respect to the gravitational null cone). 

Postulate (1) from above is equivalent to making a choice of an equation of state, which is a function that 
expresses p in terms of r\ and p. An equation of state is not necessarily a fundamental law of nature, but can 
instead be an empirical relationship between the fluid variables. In this article, we consider the case of an 

irrotational, barotropic fluid under the equation of state p = <? s p, where < c s < and according to (I3.0.9|) . 
the constant c s is the speed of sound. A barotropic fluid is one for which p is a function of p alone. Because r\ 
plays no role in the analysis of such fluids, this quantity is absent from the remainder of our article. As discussed 
in the following section, these assumptions imply that the tensor r/j. Jui * defined in ([3.0.5P is equal to the tensor 

^{scalar) defined in (I3.3.7p . and that the equations (|3.0.3p are equivalent to (I3.3.6p . a single quasilinear wave 
equation for a scalar function $ (see Remark ll.O.ip . As a consequence, it follows that Tj"^* also satisfies the 
weak, strong, and dominant energy conditions. 



3.1. Irrotational fluids. In this section, we introduce the notion of an irrotational fluid. We assume that the 
fluid is barotropic, be we do not yet impose the particular equation of state p = c 2 s p. The scalar formulation of 
an irrotational fluid in a curved spacetime goes back to at least 1937 |33| . but in this article, we use modern 
terminology found e.g. in [10]. We begin by defining the quantity ^fa, which is known as the enthalpy per 
particle: 

(3.1.1) V^=^, 

n 

where a > 0. We note for future use that we can differentiate each side of (13.1.ip with respect to y/a, and use 
(I3.0.7P and the chain rule to conclude that 

In the above formula we are viewing p as a function of a. 
We also introduce the following one form: 

(3.1.3) /?/=-V^V, (^ = 0,1,2,3). 

The fluid vorticity v is then defined to be d(3, where d denotes the exterior derivative operator. In local 
coordinates, we have that 

(3.1.4) V = d/A-^, (H,v = 0,1, 2, 3). 

An irrotational fluid is defined to be a fluid for which Vuu vanishes everywhere. In this case, by Poincare's 
lemma, there locally exists (see Remark ll.O.ip a scalar function <J> such that 

(3.1.5) ^ = 9^, (^ = 0,1,2,3), 
which implies that the four velocity is connected to d& via the equation 

(3.1.6) u^-a-^d^, Gu = 0,l,2,3), 
and the enthalpy per particle is connected to d& via 

(3.1.7) a = -(D a $)(D a <S>) = -g a P(D a *)(Df,$). 

We now show that under the assumption of irrotationality, the equations (13.0.6ap reduce to a single quasilinear 
wave equation for <!>. We begin by identifying the Lagrangian for the purported wave equation with p : 



(3.1.8) 



JSf = JSf(o-) = p. 
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From (f3X2D and flgTHJ) , it follows that 

(3.1.9) " ^ 



<9cr a/ct p + p ' 

We also recall that the Euler-Lagrange equation for $ corresponding to the Lagrangian (|3.1.8p is 



(3.1.10) D a 



da 



0. 



Using (|3.1.6p and (|3.1.9p . we conclude that for an irrotational fluid, (|3.1.10p is equivalent to the continuity 
equation D v {nu v ) = from (|3.0.6bp . 

To show that the remaining equations (|3.0.6ap follow from (|3.1.10p , we first recall that the energy- momentum 
tensor for a Lagrangian scalar-field theory can be expressed as 



ft / 

(3-i.li) Tjf^ = -2— + 9u3 ,se, 



and that if d& is a solution to (|3.1.10p . then j^ ca ' ar ) [ s symmetric and divergence free: 

(3.1.12) = T U », (/*,*/ = 0,1, 2, 3), 

(3.1.13) D u T£ alar) =0, 0* = 0,1, 2, 3). 
For future use, we remark that if d& is a solution to the inhomogeneous equation 



(3.1.14) D c 
then 



oa 



+ Id<s, = 0, 



(3.1.15) D v T {Llar) = -2/8*^"*, (H = 0, 1, 2, 3). 
In the case of the Lagrangian (|3.1.8p , one can check that (|3.1.1ip is equivalent to 

(3.1.16) Tfr calar) = 2^-(d^)(d u $) + g^. 
Using (J3JL6D, (f3~L8l) . (13.1.111) . and (13.1.161) . we have that 

(3.1.17) T(f a ' ar) = ^=(d^){d^) + g^p = ( p + p) UflUu+ pg^. 



By examining (|3,0.5p . we observe that j^ ca ' ar ) = Tpl lmd \ To summarize, we have shown that if <9$ is a 
solution to ([3.1. 10p . and if p,n,u, and p are defined through <9<1? via equations (|3.1.8|) . (|3.1.2|) . (|3.1.6p . and 
(|3.1.9p respectively, then it follows all 5 equations from (|3.0.6ap - ()3.0.6bp are necessarily satisfied. Thus, it 
follows that for an irrotational fluid, the entire content of the Euler equations is contained in the single scalar 
equation (|3.1.10p . 

We conclude with a summary of the constraints that d& and £?(o~) must satisfy in order to have an irrotational 
fluids interpretation. First, since a > 0, (|3.1.7p implies that 

(3.1.18) - 5 ^(5 Q $)(%$) > 0. 
Similarly, since p is non-negative, (|3.1.2p implies that 

(3.1.19) JSf(cr) > 0. 

Next, (|3.1.2p and (|3.1.8p together imply that n = J^=, so that the requirement n > is equivalent to 
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(3.1.20) ^ > 0. 



Furthermore, since (|3.0.7p implies that ^-^(if/cr) = V&j^ ~~ P = Pi an d since p > 0, we must have 

< 3121 » 5^(1) 

Finally, since (j3~CT8|) implies that ^ > 0, and since (j3"X2]) and the chain rule imply that ^ = j^/^= 
n /(*7?(|^))' we must have 



(3.1.22) -—2 > 0. 

3.2. The initial value problem for the irrotational Euler-Einstein system. In this section, we discuss 
various aspects of the initial value problem for the Einstein equations, including the initial data and the notion of 
the maximal globally hyperbolic development of the data. We assume that we are given a Lagrangian if = if (a) 
that is subject to the constraints discussed at the end of the previous section; this is equivalent to fixing an 
equation of state p = p(p). We remark that the discussion in this section is very standard, and we provide it 
only for convenience. 

3.2.1. Summary of the irrotational Euler-Einstein system. We first summarize the results of the previous sections 
by stating that the irrotational Euler-Einstein system is the following system of equations: 

(3.2.1a) R ^-l Rg ^ + Ag^=T^ calar \ ( M , v = 0, 1, 2, 3), 



(3.2.1b) D, 



<9if 



0. 



da 

where if = if (a), a = -(D a $)(D a $), and T { ; u calar) = 2^(0 M *)(&$) + g^ . 

3.2.2. Initial data for the irrotational Euler-Einstein system. The initial value problem formulation of the Ein- 
stein equations goes back to the seminal work [8] by Choquet-Bruhat. Initial data for the system ([3.2. la|) - 
(|3.2.1bp consist of a 3— dimensional manifold £ together with the following fields on S : a Riemannian metric 
g, a covariant two-tensor K, and functions (9$,^. By d<&, we mean that these three components, which are 
required to be globally defined on S, should be the spatial gradient of functions $ that are locally defined on £ 
(see Remark ll.O.ip . 

It is well-known that one cannot consider arbitrary data for the Einstein equations. The data are in fact 
subject to the following constraints, where D is the covariant derivative induced by g : 

(3.2.2a) R - K ab K ab + (g ah K ab f = 2T^ calar) \^ 

(3.2.2b) D a K aj - g ab D 3 K ab = T^ calar) \^, (j = 1, 2, 3). 

We remark that in the case of the equation of state p = c 2 s p, the calculations in Section 13.31 will imply that 
T ( o Ca/ar) |s = 2<r s ^ 2 - ^a s+1 and T ( f aZar) | E = 2^6^, where a = & - g ab {d a $){d b 3>). 

The constraints (13.2.2aj) - (I3.2.2bp are manifestations of the Gauss and Codazzi equations respectively. These 
equations relate the geometry of the ambient Lorentzian spacetime (M,g) (which has to be constructed) to 
the geometry inherited by an embedded Riemannian hypersurface (which will be (S, g) after construction). 
Without providing the rather standard details (see e.g. |11|). we remark that they are consequences of the 
following assumptions: 

• S is a submanifold of the spacetime manifold M 

• g is the first fundamental form of S 
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• K is the second fundamental form of £ 

• The irrotational Euler-Einstein system is satisfied along £ 

• We are using a coordinate system (x = t^x 1 ,x 2 ,x 3 ) on Ai such that £ = {x G | i = 0}, and along 
£, 5oo = -1, 90j = 0, gjk = 9jk, d t gjk = 2K jk , (j, k = 1, 2, 3), <9$ = 5$, and = 

We recall that g is the Riemann metric on £ defined by 

(3.2.3) g\ x (X, Y) = g\ x (X, Y), VX, Y G T x £, 
and that K is the symmetric tensorfield on £ that can be expressed as 

(3.2.4) K\ X (X, Y) = g\ x (D x N, Y) = g\ x (D Y N, X) = ^(D^g)\ x (X, Y), VX, Y G T x £, 

where N is the future-directed normal to £ at x, and D is the covariant derivative induced by g. 

3.2.3. The definition of a solution to the irrotational Euler-Einstein system. In this section, we define the notion 
of a solution to the Einstein equations. We begin with the following definition, which describes the maximal 
region in which a solution is determined by its values on a set S. 

Definition 3.2.1. Given any set S G Ai, we define T>(S), the Cauchy development of S, to be the union 
T>(S) = T> + (S) U T>~(S), where D + (S) is the set of all points p G Ai such that every past inextendible causal 
curve through p intersects S, and T>~(S) is the set of all points p G Ai such that every future inextendibld^l 
causal curve through p intersects S. 

We also rigorously define a Cauchy hyper surf ace: 

Definition 3.2.2. A Cauchy hypersurface in a Lorentzian manifold Ai is a hypersurface that is intersected 
exactly once by every inextendible timelike curve in Ai. 

It is well-known that if £ C Ai is a Cauchy hypersurface, then P(£) = Ai (see e.g. |27|). 

Given sufficiently smooth initial data (£, g, K, d<&, as described in Section [3.2.2[ a (classical) solution to 
the system ()3.2.1ap - ()3.2.1bp is a 4 dimensional manifold Ai, a Lorentzian metric g, the spacetime derivatives 
9$ of a locally defined function $ (see Remark ll.O.ip . and an embedding^ £ C Ai subject to the following 
conditions: 

• g is a C 2 tensorfield and and d& is a C 1 tensorfield 

• £ is a Cauchy hypersurface in (Ai , g) 

• g is the first fundamental form of £ 

• K is the second fundamental form of £ 

• D^Q = *S>, where N is the future-directed normal to £. 

The triple (Ai,g,d&) is called a globally hyperbolic development of the initial data. 

3.2.4. Maximal globally hyperbolic development. In this section, we state a fundamental abstract existence result 
of Choquet-Bruhat and Geroch [7], which states that for initial data of arbitrary regularity, there is a unique 
"largest" spacetime determined by it. The following definition captures the notion of this "largest" spacetime. 

Definition 3.2.3. Given initial data for the irrotational Euler-Einstein system (|3.2.1ap - (|3,2.1bp . a maximal 
globally hyperbolic development of the data is a globally hyperbolic development (Ai,g, <9<3?) together with an 
embedding i : £ — > Ai with the following property: if (Ai' , g' , (?$') is any other globally hyperbolic development 
of the same data with embedding t' : £ — > Ai', then there is a map ip : Ai' — > Ai that is a diffeomorphism onto 
its image such that ip*g = g' , ip*d& = d&, and if) o i! = i. Here, ip* denotes the pullback by tp. 

Before we can state the theorem, we also need the following definition, which captures the notion of having 
two different representations of the same spacetime. 

22 Under the above assumptions, it follows that at every point i£S, iV^ = (1, 0, 0, 0). 

2 ^A curve 7 : [so,s mal ) — ¥ M is said to be future-inextendible if lim s _^ Smox does not exist. Past inextendibility is defined in an 
analogous manner. 

24 Here, we suppress the embedding t : E — > M, and identify E with t(E). 
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Definition 3.2.4. The developments (A4,g,d<&) and (M' , g' , d&) are said to be isometrically isomorphic if 
the map ip from the previous definition is a diffeomorphism from M to M! . 

We now state the aforementioned theorem. 

Theorem 3.1. (Choquet-Bruhat and Geroch, 1969) Given any initial data for the irrotational Euler-Einstein 
system (|3.2.1aj) - (|3.2.1b|) . th ere is a maximal globally hyperbolic development of the data which is unique up to 
isometry. 

3.3. Calculations for the equation of state p = c 2 s p. For the remainder of this article, we restrict our 
attention to equations of state of the form 

(3.3.1) P = c 2 s p, 

where by equation (|3.0,9p . c s is the speed of sound. As mentioned in the introduction to this article, our stability 
results are limited to the following parameter range: 

(3.3.2) < c s < Ji. 

Equations (|3.0.7p . (|3.1.8p . ()3.1.9|) . and f)3.3. 1[) imply that there exist constants C > and C > such that 

(3.3.3) p = Cn 1+C ° = Ca s+ \ 
where 

1 - <?. 2 _ i 



(3-3.4) .= -5^, <-2s + l- 

Choosing a convenient normalization constant, it follows that under the equation of state (|3.3.1|) . the Lagrangian 
(|3.1.8p is given by 

(3.3.5) JS? = (a + 

Furthermore, in the case of the Lagrangian (|3.3.5j) . the Euler-Lagrange equation ([3.1.10P is 

(3.3.6) D a {o s D a <&) = 0, 
while the energy-momentum tensor is given by 

(3.3.7) T(r Zar) = 2a s (d^){d v <S>) + g uu (s + 1)~ V s+1 . 

For future reference, we record here the following two identities, which follow easily from ([3.3.7P : 



(3 3 8) rp(scalar) SLi 1 g a P'J , v 



(scalar) 4£ f a/3rp(scalar) 



2(1 -s) 

S + l 



(3.3.9) rp(scalar) _ ^(scalar) ^ = 2a s (d,<Z>)(d u <t>) + (_!_) a^ 1 ^, (/i, V = 0,1, 2, 3). 

3.3.1. Summary of the irrotational Euler-Einstein system under the equation of state p = (? s p. To summarize, 
we note that under the equation of state p = c^p, the irrotational Euler-Einstein system comprises the equations 

(3.3.10a) iV - kg^ - T^ ca/ar) + ±r( flcaior > g» v = 0, (p, u = 0,l, 2, 3), 

where T%° alar) = 2a s (d fl <l>){d L ,<l>)+g fM ,(s + iy 1 a s+1 is as in (f3TB~7T) . together with (f3~3~6]) . the equation of motion 
for an irrotational fluid: 



(3.3.10b) 



D a (a s D a $) = 0. 



Stability of Irrotational Euler-Einstein 

16 

4. Background Solutions 

Our main results concern the future-stability (with respect to irrotational perturbations) of a class of back- 
ground solutions (JO, oo) x T 3 ,g,d(j)) to the system (|3.3.10ap - (|3.3.10bp . These background solutions, which 
are of FLRW typc3, physically represent the evolution of an initially uniform quiet fluid in a spacetime that is 
undergoing rapid expansion. To find these solutions, we follow a procedure outlined in chapter 5 of |35j which, 
under appropriate ansatzes, reduces the Euler-Einstein equations to ODEs. Although our goal is to find ODE 
solutions to the irrotational equations (13.3. lOaj) - (|3.3.10bj) . the procedure we follow will produce ODE solutions 
to the full Euler-Einstein system fj3.0. lf> + (I3.0.5p . However, these ODE solutions will turn out to be irrota- 
tional. Thus, as discussed in Section 13. 1\ these solutions can also be interpreted as solutions to the irrotational 
system. We remark that the derivation of these solutions is very well-known, and that we have provided it only 
for convenience. 

4.1. Derivation of the background solution. To proceed, we first make the ansatz that the background 
metric g is of the form 

3 

(4.1.1) g = -dt 2 + a 2 (t)^2(dx i ) 2 , 

i=l 

from which it follows that the only corresponding non-zero Christoffel symbols are 

(4.1.2) f/ fc = f k °j = aaS jk , f/ = f^- = (j, fc = 1,2,3), 
where 

(4.1.3) co^^, 

a 

and a = f ^a. Using definitions (|3,0,2bp and (|3,0.2cp . together with (|4.1,2p . we compute that 



(4.1.4a) 


Rqo ~ 


2^500 


(4.1.4b) 


R 0j - 


1 


(4.1.4c) 


Rjk - 




We then 


assume that p = 


p(t),p 



), (j = 1,2,3), 

-(2ad + a 2 )5j k , (j, k = 1, 2, 3). 

and = (1,0,0,0), which implies that <9<J> = (dt$(t), 0, 0, 0). We 
also assume that the equation of state ()3.3. 1[) holds. Inserting these ansatzes into the Einstein equations (|3.0.ip 
+ (|3.0.5p . we deduce (as in |35| ) the following equations, which are known as the Friedmann equations in the 
cosmology literature: 

(4.1.5a) pa 3 ^ = pd 3 ^ d = f k, 



A p A k 

(4.1.5b) a = a\ — + ^ = a\ — + 



3 ' 3 1 3 ' 3a 3 (!+ c i) ' 

where the positive constant p denotes the initial (uniform) energy density, and the positive constant d is defined 

by d = f a(0). Observe that the rapid expansion of the background spacetime can be easily deduced from the 

ODE ()4.1.5bp . which suggests that the asymptotic behavior a(t) ~ e Ht ,H = f ^/A/3. A more detailed analysis 
of a(t) is given in Lemma 14.2.1 i 

With of = f -g af3 {d a $){dp$) = (<9 t $) 2 , we use (|5T51l . ([5X4} . (l3~33]l . and (|4.1.5al) to compute that 



Technically, the term "FLRW" is usually reserved for a class of solutions that have spatial slices diffeomorphic to S 3 , R 3 , or 
hyperbolic space (see e.g. |35|L 
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clka^ 1 ^ 



(4.1.6) c' s Ka = c' s p = p 

The equalities in f|4. 1 .6j) imply that 
(4.1.7) 
where 



2c s ~(l+cj)/(2^) _ 2c s (g^q+M/tc 2 .) 

1 + cf l + cV 0t > 



(4.1.8) 
(4.1.9) 
(4.1.10) 



T dcf /_ S + l \V(2-+2), 



2s + 1. 
0(t)=ln (a(t)), 

def 3 



Remark 4.1.1. 0(f) has been introduced solely for cosmetic purposes. 

For future use, we note the following consequences of the above discussion: 



(4.1.11a) 
(4.1.11b) 



3oj 2 — A = p = ke 



2S + 1W 
s + lJ 



-3d) = 3u? — A + k( e- 2 * s +W = 3k( ^±1 1 « 

V2s + 1/ \2s + l' 



2x(s+i)n 



3a 



s + l 



4.2. Analysis of Friedmann's equation. In the following lemma, we analyze the asymptotic behavior of 
solutions to the ODE (14.1.5bj) . 

Lemma 4.2.1. Let a, k,P > 0, and let a(t) be the solution to the following ODE: 

d 



(4.2.1) 



'A k 

-a = a\l 1 5-. 

dt V 3 3a p ' 



o(0) = a. 



Then with H = y/A/3, the solution a(t) is given by 



(4.2.2) 



a(t) 



and for all integers N > 0, there exists a constant CV > such that for all t > 0, with 

2/P 



def I 1 



.4 

(4.2.3a) 
(4.2.3b) 



3W 



+ + a 



P/2 



, we have that 



(l/2) 2 ' p ae m < a(t) < Ae Ht , 



< C N e 



-PHt 



Furthermore, for all integers N > 0, there exists a constant Cn > such that for all t > 0, with 

def a 



(4.2.4) 
we have that 



(4.2.5a) 
(4.2.5b) 



H < u(t) < \ H 2 + 



dt N 



(t) - H) 



3a p: 
< C N e- pm . 
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Remark 4.2.1. Because of equation ()4.1.5bp . we will assume for the remainder of the article that P = 3(l + c^). 
Proof. We leave the elementary analysis of this ODE to the reader. □ 

5. The Modified Irrotational Euler Einstein System 

In this section, we introduce our version of wave coordinates. We then use algebraic identities that are valid 
in wave coordinates to construct a modified version of the irrotational Einstein equations, which are a system of 
quasilinear wave equations containing energy-dissipative terms. Next, to facilitate our analysis in later sections, 
we algebraically decompose the modified system into principal terms and error terms. Finally, we show that 
the modified system is equivalent to the un-modified system if the Einstein constraint equations and the wave 
coordinate condition are both satisfied along the Cauchy hypersurface E. 

5.1. Wave coordinates. To hyperbolize the Einstein equations, we use a version of the well-known family 
of wave coordinate systems. More specifically, we use a coordinate system in which the contracted Christoffel 
symbols = f g a @T ^ a of the spacetime metric g are equal to the contracted Christoffel symbols = f g a ^^a(S 
of the background metric g. This condition is known as a wave coordinate condition since = if and only 
if the coordinate functions I are solutions to the wave equation g^D^D^ + 1^ = 0. Using fliXT]) and 
(|4.1.2p . we compute that in such a coordinate system 



(5.1.1) r" = r> = 3u<, r M = 9lia r a = 3^, (a* = 0,1,2,3), 

where co(t) is defined in (|4.1.3|) . 
We now introduce the tensorfields 



(5.1.2a) P" = 3a;$, P M = 3u; ff o M , (a* = 0, 1,2,3), 

(5.1.2b) Q» d £p»- T », Q» = P»-T^ (a* = 0,1, 2, 3). 

The idea behind wave coordinates is to work in a coordinate system in which = 0, so that whenever it is 
convenient, we may replace with 3lu5q (and vice-versa) without altering the content of the Einstein equations. 
The existence of such a coordinate system is nontrivial, and it was only in 1952 that Choquet-Bruhat [8j first 
showed that they exist in general. With this idea in mind, we define the modified Ricci tensor R pv by 

(5.1.3) R pu = R pu + ^ {D^Qv + D V QJ) 

1 1 

where 



(5.1.4) h g = g ap d a dp 

is the reduced wave operator for the metric g. 

We now replace the R pv with R pv in (|3.3.10ap . expand the covariant differentiation in ()3.3.10bp . and add 
additional inhomogeneous terms I pv and Iq§ to the left-hand sides of (|3.3.10ap and ()3.3.10b| ) respectively, 
arriving at the modified irrotational Euler-Einstein system: 

(5.1.5a) iV - A 9fa/ - T% calar ^ + ^ scala ^g pv + I pv = 0, (p, u = 0,l, 2, 3), 

(5.1.5b) [ag a/3 - 2s(d a <5>){d (3 <S>)]d a d [S <S> - aT a d a ^ + 2sF aXI3 (d a $)(d x $)(dp<f>) + J w = 0. 



'The £ M are scalar-valued functions, despite the fact that they have indices. 
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(5.1.6a) I 00 = -2wQ° = 2u(T° - 3u), 



(5.1.6b) I 0j = Ijo d = lujQj = 2uj(3(jjg 0j - Tj), (j = 1, 2, 3), 



(5.1.6c) / Jfe = / ifc d 4 f 0, (j,k = 1,2,3), 

(5.1.6d) Ig* = -ag a PQ a d p $> = aT a d a <t> - 3ujad t $. 

We have several important remarks to make concerning the modified system (15.1.5ap - (I5.1.5bj) . First, because 
the principal term on the left-hand side of (|5.1.5a|) is — ^Ugg^u, the modified system comprises a quasilinear 
system of wave equations and is of hyperbolic character. Second, the gauge terms 1^,10$ have been added 
to the system in order to produce an energy dissipation effect that is analogous to the effect created by the 
3(<9jt>) 2 term on the right-hand side of the model equation (jl.2.ip . These dissipation-inducing terms play a key 
role in the global existence theorem of Section QTJ Finally, in Section 15.51 we will show that if the initial data 
satisfy the Gauss and Codazzi constraints (|3.2.2ap - ()3.2.2bj) . and if the wave coordinate condition Q^\t=o = is 
satisfied, then Q^, I^u, Id$ = 0, and R^ u = R^ u ; i.e., under these conditions, the solution to f j5. 1 .5a|) - (I5.1.5bj) 
is also a solution to the irrotational Euler-Einstein system (I3,3.10aj) . 

5.2. Summary of the modified irrotational Euler-Einstein system for the equation of state p = c 2 s p. 
For convenience, we summarize (with the help of (I3.3.9P ) the results of the previous section by listing the 
modified irrotational Euler-Einstein system: 

(5.2.1a) Roo + 2ljT - 6oj 2 - Ag 00 - 2o s (d t <$>) 2 - (^y)^ +1 5oo = 0, 

(5.2.1b) R 0j + 2co(3cjg oj - Tj) - Ag 0j - 2a s - (^tJ^ + V = 0, (j = 1, 2, 3), 

(5.2.1c) R jk - Ag jk - 2a s (d J <S>)(d k <S>) - = , (j, k = 1,2,3), 

(5.2.1d) [ag af3 - 2s{d a <S>){d p $)]d a d^ - 3uad t <f> + 2sr aA/3 (d a $)(d A $)(d /3 $) = 0. 



5.3. Construction of initial data for the modified system. In general, we may consider arbitrary data for 
the modified equations (|5.2.1ap - (|5.2.1dp . However, if the solution of the modified system is also to be a solution 
of the Einstein equations (j3.2. lap . then we cannot choose the data arbitrarily. In this section, we assume that 
we are given initial data (£, g, K, d<&, VP) for the irrotational Euler-Einstein equations (|3.3.10ap - (|3.3.10bp as 
described in Section [3.2.21 In particular, we assume that it satisfies the constraints (|3.2.2ap - (|3.2.2bp . We will 
use this data to construct initial data for the modified equations that lead to a solution (Ai,g, 9<f>) of both the 
modified system and the Einstein equations; recall that a solution solves both systems <^=^ = 0, where 
is defined in (|5.1.2bp . 

To supply data for the modified equations, we must specify along £ = {t = 0} the full spacetime metric 
components ^ v | f=0 , their (future-directed) normal derivatives dtg^ u \t=07 (/**> v = 0, 1, 2, 3), the spatial derivatives 
9$|i = o of the fluid potential, and the normal derivative dt 3>|t=o of the fluid potential. To satisfy the requirements 

• £ = {t = 0} 

• g is the first fundamental form of £ 

• dt is normal to £ 

• K is the second fundamental form of £ 

• a$| s = 8$ 

° dcf 

• Dfi&\Y, = ^? , = differentiation in the direction of the future-directed normal to £, 
we set (for j, k = 1, 2, 3) 
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(5-3.1) goo\t=o — — 1, 9oj\t=o — 0, 9jk\t=o — 9jk, 

(5.3.2) d$\ t=0 = 8$, d t <f>\ t = = if, d t g jk \ t=0 = 2K jk . 

Furthermore, we need to satisfy the wave coordinate condition Qu\t=Q = 0, (p = 0, 1, 2, 3). To meet this need, 
we first calculate that 

(5.3.3) r | t=0 = ~{d t g o)\t=o - 9 ab K ab , 

(5.3.4) r,-| t=0 = -d t g 0j \ t=0 + X -g ab (2d a g bj - dfg ab ), (j = 1, 2, 3). 

Using (|5.3.3p and ()5.3.4p . the condition Q^\t=o = is easily seen to be equivalent to the following relations, 
where oj is defined in f|4. 1.3|) : 

(5.3.5) d t g 00 \ t=0 = 2(-3w| fc=0 g o\t=o -g ab K ab ) = 2(3w(0) - g ab K ab ), 

-l 

(5.3.6) d t g 0j \ t=0 = -3uj\ t =o goj\t=o +^g ab (2d a g bj - djg ab ) = g ab (d a g bj - ^djg ab ^ , (j = 1, 2, 3). 

o 

This completes our specification of the data for the modified equations. 

5.4. Decomposition of the modified irrotational Euler-Einstein system in wave coordinates. Nat- 
urally, the key to our global existence theorem involves a careful analysis of the nonlinear terms. In order to 
better see their structure, we dedicate this section to a decomposition of the modified system (15.2. lap - (15.2. lcl|) 
into principal terms and error terms, which we denote by variations of the symbol A. The estimates of Section 
[9]will justify the claim that the A terms are in fact error terms. We begin by recalling the previously mentioned 
rescaling hj k of the spatial indices of the metric: 

(5.4.1) h jk = e- 2Q g jk , (j, k = 1, 2, 3). 

The decomposition is captured in the next proposition. 

Proposition 5.4.1. (Decomposition of the modified equations) The equations (|5.2.1ap - (j5.2.1c[) can be 

written as follows: 

(5.4.2a) Dgigoo + 1) = 5Hd t g 00 + 6# 2 (<? o + 1) + A 00 , 

(5.4.2b) n g90j = md t g 0j + 2H 2 g 0j - 2Hg ab T ajb + A 0j , (j = 1, 2, 3), 

(5.4.2c) ti g h jk = 3Hd t h jk + A jk , (j,k = 1,2,3), 

where H = A/3, and the error terms A^ u are defined in (I5.4.4ip - (I5.4.4kp below. 
Furthermore, the fluid equation (|5.2.1dp can be written as 



(5.4.3) 



□ m $ = xujd t <$> + A a $, 
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where uj(t) is defined in (|4.1.3p , mP- v (u,u = 0,1,2, 3) is the reciprocal acoustical metric, 

□m = —d 2 + m ab d a d b + 2m 0a dtd a is the reduced wave operator corresponding to m^ u , the error term Aq$ is 
defined below in (|5.4.4hD . 

(5.4.4b) m 00 = -1, 

9 : 

(1 + 2s) + A (m) ' 



(5.4.4c) r^ fc - — ™ 



(5.4.4d) m 
(5.4.4e) 



g 00 g°J(l + 2s) + 2e n [(s + l)g°V Q + sg°°g a i] z a + e 2 "( g y ft + 2sg^g 0b )z a z b 



(l + 2s) + A (m) 



(5.4.4f) 

A (m) = (1 + 2s)(<? 00 + l)(g 00 - 1) + 2(1 + 2s)e n g 00 g 0a z a + e 2n (g 00 g ab + 2sg 0a g 0b )z a z b , 
(5.4.4g) 

A M = (5°° + 1)9* + 2s 5 °^ 0fc + 2e VV a + 2 + e 2 W" + 2sg^ g bk )z a z b , 

(5.4.4h) 



= Loi 7T"T~9 \ _i_ A 

I (1 + 2s) + A (m) 



1 



3a;( 5 00 + l)dt$ + 6weV a ^$ + (3 - 2s)we 2 V 6 *a*b3t$ 



(1 + 2s) + A (m) 

+ 2s(a°°^$ + [A$> + AgJ + A°^]d a <D + e 2Q [A^ + A$ + A^]^5t* + e m A^z^) |, 
i/ie A^pj (//, a, z/ = 0, 1, 2, 3j are defined below in Lemma \5.4-5\ 

(5.4.41) ^Aqo = A A00 + Ac,oo + [f(d$) - - (g 00 + l)/(0$) - ^-(<7oo + 

+ l{u - H)8 t g m + 3(^ 2 - H 2 )(g 00 + 1), 
(54.4j) ^A 0j = A Aoj + A c , 0j + *~ 1 a s+1 g oj - 2a s (9^)(5,cl>) - -f-^+i^. 

+ !( w _ tf)^ 5oj + (u 2 - H 2 )g 0j - (w - H)g ah T ajb , 
(5.4.4k) iA ifc = e" 2fi A Aifc - ( 5 00 + l)a s+1 h jk - 2u;g 0a d a h jk - 2e" 2 V (d^)(d k ^) 



-a s+1 )h jk + l(u-H)d t h jk , 



(5.4.41) f{d$) = f 2a s (d t $) 2 — a s+ \ 

s + 1 

the A A ^ u are defined in (|5.4.10ap - (|5.4.10cp . and Aq^q, A(j,oj are defined in (|5.4.12ap - (|5.4.12bp . In the above 
expressions, quantities associated to the background solution of Section^ are decorated with the symbol ~ . 

Remark 5.4.1. Note that there are two metrics in this system of PDEs: the inverse spacetime metric g^ u , and 
the reciprocal acoustical metric m^" ' . 
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Remark 5.4.2. As discussed in Section \1.2\ the variable Zj has been introduced to facilitate our analysis of the 
ratio of the size of the spatial derivatives of <1> to its time derivative. 

Proof. The proof involves a series of tedious computations, some of which are contained in Lemmas 15.4.21 - [5.4.51 
below. We sketch the proof of (|5.4.2c|) and leave the remaining equations to the reader. To obtain ()5.4.2cp . we 
first use equation (j5.2. lc|) . Lemma [5.4.21 and Lemma 15.4.31 to obtain the following equation for hj k = e~ 2Q gj k : 



(5.4.5) n g h jk = 3ood t h jk - 2{g m + l){d t u)h jk - 4ujg 0a d a h jk + 2 [ e - 2n A AJk - 2e~ m a s {d^){d k ^)] 

+ 2{(3w 2 - A + d t u)h jk - j^Y<T s+1 h jk }. 

We now use (|4.1.11ap - (|4.1.11bp to substitute in equation f)5 .4. 5[) for 3oj 2 — A and dtui in terms of a, arriving at 
the following equation: 



(5.4.6) U g h jk = 3u>d t h jk + 2(g 00 + l)a s+1 h jk - 4wg 0a d a h jk + 2[e~ 2n A Adk - 2 e - 2 V(0 j $)(0 fc *)] 

+ 2^-{a s+1 -a s+l )h jk . 
s ~\~ ± 

Equation (I5.4.2c[) now easily follows from (15.4.6|) . We remark that the proofs of (I5.4.2ap and (15.4.2b|) require 
the use of Lemma I5.4.4} and that the proof of (15.4.3|) requires the use of Lemma 15.4.51 To obtain (|5.4.3p , it is 
also useful to note that — (dt&) 2 {(l + 2s) + A/ m )j is the coefficient of the differential operator d 2 in equation 
(15.2. ldp . □ 



We now state the following four lemmas, which are needed for the proof of Proposition 15.4.11 
Lemma 5.4.2. The modified Ricci tensor from (|5.1.3p can be decomposed as follows: 

13 3 

(5.4.7) R^ v = --U g g^ v + -{go^dyuj + govd^u) + -ujd t g^ v + A^, {fi,u = 0, 1,2,3), 

where 

(5.4.8) A^ = g a Pg* x l{d a g VK ){d p g^)-T avK T^ x ], (/x,i/ = 0, 1,2,3). 

Proof. These computations are carried out in Lemma 4 of [28]. □ 

Lemma 5.4.3. The term A^ u ([i,v = 0, 1,2,3) defined in ([5.4.8P can be decomposed into principal terms and 
error terms A A ^ V as follows: 



(5.4.9a) 


^oo 


= 3lo 2 - ug ab d t g ab + 2ug ab d a g 0b + A A00 , 






(5.4.9b) 




= 2u J g 00 d t g 0j - 2cj 2 g 0O g Oj - u}g w d j3ao + cog ab T ajb + A A , j, 


(J = 1, 


2,3), 


(5.4.9c) 


Aj k 


= 2ug 00 d t g jk - 2oj 2 g 00 g jk + A Ajk , 


U,k = 


1,2,3) 
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where 
(5.4.10a) 

Aaoo = (g 00 ) 2 {(d t goo) 2 - (r oo) 2 } + g 00 g 0a {2(dtgoo)(d t goa + d a g o) - 4r oor oa} 
+ g 00 g ab {(d t goa)(d t g ob ) + (d a goo)(d b g o) - 2roo o r o&} 
+ g 0a g 0b {2(d t g 00 ){d a g 0b ) + 2(d t g 0b )(d a g 00 ) - 2T 000 T a0b - 2r o 6 r oa} 
+ g ab g 0l {2{d t g 0a )(d mb ) + 2{d b g m )(d a g m ) - 4r 00a r /0b } 
+ g ab g lm (d a g i)(d b go m ) + \g lm { g ab d t g a i - 2ub\ )(d b g 0m + d m9ob ) 

e 2n g ab dthal _2ugMg m 

-]g ab g lm (d a goi+d ig0a )(d b g 0m + d m g ob )-]( g ab d t9al - 2u5 b )( g lm d t g bm - 2ub\ ), 
4 4 ^ v ' v v ' 

e 2n g a,b dthal _2ugOb gol e 2n g lm dthbm _2ojgOlg ob 

(5.4.10b) 

A A oj = (5°°) 2 {(5^oo)(^5oi) - roooFoio} 

+ g 00 g 0a [(dtgoo)(d t gaj + d a g 0j ) + (d t goj)(d t goa + d a g 00 ) - 2T ooT ja - 2r j rooa} 
+ g 00 (g ab d t g bj - 2^ ) (d t g 0a - ld a9oo ) + \g m g ab {d a g m ){d b g 03 + d jgob ) 
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e 2n g ab dthbj-9 0a 90j 
„0a 06 



2 JU 7 2 

_l_ g Ua g w Udtgoo)(d a g bj ) + (dtg ob )(d a g 0j ) + {d a g m ){d t g bj ) + (d a g ob )(d t goj) 

— rooor a jfe — 2roo6roj a — r a obrojo| 

+ g ab g 0l [(dtgoa)(dig bj ) + (dig 0a )(dtg b j) + (d b g o)(d a gij) + (d b goi){d a g j) - 2T 00a Ti jb } 

- g ab g° l [(dig a + d a g i)T 0jb - ^(d t gi a ){d b g j - <% fc)j 

+ ojg 0a (d t g aj - 2ojg aj ) + \g m { g ab d t g la - 2uSf )d t g bj 
s „ ' 2 » v ' 

e 2n dth a3 e™g" b dth la -2wg 0b g 0l 

+ g ab g lm {(d a g i)(d b g mj ) - ^(d a g ol + d ig0 a)r bjm } + \g ab ( g lm d t gia -2u;5™) T bjm , 

e 2n g l mdthla _ g 0,n g0a 

(5.4.10c) 

&A,jk = (g 00 ) 2 {(dtgoj)(dtgok) -r 0j0 r oko } 

+ g 00 g 0a {(dtgoj)(dtg a k + d a90 k) + (dtgok)(d t g a j + <9 a 5f j) - 2r j r fca - 2r fcor ja} 

+ g°°g ab {(d a goj)(d b gok) - ^(d a g j - djgo a ){d b g ok - d k g ob )} 

-lg 00 {( g ab d t g aj -2uj5 b )(d b g ok - d k9ob ) + ( g ab d t g bk - 2u5 a k )(d a9oj - % 0a )| 
2 v v > „ ' > 

e m g ab dtha ._ 2 Lug» b g 0] e 2n g« b d t h bk -2u)g 0a g 0k 

+ u J g 00 ( g bk g ab ^-5 a k ) d t9aj + ±g 00 ( g ab d t g a] - 2a^ ){ d t g bk -2ug bk ) 

-gok9 0a e 2n g ab d t h a] -2ujg° b g 0j e 2fl dth bk 

+ g 0a g oh {(dtgoj)(d a gbk) + (d t g b j)(d a gok) + (d a goj)(d t g b k) + (d a g b j)(d t gok) 

— ^OjO^akh — 2T()j b To ka — T a j b T() k o^ 

, aft 01 ( 
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Proof. These computations are carried out in Lemma 5 of [28]. □ 



Lemma 5.4.4. The sums Aqq, + ioo an d ^Oj + loji (j = 1)2,3), can be decomposed into principal terms and 
error terms as follows, where Ioo,-^Oj are defined in (|5.1.6ap - (j5.1.6b|) ; Aqo,Aoj are defined in ([5.4. 80 : and 
Aa,oo, A A ,oi are defined in ()5.4.10a|) - (|5.4.10b|) : 



(5.4.11a) A 00 + 2wr° - 6oj 2 = cjd t g 00 + 3cj 2 (g 00 + 1) + 3w 2 #oo + A Aj0 o + A Ci00 , 

(5.4.11b) A 0j + 2u(3ujg 0j - Tj) = 4io 2 g 0j - ujg ab T ajb + A Afij + A c , 0j , (j = 1, 2, 3), 



where 



(5.4.12a) 

Ac.oo = -6( 5 oo)-^ 2 {(5oo + I) 2 - 9° a 90a} - u(g°° + 1)( f%9ab ~ 6a; ) + 2u,( 5 00 + l) 5 ab a a <7o b 

e^g^d t h ab ~2ujg^g 0a 

+ a;( 5 00 + - l)d t9m + 2^ 00 5 0a (r 0a o + 2r 00a ) + ^g 0a g ob T 0ab + 2^°^, 

(5.4.12b) 

A c , 0j = 2uj 2 (g 00 + l)g 0j - 2ujg 0a \( d t g aj - 2ug aj ) + d a g 0j - d j90a }, (j = 1, 2, 3). 



Proof. These computations are carried out in Lemma 6 of [25]. □ 



Lemma 5.4.5. The fully raised Christoffel symbols Y^ au (\i,ol,v = 0,1,2,3) can be decomposed into principal 
terms and error terms A^T as follows: 



(5.4.13a) 


p000 


a 000 
- A (r) , 






(5.4.13b) 


poo 


_ F 00i _ ajOO 

-^(r)) 


(i = 


1,2,3), 


(5.4.13c) 




- A 0j0 

- ^ (r) , 


(i = 


1,2,3), 


(5.4.13d) 




= r fc J'° = -ugi* + A{£*, 


U,k 


= 1,2,3), 


(5.4.13e) 


pjOfc 




(j,k 


= 1,2,3), 


(5.4.13f) 


pijfc 


_ Aijk 
- ^ {r) , 


(hJ, 


k = 1,2,3) 



where 
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Aj$ = \{g™?d t9m + \g m 9 0a g \d t g ah + 2d a g 0b ) 



+ \(g 00 ) 2 g 0a (2d t g 0a + d a g 00 ) + \g 0a 9 0b g 0l d a g bl , 
(5.4.14b) A$° = ^( 5 00 ) V^soo + 5 ai ^5oo) + g 00 g oj g a0 d t g 0a + g 00 g aj g ob d a g ob + \g a3 g ob g ol d a g b i, 

(5.4.14c) Ajg = i(g 00 ) V^Soo + - <? OJ 'd <?oo) 

+ <?°°(<? a V^ afe + rt^oo + <? a V^a<7o& - g a 'g ob d a9ob ) 
+ g 0a g bj g ol d a g bl - \g Qa g bj g 01 d b9ah 

(5.4.14d) A°f* = ^°°{<?V fe %00 + g° j 9 ak d a900 + g aj g 0k (2d t g 0a - daSoo) + sV^SOa - daSoft)} 

+ ^{sV^aSW + g° j g bk (d a gob + 5 b50 a - d t9ab ) 

+ g bj g° k (d t g ab + d a g ob - d b g 0a ) + g bj g lk (d a g bl + d l9ab - d b g al )} 

+ \g m { g aj g bk d t9ab -2ujg^ ) + u(g 00 + l)g^ k , 

g bk (e 2n 9 aJ dth ab -2u,g°ig 0b ) 

(5.4.14e) A$ k = l -g m {g^g 0k d t9m + g 0j g ak d a9oo + g aj g ok d a g 00 + g aj g bk (d a g ob + d b g 0a )} 
+ lg 0a (g 0j g 0k (2d t g a - d a g o) + g 0j g bk (d t 9ab + d&5oa - d a g ob ) 



+ g bj g ok (d t gab + d b g 0a - d a g ob ) + g bj g lk (d b g a i + dig ab - d a gu)} 



1 oo ( ffi^dtgrt-tog* ) - u(g 00 + l)g^ k , 

2 - v ' 



gbk( e 2Sl g aj dthab _2ug°ig 0b ) 

(5.4.14f) A$ = ^VV&Soo + ^(<7 a W* + <7°W* - 0°W*) d a9oo + g 0i g aj g 0k d t g Oa 

+ \ {g ai g bj g oh + g m g 0j g bk + g H g oj g ak + g 0l g bj g ak - g M g aj g ok - g 0i g aj g bk ) d a g ob 

+ \ {g 0l g aj g bh + g ai g bj g ok - g at g 0j g bk ) d t g ab + \ (g m g b3 g lk + g u d 3 g ak - g bl g ai g lk ) d a9bl . 

Proof. The proof is again a series of tedious computations that follow from the formula 

TU .au = y^g^ g ^(d K g X(7 + d a g K \ - d x g KtJ ). □ 

5.5. Classical local existence. In this section, we discuss classical local existence results for the modified 
system of PDEs (|5.4.2ap - (|5.4.3p . The theorems in this section are stated without proof; they can be proved 
using standard methods described in e.g. |19[ Ch. VI], |34[ Ch. 16]. Also see \2&\ Proposition 1]. 

Theorem 5.1. (Local existence) Let N > 3 be an integer. Let g^ u = gu,v\t=o, 2K^ U = (Ot^^) |t=o, (f-i, v = 
0,1,2,3), d<& = 0$|t=o, and ^> = (dt&)\t=o be initial data (not necessarily satisfying the Einstein constraints) 
on the manifold £ = T 3 for the modified system (|5.4.2a[) - (15.4.31) satisfying (for j, k = 1, 2, 3) 

(5.5.1a) gao + l€H N+1 , goj G H N+1 , Bg jk G H N , 

(5.5.1b) K 00 eH N , K 0j eH N , K jk -uj(0)e 2n ^g jk £H N , 

(5.5.1c) <9<l> G H N , 9 - # G H N , 
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where ^ > is a constant. Assume further that there are constants C\,C 2 > 1 and C3 > such that 

(5.5.2) C^5 ab X a X b < g ab X a X b < C 1 5 ab X a X b , V(X\X 2 ,X 3 ) G R 3 , 

(5.5.3) C^d ab X a X b < m ab X a X b < C 2 5 ab X a X b , y(X 1 ,X 2 ,X 3 ) G M 3 , 

(5.5.4) goo < "C 3 , 

where the m? k , (j,k = 1,2,3) are the spatial components of the reciprocal acoustical metric (see (|5.4.4cp ). Then 
these data launch a unique classical solution [g^ v ,d^), (/U, v = 0, 1,2,3), to the modified system existing on an 
interval (T_,T+), with Tl < < T+, such that 

(5.5.5) 9flu G C»- 1 ((T-,T+) x T 3 ), d u $ G Cf" 2 ((T_ , T+) x T 3 ), 

such that goo < 0, and such that the eigenvalues of the 3x3 matrices g^ and m? k are uniformly bounded below 
from and from above. 

The solution has the following regularity properties: 

(5.5.6a) 

<7oo + l G C°((T.,T + ),H N+1 ), g 0] G C°((T_, T + ), J^ +1 ), G C°((T^,T + ), H N ), 

(5.5.6b) 

3 t <?oo G C%(T-,T+),H N ), d t <? ; G C°((T-,T+),H N ), d t g jk - 2u(t)g jk G C°((T_, T+), 

(5.5.6c) 9$ G C°((T-,T + ),H N ), e^dt® - * G C°((r_, T+), H N ). 

Furthermore, g^ is a smooth Lorentz metric on (T_,T + ) x T 3 , and the sets {t} x T 3 are Cauchy hypersurfaces 

in the Lorentzian manifold (M ="^(T1,T + ) x T 3 ,^) fort G (T_,T+). 

In addition, there exists an open neighborhood O of {g^ y ,K^ u ,d^,^!) such that all data belonging to O 
launch solutions that also exist on the interval (T_,T + ) and that have the same regularity properties as (g,d&). 
Furthermore, on O, the map from the initial data to the solution is continuousW\ 

Finally, if, as described in Section 15.31 the data for the modified system are constructed from data for the 
irrotational Einstein- Euler system satisfying the constraints (|3,2,2ap - (|3.2.2bp on an open subset S C T 3 , and 
if the wave coordinate condition Q^\s = holds, then (g^^d^) is also a solution to the un-modified equations 
(|3.2.1ap - (|3.2.1bp onD(S), the Cauchy development of S. 

Remark 5.5.1. The hypotheses in Theorem 15.11 have been stated in a manner that allows us to apply to it 
initial data near that of the background solution of Section dJ Furthermore, we remark that the assumptions 
and conclusions concerning the metric components would appear more natural if expressed in terms of the 

variables hjk = f e~ 2n gjk', these rescaled quantities are the ones that we use in our global existence proof. 

Remark 5.5.2. The fact that {g^ p ,d^) is also a solution to the modified equations if the constraints and the 
wave coordinate condition |s = are satisfied is discussed more fully in Section 15.61 

In our proof of Theorem lll.il we will use the following continuation principle, which provides criteria that are 
sufficient to ensure that a solution to the modified equations exists globally in time. See |19[ Theorem 6.4.11] 
for the ideas behind a proof. 

Theorem 5.2. ( Continuation principle) Let T max be the supremum over all times T + such that the solution 
{g^ v ,d^), (u, v = 0,1,2,3), exists on the interval [0,7+) and has the properties stated in the conclusions of 
Theorem I5.il Let m^ v be the reciprocal acoustical metric, which is defined above in Proposition \5.4- 1\ Then if 
Tmax < 00 j one °f the following four possibilities must occur: 

(1) There is a sequence (t n ,x n ) G [0,T max ) x T 3 such that lim n _ 5 . 00 goo(t n , x n ) = 0. 



By continuous, we mean continuous relative to the norms on the data and the norms on the solution that are stated in the 
hypotheses and above conclusions of the theorem. 



Stability of Irrotational Euler-Einstein 

27 

(2) There is a sequence (t n ,x n ) £ [0, T max ) x T 3 such that the smallest eigenvalue of the 3x3 matrix 
9jk(tn, x n) converges to as n — > oo. 

(3) There is a sequence (t n ,x n ) G [0, T max ) x T 3 such that the smallest eigenvalue of the 3x3 matrix 
m^ k (t n ,x n ) converges to as n — > oo. 

(4) \im t -+ Tmax sup < T < t | J2l,u=o I Wlc* + £)U=o ll^lc, 1 + Ej=o \\ d ^h°° | = oo. 
Similar results hold for an interval of the form (T m i n , 0] . 

Remark 5.5.3. If one the first two possibilities occurs, then the hyperbolicity of equations (I5.4.2ap - (j5.4.2cj) 
breaks down. Similarly, if the third possibility occurs, then the hyperbolicity of equation (15, 4. 3ft breaks down. 

5.6. Preservation of the wave coordinate condition. In Section [5~3l from given initial data for the Einstein 
equations, we constructed initial data for the modified equations that in particular satisfy the wave coordinate 
condition along the Cauchy hypersurface S; i.e., Q^\t=o = 0. As mentioned in the statement of Theorem 15.11 
these data launch a solution of both the modified equations and the Einstein equations. As mentioned in Section 
[5j this fact follows from the fact that = in T>(T,). In the next proposition, we prove this fact. 

Proposition 5.6.1. (Preservation of wave coordinates) Let (S = {x G M. \ t = 0}, g~jk, Kjki d<&, , 
(j,k = 1,2,3), be initial data for the Einstein equations (|3.2.1ap - (|3.2.1bp that satisfy the constraints (|3.2.2ap 
- (I3.2.2bp . Let (g^u\t=o,dtg^ u \t=o^d^^)j (/•*) v = 0,1,2,3), be initial data for the modified equations (j5.1.5aj) - 
(I5.1.5bp that are constructed from the data for the Einstein equations as described in Section \5. 3[ In particular, 
we recall that the construction of Section [5.31 leads the fact that Q^\t=o = where is defined in (|5.1.2bp . 
Let (A4 , g^ v , d^) be the maximal globally hyperbolic development of the data for the modified equations. Then 
Q M = in £>(£). 

Remark 5.6.1. Although we assume that the equation of state satisfies the assumptions of Section O we do 
not assume here that it is of the form p = c 2 s p. 

Proof To prove Proposition 15.6.11 we first compute that for a solution of the modified equations, it follows that 
is a solution to the hyperbolic system 

(5.6.1) g a PD a DpQ„ + R«Q a + 2g^D a I^ - g^D^ = -4/ M D^, = 0, 1, 2, 3). 

where and Ig$, which depend linearly on the Q^, are defined in (I5.1.6ap - (I5.1.6dp . More specifically, 
to obtain equation (I5.6.ip . apply D u to each side of equation (I5.6.2P below, and use the Bianchi identity 
D u {Rilv — \Rdiiv) = 0, equation (I3.1.15p . and the curvature relation D fJi D a Q a = D a D^Q a — R^Q a . 

Since (|5.6.ip is a system of wave equations and is of hyperbolic character, the fact that = in T>(T,) would 
follow from a standard uniqueness theorem for such systems (see e.g. [19} Ch. VI], |34[ Ch.16]), together with 
the knowledge that both Q/Js = and dtQ^\s = hold. However, in constructing the data for the modified 
equations, we have already exhausted our gauge freedom. Although the construction of Section [5.31 has led to 
the condition Q/Je = 0, it seems that we have no way to enforce the condition dtQ^T, = 0. The remarkable fact, 
first exploited by Choquet-Bruhat in jS], is that the assumption that the original data for the Einstein equations 
satisfies the constraints (|3.2.2ap - (|3.2.2bp . together with the assumption = 0, automatically imply that 

dtQ^T. = 0. The remainder of the proof is dedicated to proving this fact. 

First, using definition (|5.1.3|) . we compute that for a solution of the modified equation (|5.1.5ap . the following 
identity holds: 



(5.6.2) 

R»v - \R9„u + Ag^ - T^ ca/ar) = ~\ {D»Qu + D U Q^) + \(D a Q a ) gi , u - V + ^I*^, (ji, " = 0, 1, 2, 3). 

The left-hand side of (|5.6.2p is simply the difference of the left and right sides of the Einstein equations 
(I3.2.1a|) . Since the initial data (£, g, K, (9$, 4/) for the Einstein equations are assumed to satisfy the constraints 
(I3.2.2a|) - (I3.2.2bp . as described in Section [3.2.21 it follows that the Einstein equations are satisfied along E; i.e., 
the left-hand side of (I5.6.2P vanishes at t = 0. Furthermore, since <3/Jj=o = 0) it follows from definitions (I5.1.6ap 



Stability of Irrotational Euler-Einstein 

28 

- (|5.1.6dp that I^, u \t=o = 0, (//, v = 0, 1,2,3), and Id$\t=o = 0. Using these facts and (|5.6,2p . we conclude that 
the following equation holds: 

(5.6.3) -\{D^Q u + D u Q fl )\ t=0 + ^(D a Q a )\ t=0 g^\ t=0 = 0, (ji, v = 0, 1, 2, 3). 

Now let iV^ = 8q denote the future-directed unit normal to S, and let = <5J, a vectorfield tangent to X. 

Contracting (|5.6.3|) with N^X U , setting i = 0, using the facts that Q^|t=o = and (X^d u Q ^^=0 = 0, we have 
that 

(5.6.4) = -^N^X^^dnQv + d u Q^)\ t =o = dtQj\t=o, (j = 1,2,3). 

We also contract (15-6.31) with N^N U set t = 0, and use (15. 6. 4ft . the fact that Q^\t=o = 0, and the fact that 
g 00 \t=o = -1, obtaining 

(5.6.5) d t Q \t=o = 0. 

From (|5.6.4p and (|5.6.5p . we conclude that the data for the system (|5.6.ip are trivial. 

□ 

6. Norms and Energies 

In this section, we define the Sobolev normiQ and energies that will play a central role in our global existence 
theorem of Section [TTJ Let us make a few comments on them. First, we remark that in Section [TOj, we will 
show that if the norms are sufficiently small, then they are equivalent to the energies; i.e., the energies can be 
used to control Sobolev norms of solutions. Next, we recall that the background solution variable <9<f> satisfies 
dt& = , §<& = 0, where ^> > is the constant defined in (14.1.8p , The quantity Sg^-N, which is introduced 

below in (|6.1.2ep . measures the difference of the perturbed variable ((?(<!>, d&) from the background (dt$, 0). We 
also follow Ringstrom by introducing scalings by e a ^, where a is a number, in the definitions of the norms and 
energies. The effect of these scalings is that in our proof of global existence, a convenient and viable bootstrap 
assumption to make for these quantities is that they are of size e, where e is sufficiently small. Finally, we 
remark that the positive number q that appears in this section and throughout this article is defined in (I8.2.2P 

below, and we remind the reader that hj k = e~ 2n gjk, (j, k = 1,2,3). 

6.1. Norms for g and d&. In this section, we introduce the weighted Sobolev norms that will be used in 
Section [9] to estimate the terms appearing in the modified equations. The weights are designed in order to make 
the bootstrap argument of Section [11] easy to close. 

Definition 6.1.1. We define the norms S go0+1 . N (t), Sg 0m . N (t), S h ^. N {t), S g . N {t), Sg^. N (t), S N (t), S goo+1 - N (t), 
Sgo*;N(t), S htt . N (t), S g]N (t), and S d $ ;N (t) as follows: 

(6.1.1a) S m+1 . N = e^\\d t g 00 \\ H N + e^\\g o + + ^"^H^ooH^, 

3 

(6.1.1b) ^„^;7sr 5^ (e^- 1 )"!!^^^!!^- -K e<^- 1 >^|| i 7o J -||^^ -H e 2 >^ 1 1 <9^ 0j - 1 1 ) , 

3=1 
3 

(6.1.1c) S htt . N = £ (e« Q \\d t h jk \\ H N + \\8h jk \\ H N-i +e^ n \\Bh jk \\ HN ), 

j,k=l w 

absent if ^=0 

(6.1. Id) S^. N d ^ \\e* u d t $ - n HN + e^ a \\8^\\ Hlfi 

(6.1. le) 



'Technically, Sq^-.n is a norm of the difference of <9$ and the derivatives 9$ of the background potential. 
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(6.1.2a) S goo+1 . N (t) = sup S +1n (t), 

0<T<t 

(6.1.2b) S gat - N (t) = sup S n (t), 

0<T<t 

(6.1.2c) S htt . N {t) = sup S htt . N (r), 

0<T<t 

(6.1. 2d) Sg-N = S go0+ i-N + S g0t -N + Sfi**;N , 

(6.1. 2e) Sdn>;N(t) = sup Sg^. N (r), 

0<T<t 

(6.1.2f) Sat = f Sg-N + Sg^.pf. 

6.2. Energies for the metric g. 

6.2.1. The building block energy for g. The energies for the metric components will be built from the quantities 
defined in the following lemma. They are designed with equations (|5.4.2ap - (|5.4.2c|) in mind. 

Lemma 6.2.1. Let v be a solution to the scalar equation 

(6.2.1) a g v = aHd t v + PH 2 v + F, 

where n g = g XK d\d K , a > and /3 > 0. For any constants 7 > 0, 5 > 0, we define 

(6.2.2) £^ S) [v, dv] = f X - J^{-g™(d t vf + g ab {d a v){d b v) - 2 1 Hg 00 vd t v + 5H 2 v 2 } d 3 x. 

Then there exist constants i] > 0, C > 0, 5 > 0, and 7 > 0, with 7] and C depending on a, f3, 7 and S, such that 

(6.2.3) \g 00 + M<v 

implies that 

(6.2.4) £l S )[v,dv] >C f (d t v) 2 + g ab (d a v)(d b v) + C h) v 2 d 3 x, 

where C( 7 ) =0 1/7 = and C( 7 ) = 1 if j > 0. Furthermore, if (3 = 0, i/ien 7 = 5 = 0. Finally, we have that 

(6.2.5) j t (£f J>6) [v, dv]) < -r]HS 2 lS) [v, dv] + J { - + 7 H«)F + A £ .^ S) [v, dv]} d 3 x, 
where 

(6.2.6) A £ . M [v,dv] = -"/H(d a g ab )vd b v - 2 1 H(d a g 0a )vd t v - 2 1 Hg* a {d a v)(d t v) 

- (d a g 0a )(d t v) 2 - (d a g ab )(d b v)(d t v) - ^(d t g 00 )(d t v) 2 

+ (±d t g ab + ug ab ^j (d a v)(d b v) + {H- u)g ab (d a v)(d b v) 

- 7 H(d t g 00 )vd t v - 7^(5°° + l)(d t v) 2 . 

Proof. A proof based on a standard integration by parts argument (multiply both sides of equation (|6.2.ip by 
— (dtv + jHv) before integrating over T 3 )) is given in Lemma 15 of |28| . In particular, we quote the following 
identity: 



30 



Stability of Irrotational Euler-Einstein 



(6.2.7) -^( £ ( 7 ,5)h dv}) = j T3 {-(<*- 7)#(<9^) 2 + (6-0- ja)H 2 vd t v - f3 7 H 3 v 2 

- (1 + 7 )Hg ab (d a v)(d b v) - (d t v + 1 Hv)F + A £ . { ^ s) [v, 3v]} d 3 x. 

□ 

6.2.2. Energies for the components of g. In this section, we will use re-scaled versions of energies of the form 
(|6.2.2p to construct energies for the components of g. 

Definition 6.2.1. We define the non-negative energies E goo+1 . N (t), E goo+1;N (t), E gQt . N (t), E g0m , N (t), E htt . N (t), 
Eh,,;N(t), E_ g . N (t), and Eg.tf(t) as follows: 

(6.2.8a) E_l Q+1]N d =< £ e 2qQ £? JOO , Soo) [d S (9oo + l),d(d s9oo )}, 

\5\<N 

(6.2.8b) E^. N d = f E e2(9 " 1)Q 4o.,^)[^Oi, d(d s90j )], 

\5\<Nj=l 

(6.2.8c) ^L;JV= E { E e 2 " n £f m [0,d(d d h, k )} + ^ j C(} H 2 (d d h jk {r)f d 3 x}, 

\d\<N j,k=l ^ T3 
(6.2.8d) E 2 g;N = K 2 goo+1 . N + E 2 ( ^. N + E_l tt . N , 

(6.2.9a) E goo+1 . N (t) = sup E +in (t), 

0<T<t 

(6.2.9b) E go ,. N (t) = sup E n (t), 

0<T<t 

(6.2.9c) ^« ; jv(t) d =^„ ;j v(r), 

(6-2. 9d) E 2 . N A = E goo+1 . N + E 2 Qt . N + El mm . N , 

where 

(6.2.10a) h jk = e~ 2n g jk , (j, k = 1,2,3), 

(6.2.10b) c a = 0, if|a| = 0, 

(6.2.10c) cs = l, if|a|>0, 

and (700, o~oo), (7o* 5 6o*)i an d (7**, = (0, 0) are the constants generated by applying Lemma [6.2,ll to equations 
(|5.4.2ap - (|5.4.2c|) respectively. 

In the next lemma, we provide a preliminary estimate of the time derivative of these energies. 

Lemma 6.2.2. Assume that (g^, d^Q), (fJ,,v = 0,1,2,3), is a solution to the modified equations ()5.4.2a|) - 
(j5.4.2cj) . and let E_ gm+1 . N , E_ g0t . N , and E_ htt . N be as in Definition \6.2.1[ Let [n g , 83] denote the commutator of 
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the operators U g and dg. Then the following differential inequalities are satisfied, where Ag.^^]/, d(-)] is defined 
in (|6.2.6p . and the constants (700, 5qq), (70*, $o*), and (7**, 5**) = (0,0) are defined in Definition \6.2.1\ 



(6.2.11a) 

J t (E 2 goo+1 , N ) < (2q - V0O )HE 2 goo+1 . N + 2q{u - H)E 2 gQQ+1]N 



\5\<N 

(6.2.11b) 
d_ 

dt 



[ e 2qU {dtd s (goo + 1) + looHd s (g 00 + l)}{d s A 00 + [n g , d d ](g 00 + l)}d 
S|<iV ^ T3 

V / e 2qQ A £ . i7QOtSoo) [d a (goo + l),d(d s goo)]d 3 x, 



tM^n) < [2(9 - 1) - Vo*]HE 2 g0t . N + 2( 9 - l)(w - ^)^ . ;JV 



3 

' E E /e 2(9 ~ 1)n {dt^+7o*^^ 

\5\<Nj=l Jj3 
3 . 

+ EE e 2 ^ 1 )°A £ . (70 , i ^ ) [^,a(a S 50i)]d 3 x, 
laKJV .7=1 ,/T3 



|a|<JV j 

(6.2.11c) 

d Ml^, N ) < (29 - V **)HEl^. N + 29( W - ff)^L.;JV 



dt 



3 

V V / e 2gn {dtd 3 h jk + 1** Hd d h jk }{d$A jk + [n g ,d s ]h jk } d 3 x 

/T3 v —v-~'' 

;|<Arj,fc=i JJ1 
3 

^ ^ / e 2 ^A £;(0i0) [0,9(9 cf V)]^+ E / H 2 (d 3 d t h jk )(d s h jk )d 3 x. 

*\<N j,k=l 3 l<\S\<N Jt3 



Proof. Lemma 16.2.21 follows easily from definitions (|6.2.8ap - (|6.2.9d|) . and from (|6.2.5|) . □ 



The following corollary follows easily from Lemma 16.2.21 definitions (|6.2.8ap - (|6.2.8cp . and the Cauchy- 
Schwarz inequality for integrals. 
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Corollary 6.2.3. Under the assumptions of Lemma 1 6. 2. 21 we have that 



Ml )0+1 . N ) < (2q - 7 ]00 )HE 2 gQO+1 . N + 2q(u - H)E 2 gQ0+1 . N + S goo+1 . N e« Q \\ 



6.2.12a) 

d_ 
~dt' 

+ S 900 +i;N E e 9U \\lu g ,d 3 }(9oo + l)\\ L *+ E e29Q ll A £;(7Qo,5oo)[^(5oo + 1), d(d sgoo )]\\ L i 

\a\<N \a\<N 

6.2.12b) 

d 



■£(ML;n) < [2(? " 1) " %*]#^o,;iV + 2 (9 -!)("- #)#^ 



+ 2^5 90t; ^^e^ 1 ^||^r a , 6 ||^+^^^e^ 1 ^||A ,||^ 
j=i j'=i 

3 3 

+ E E e(9_1)n |l[^'^boillL 2 + E E e%_1)n H A *;(Wo*) [%30i, lb , 

|cf|<JVi=l |a|<JVj'=l 

(6.2.12c) 

d 3 
j t (Kl^. N ) < (2q - i M )HE 2 htt . N + 2g(a; - H)E_l^. N + S h ^. N £ ||A jfc ||^ 

i,fc=i 

3 3 

+ 5^ ;jV E E ll[° fl ,W ifc || i3 + E E II A £;(o,o) [0, S(^- fe ) + i? 2 e-^_.^, 

\S\<Nj,k=l \a\<Nj,k=l 
where the norms S_ goo+1 . N , S_ gQt . N , S_ h ^^. N are defined in Definition \6.1.1[ 

6.3. Energies for the fluid variable d&. In this section, we define the energies that we will use to study the 
irrotational fluid equation (15. 4, 3ft . We begin by stating their definitions. 

Definition 6.3.1. Let ^ be the positive constant defined in (|4,1.8p . Then we define the following non-negative 
energies Eq§ N (t) , Eg^.^it) for <9<£ as follows: 

(6.3.1a) E 2 ^ d ^ \ [ (e^dt® - *) 2 + e 2M m ab (d a $>){d h <S>) d 3 x, 

(6.3.1b) El^. N ^E 2 + E If e 2 ^(d t d^) 2 + e 2M m ab (d a d^)(d b d^)d\, (N>1), 

— 2 Ifs 

1<\3\<N J 11 

(6.3.1c) E d ^. N (t) = sup Eq^. n (t). 

0<T<t 

In the next lemma, we provide a preliminary estimate of the time derivative of the fluid energies. 
Lemma 6.3.1. Let d& be a solution to 

(6.3.2) Dm$ = >cu{t)d t § + A 9 $, 
where 

(6.3.3) D m = J -«9 2 + m afe d a d b + 2m 0a d t d a 

is the reduced wave operator corresponding to the reciprocal acoustical metric m^ u , (fJ,,u = 0, 1,2,3), and 
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Let [n m ,dci] denote the commutator of the operators n m and d$. Then 
(6.3.5) 



i(M*-,N) = -f (d a m° a )(.e* n d t <t>-y) 2 d 3 x - V f e 2M (d a m 0a )(d t d d ^) 2 d 3 x 

- I e M {d a m ab ){e M d t <5> - ^>){d b <$>) d 3 x - V / e 2xQ (a 6 m a6 )(^a 5 $)(a 6 a s $) d 3 x 
Jt3 i<\s\<n JtZ 

- I e* n (A d<s> )(e xn d t <5>-^)d 3 x- V f e 2>tQ {d s A 8 ^)d t d^ d 3 x 

- Y, f e 2M {[n m ,d d \$)d t d ci <5>d 3 x + \ Y, I e 2M (d t m ab + 2xu:m ab )(d a d^)(d b d^)d 
l<\a\<N Jj3 \a\<N JTA 

Proof. This is a standard integration by parts lemma that can be proved using the ideas of Lemma 16.2.11 We 
provide a sketch of the proof. We begin by differentiating under the integral in the definition of Eg$. N to 
conclude that 



3 x. 



■3.6) ±(^. N )= V ( {d d {e M d^-^}d t d d {e M d^) + e 2M m ab {d t d a d^){d b d^)d 3 x 

+ Y I \e 2 * n (^d t m ab + ^ojm ab yd a d^)(d b d^)d 3 x. 



\d\<N' 

For each fixed a, we will now eliminate the highest derivatives of $ in f|6.3.6j) (i.e., the derivatives of order- 
Id?! +2). To this end, we first differentiate equation (|6.3.2p with d d and multiply both sides of the equation by 
e x ^, which allows us to express the resulting equality as 

(6.3.7) -dtdsie^dt* - *) + e* U m ab d a d b d s $ = e M d d A d „ + e M [u m , d s ]$ - 2e* n m 0a d t d a d s $ . 

We then multiply both sides of (|6.3.7p by —d$(e dt& — ^), integrate over T 3 , and integrate by parts. Inserting 
the identity that arises into (|6.3,6p . we arrive at (|6.3.5p . □ 

We now state the following corollary, which follows easily from definition (|6.1.1dp . Lemma 16.3.11 and the 
Cauchy-Schwarz inequality for integrals. 

Corollary 6.3.2. Under the hypotheses of Lemma 1 6. 3. 1\ we have that 

d 3 

(6.3.8) T t ^;N) < S%*;N\\dam 0a \\ L °° +^. JV ^e n ||5 m a6 || i oc + S S9 . N e* l \\A a *\\ H * 

6=1 

1 3 

+ iW E e ^ll[flm,9 5 ]$|| L2 + -51$.^ Y e 2xn \\d t m ab + 2um ab \\ L ^ 

l<\a\<N a,b=l 

+ (x-1)ljY I e 2M m ab (d a d^)(d b d^)d 3 x, 

where S_q$-n * s defined in ()6.1.1dD . 
6.4. The total energy E^y. 

Definition 6.4.1. Let E g .^ and Eq$-n be the metric component and fluid energies defined in (I6.2.9dp and 
(16.3. lcp respectively. We define Ejv, the total energy associated to (g^t,, c^<3?), (/i, v = 0, 1, 2, 3), as follows: 



(6.4.1) 



Eat — Eg-N + Eqq-n. 
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7. Linear- Algebraic Estimates of g^ u and g^ v , (/j,, v = 0, 1,2, 3) 

In this section, we provide some linear- algebraic estimates of g^ u and g^ v '. In addition to providing some 
rough L°° estimates that we will use in Sections [9] and [TOl the lemmas that we prove will guarantee that g^ is 
a Lorentzian metric. We remark that we already made use of this fact in our statement of the conclusions of 
Theorem 15.11 The estimates are based on the following rough assumptions, which we will upgrade during our 
global existence argument. 

Rough Bootstrap Assumptions for g^ : 
In this section, we assume that there are constants g > and c\ > 1 such that 



(7.0.1a) 
(7.0.1b) 

(7.0.1c) 



boo + 1| < g, 
c^5 ab X a X b < e~ m g ab X a X b < Cl 5 ab X a X b , 



V(X\X 2 ,X 3 ) G M 3 , 



^2 \90a\ 2 < V c i le 



l p 2(l-q)0 



a=l 



For our global existence argument, we will assume that g = g m i n , where g m in is defined in Section I8.2i 

Lemma 7.0.1. Let g^ v be a symmetric 4x4 matrix of real numbers. Let (<7t,)jfc be the 3x3 matrix defined by 
(g\i)jk = 9jki let {9^ l y k be the 3x3 inverse of (g\,)jk- Assume that goo < and that {g\,)jk is positive definite. 

Then g^ u is a Lorentzian metric with inverse g^", g 00 < 0, and the 3x3 matrix (g^y k defined by (g^Y k d = gi k 
is positive definite. Furthermore, the following formulas hold: 



(7.0.2a) 
(7.0.2b) 

(7.0.2c) 

where 



,00 



900 



d 2 ' 



9oo 



9 -^{g^f b X a X b < (g#) ab X a X b < (g^XaX,, 



-W— (#b 1 ) aJ 90a, 

d z - goo 



i°3 



V(Ii,I 2 ,I 3 )eK 3 , 
(j = 1,2,3), 



-l\ab. 



(7.0.3) d = {g^ ^rgoagob- 

Proof. Lemma 17.0.11 is a combining of Lemmas 1 and 2 of [28] . 



□ 



Lemma 7.0.2. Let g^ v be a symmetric 4x4 matrix of real numbers satisfying (|7.0.1ap - (|7.0. lcj) . where S7 > 
and < q < 1. Then g is a Lorentzian metric, and there exists a constant go > such that < g < go implies 
that the following estimates hold for the its inverse g^ v : 



(7.0.4a) 
(7.0.4b) 



(7.0.4c) 
(7.0.4d) 



\9 00 + 1| < 4r?, 



\ a=l 



2 < gc^e'™ 



k J2\90a\ 
\ a=l 



\g 0a g 0a \ <2c ie - 2n ^2\go a \ 



a=l 



^S ab X a X b < e 2n g ab X a X b < ^-5 ab X a X b , 



2 

3^r j ----- 2 

Proof. Lemma 17.0.21 is proved as Lemma 7 in [28 



y(x u x 2 ,x 3 ) g 



□ 
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8. The Bootstrap Assumption for Sn and the Definition of N,r] min and q 

In this short section, we define the quantities N,rj m i n , and q. We then introduce some assumptions that will 
be used in our derivation of the estimates of Sections [9] and [10] 

8.1. The definition of N and the assumption Sjy < e. For the remainder of the article, we will assume 
that N is an integer subject to one of the following requirements: 

(8.1.1) N > 3, (this is large enough for all of our results except some of the conclusions of Theorem 112. ip . 

(8.1.2) N > 5, (for the full results of Theorem QUE] to be valid). 

We require N to be of this size to ensure that various Sobolev embedding results are valid; see also Remark 
112.0.21 

In our global existence argument, we will make the following bootstrap assumption: 

(8.1.3) S N [g,d$]<e, 

where Sat is defined in (]6.1.2f|) . and e is a sufficiently small positive number. Above we have written Sjv[<7, 93>] 
to emphasize the dependence of Sat on (g,d&). Observe that Syv[g, d$] = 0, where g,d& is the background 
solution from Section HJ i.e., Sjv[<7, d$] measures how much (g, «3>) differs from the background solution. 

8.2. The definitions of r\ m i n and q. 

Definition 8.2.1. Let us apply Lemma [6.2.1l to each of the equations ([5.4. 2aft - (|5.4.2cp . denoting the constant 
rj produced by the lemma in each case by 7/00)%*! an d respectively. Furthermore, let rjo be the constant 
from Lemma [7.0.21 We now define the positive quantities (recalling that <^ x <C 1 when <C c s <C \J 1 /3) Tj m { n 
and q by 

dcf 1 

(8.2.1) rj min = -minjl, rjo, Voo, %*,??**}, 

o 

(8.2.2) q = f -min{j7 min , x, 1 - x}. 

We remark that r/ m j n and q have been chosen to be small enough so that the bootstrap argument for global 
existence given in Section [11.21 will close. In particular, inequality (|7.0.4ap . with r\ < r] m i n , guarantees that the 
energies £^^[-,d(-)] for solutions to (|5.4.2ap - (|5.4.2c|) have the coercive property (|6.2.4p . 

Remark 8.2.1. By Sobolev embedding and Lemma |7.0.2[ if e is small enough, then the assumption Sat < e 
implies that there exists a constant C > such that 



(8.2.3a) \g 00 + 1| < Ce, 

3 

(8.2.3b) £> ;| < Cee^-W. 

i=i 

Therefore, if e is sufficiently small, the inequalities ( 17.0. lap and (17.0. lcp are an automatic consequence of 
(|7.0.1bp and the assumption Sn < e. 

9. Sobolev Estimates 

In this section, we use the bootstrap assumptions of Sections [7] and [8] to derive estimates of all of the terms 
appearing in the modified equations (|5.4.2ap - ([5.4. 3p in terms of the norms defined in Section 16.11 The main 
goal is to show that the error terms are small compared to the principal terms, which is the main step in closing 
the bootstrap argument in the proof of Theorem 111.11 More specifically, in Section lll.l[ the estimates of this 
section will be coupled with the energy inequalities of Corollary 16.2.31 and Corollary 16.3.21 in order to derive a 
system of integral inequalities for the solution. We divide the analysis into two propositions: Proposition 19.1.11 
provides estimates for g and c\<&, while Proposition 19.2. II provides estimates for the error terms, the reciprocal 
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acoustical metric, and for d^Q, (fi = 0,1,2,3). In particular, the latter proposition provides estimates for the 

ratio Zj == e g^— , (j = 1, 2, 3), which are crucial for closing the bootstrap argument of Theorem lll.il The main 
tools for proving the propositions are standard Sobolev-Moser product-type estimates, which we have collected 
together in the Appendix for convenience. 

. In this section, we state and 



9.1. Estimates of g^,g^ u , Qi,u = 0,1,2,3), ||e , « n ^* - tf||z,°°, and ^ 
prove the first proposition that will be used to deduce the energy inequalities of Section 111.11 

Proposition 9.1.1. Let N > 3 be an integer, and assume that the bootstrap assumptions H7.0. laH - (|7.0.1cp 
hold on the spacetime slab [0,T) x T 3 for some constant c\ > 1 and for r] = r] m in- Then there exists a constant 
e' > and a constant C > 0, where C depends on N,c\, and r] m i n , such that if and Sj\r(i) < e' on [0,T), then 
the following estimates also hold on [0, T), where hj k = e~ 2Q gj k '■ 

9.1.1a) HffooIlL^ < 2, 

9.1.1b) \\g oj \\ Loa <Ce^ n S g , N , 

9.1.1c) HSjfclU- < Ce 2n , 



9.1.2a) \\dg°°\\HM-i < Ce~ qQ S g . )N , 

9.1.2b) < 5, 

9.1.2c) \\g 00 + 1\\ HN < Ce~ qn S g]N , 

9.1. 2d) \\g 00 + l\\ Loa < Ce- qn S g]N , 

9.1.2e) \W k \\L°° < Ce~ 2Q , 

9.1.2f) \W k \\ H N-! < Ce- 2n S g . N , 

9.1.2g) < Ce~ 2n Sg., N , 

9.1.2h) \\g 0j \\ H N <Ce^ 1+ ^ n S g -, N , 

9.1.2i) \\g 0j \\cl<Ce^ 1+ ^ n S g , N , 

9.1.3a) \\d t g jk - 2cjg jk \\ H N < Ce^ n S h „. N , 

9.1.3b) \\d t g jk - 2ojg jk \\ c i < Ce^ n S htt]N , 

9.1.3c) \\dtg jk \\cl ~ Ce ™' 

9.1.4a) \\g aj d t g ak - 2u5{\\ h n < Ce~" n S g , N , 

9.1.4b) \\g aj d t9ak - 2^|| L oo < Ce-i n Sg. N , 

9.1.5a) \\d t g jk + 2ujg jk \\ H N < Ce-^ 2+q)n S g , N , 

9.1.5b) \\d t g jk + 2^' fc || L ^ < C e -( 2 +^ S g . N , 

9.1.5c) ||^ 00 || L oo < Ce~ qQ S g , N , 

9.1.5d) \\dtg 0j \\L°° <Ce-( 1+ ^Sg; N , 

9.1.5e) \\d t g jk \\L°° < Ce~ 2n , 
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(9.1.6a) \\g ab T ajb \\ HN < Ce^ n S h „, N , 

(9.1.6b) \\9 ab T a] b\\H^ < CS h „ ;N . 

Furthermore, the following estimates for dt<& hold: 

(9.1.7a) We^dtQ - f < CSg^N, 



1 



< (^ + CS d *. >N )e* Q . 

L°° \W / 



< 917b > 8 «* 

In the above estimates, the norms Sh**-,N, Sg-,N, an d Sq^-n are defined in Definition 1 6. 1 . 1\ 

Proof. Most of these estimates can be found in the statements and proofs of Lemmas 9, 11, 18, and 20 of |28| . 
The exceptions are (I9.1.5dj) . (|9.1.5ej> . (|9.1.7aj) . and (I9.1.7bj) . For brevity, we do not repeat all of the details of 
the estimates that are proved in [28J. 

Remark 9.1.1. Throughout all of the remaining proofs in this article, we freely use the results of Lemma [4. 2. 11 
the definitions of the norms from Section [6.11 the definitions (|8.2.ip . (|8.2.2|) of r] m in an d q, and the Sobolev 
embedding result H M+2 (T 3 ) C^ // (T 3 ), (M > 0). We also freely use the assumption that Sjy, which is 
defined in (|6.1.2f| « is sufficiently small without explicitly mentioning it every time. Furthermore, 
the smallness is adjusted as necessary at each step in the proof. To avoid overburdening the paper 
with details, we don't give explicit estimates for how small S^r must be. We also remark that as discussed in 
Section 12. 4[ the constants c, C, C* that appear throughout the article can be chosen uniformly (however, they 
may depend on N) as long as Sat is sufficiently small. Finally, we prove statements in logical order, rather than 
the order in which they are stated in the proposition. 

Before beginning the proof, we observe the following Counting Principle, which provides a very useful heuristic 
guideline for the estimates concerning the spacetime metric, its inverse, and its derivatives. This tool is only 
intended to help guide the reader through the estimates; we will provide full proofs of many of the estimates. 

Counting Principle 

Consider a product which contains as factors metric components g^ u , the inverse metric components g^ v , and 
first derivatives of these quantities. If u denotes the total number of upstairs spatial metric among these 
factors, and d denotes the total number of downstairs spatial metric indices, then the expected contribution 
to the rate of growth/ decay of the H N norm of the product coming from these terms is no larger than e ^( d - u ) , For 
purposes of counting, a spatial derivative dj is considered to be a downstairs spatial index, while time derivatives 
are neutral. We remark that by these criteria, hj^, (j,k = 1,2,3), is an order 1 term that doesn't contribute to 
the decay rate. Furthermore, each factor in a product, excluding hj^ but including d a hjk, that is equal to one of 
quantities under an H N norm in definitions (|6.1.1a|) - ()6. 1 . lc[) contributes an additional decay factor of e~ qn . 

As an example, we remark that ||5 a6 r a j;,||^]v grows at a rate no larger than eS l ~ q ^ (compare with inequality 
(|9.1.6a|) ). As a second example, we note the identity g 00 + 1 = ^[(<?00 + 1) — 9° a 90a], the right-hand side of 
which allows us to conclude that \\g 00 + 1\\hn grows at a rate no larger than e~ qQ (because goo + 1 is one of the 
factors under an H norm on the right-hand side of (I6.1.1ap ). The net effect is that with only a single exception, 
the products appearing on the right-hand sides of the energy inequalities of Corollary 16.2.31 and Corollary 16.3.21 
can be bounded by e~ qQ times one or more of the norms (which, according to our assumption Sjy < e 7 , should 
be thought of as size < e'). The one exception is a term in inequality (I6.2.12bj ). namely the dangerous term 

Consequently, we have to pay special attention to it during the global existence 
argument of Section [Til see Remark 111. 1.11 

Remark 9.1.2. The Counting Principle applies only to the H N norms of products of the metric components, 
the inverse metric components, and first derivatives of these quantities. For example, one should not estimate 
the growth rate of \\d a dbgoo\\ H N-i by simply observing that there are two spatial indices downstairs. The correct 
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rate in this case, which can be discerned from definition (|6.1,lap . is e^~ q ' n . Furthermore, the Counting Principle 
is not intended to be used on the fluid quantities. 

Proof of (|9.1.1a|) - (|9.1.1cl) : The estimates (|9.1.1al) and (|9.1.1bl) follow from the definition (|6. 1.211) of S N and 
Sobolev embedding. The estimate (19.1.1 eft follows from the bootstrap assumption (17.0. f bj) . 

Proof of (|9.1.2bp . (|9.1.2ep . and the preliminary estimate \\g 0j \\ L °° < Ce~^ 1+q ^ n S g . N : The L°° estimates ()9.1.2bj) . 
(I9.1.2ep . as well as the estimate || Sf - 7 < Ce~^ 1+q ^S g -N, which we will need shortly, follow from the definition 
(|6.1.2dp of S g -]\r, Lemma I7.U.11 and Lemma I7.U.2I 

Proof of (|9.1.2ah . (I9.1.2fft . ( |9.1.2gP , (I9.1.2hft . and (I9.f.2ip : The proofs of (j9.1.2aj) . (19.1.211) . and (I9.1.2hp begin 
with the fact that d$g^ u is a linear combination of terms of the form 



(9.1.8) ^V 3 " 1 . . . g^-ig™"(d Sl g XlKl )- ■ ■ (d dn g XnKn ), 

where a.\ + ■ ■ ■ a n = a, where each |d<j| > 0. We remark that (|9.1.8p can be shown inductively, using the identity 
dxlf v = -g ^g^d X g a(3 . 

To prove (j9.1.2a|) . we must estimate (|9X8|) in I? for 1 < \a\ < N, with fx,v in (f9~l~8|) equal to 0,0. To 
this end, we first bound the terms g^ l0 g X2Kl ■ ■ ■ g*n K n-igOK n m j^oo^ anc j estimate the remaining product 
(<9qi5Aiki) • • • (da n g\ n n n ) in L 2 using Proposition IA-21 The L°° terms are bounded using (|9.1.2bp . (|9.1.2ep . and 
the estimate ||c/ 0j ||/,oo < Ce~^ 1+q ^S g -N shown above. The L 2 terms are controlled by S g -N. The only difficulty 
is keeping track of the powers of e n , which Ringstrom accomplishes inductively through a counting argument 
that is analogous to the Counting Principle stated above; the details can be found in the proof of Lemma 9 of 
[28] . The proofs of (I9.1.2fp and ()9.1.2hl) are similar. Inequalities ( |9.1.2g[ ) and (I9.1.2il) now follow from (I9.1.2fl) . 
(|9.1.2hp . and Sobolev embedding. 

Proof of (|9.1.2cl) - ()9.1.2dl) : The estimate (|9.1.2cl) follows from the identity g 00 + 1 = ^[Ooo + 1) - g° a goa], 
Corollary IA-31 with v = goo and F(goo) = ^ in the corollary, Proposition \A-5\ the definition (16. 1 .2d|) of S 9 -n, 
()9.1.2hp . and ()9.1.2ip . Inequality (j9.1.2d[) then follows from (|9.1.2cp and Sobolev embedding. 

Proof of (|9.1.3ap - (j9. 1 .3c[) : The estimate (|9.1.3a|) follows directly from the definition ()6.1.2c|) of Sh tt -N and 
the observation that dthjk = e~ 2n (dtgjk — ^gjk)- Inequality (|9.1.3bp then follows from (|9.1.3ap and Sobolev 
embedding. Inequality (|9.1.3cp now follows from (j9.1.1c|> . ()9.1.3bp . and Sobolev embedding. 

Proof of (|9.1.4ap - (|9.1.4b|) : Recall that 5{ = g 0j gok + g aj g a k- Using this fact, the estimate (|9.1.4ap follows from 
Proposition [A31 the definition (j6.1.2d|) of S g - N , and (|9.1.2ep - (|9.1.3bp . Inequality (j9.1.4b|) then follows from 
(|9.1.4ap and Sobolev embedding. 

Proof of (|9.1,5ap - ()9.1.5ep : To prove (I9.1.5ap . first note the identity dtg^ k = —g^g^dtgafi, which was mentioned 
below the statement of (|9.1.8p . We then use Proposition IA-5|. the definition (|6. 1 . 2d[) of S g -N, Sobolev embedding, 
(|9.1.2ap . (I9.1.2bh . (|9.1.2ep . (j9.1.2f|l . (I9.1.2hp . and (I9.1.4ap to conclude that 

(9.1.9) 

||<9^' fc + 2u:g^\\ HN < \\g bk {g a ^ d t9ab - 2^)||^ + \\g 0j g bk d t go b \\ H N + \\g aj g 0k d t g 0a \\ H N + \\g oj g 0k d t g 00 \\ H » 

which gives (I9.1.5ap . Inequality (19.1.5bft now follows from (|9.1.5a|) and Sobolev embedding. Inequality (j9. 1 .5e|) 
then follows from (|9. 1.2e|) and (|9.1.5bp . The proofs of (|9.1.5cp and ()9.1.5dp are similar, and we omit the details. 



Proof of (|9.1.6ap : To prove (|9.1.6ap . we first use Proposition IA-51 to conclude that 
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(9-1-10) \\9 ab T a jb\\H N ^ C{ Il5 ab ||x,°° llr'ajffcH/f at + ||(?5 ab ||^iV-i ||r aj -;,||ioo ]. 

Recalling that F a jt, = ^[d a gbj + %9aj — djg a b) and that gj k = e 2Q hj k , and using (j9. 1 .2e|) . (|9. 1 .2f|) . the def- 
initions (|6.1.2c|) and (|6.1.2f'P of Sh„- t N an d Sat, and Sobolev embedding, it follows that the right-hand side of 
(19.1.101) is bounded from above by Ce^- q "> n S h „- N . This proves (|9.1.6ap . The proof of (|9,1.6bp follows similarly. 

Proof of (|9.1.7ap - (|9.1.7bp : Inequality (|9.1.7ap follows immediately from the definition t|6.1.2e[) of Sqq-^n and 
Sobolev embedding. Assuming that Sg$-N is sufficiently small, inequality (|9. 1.7bj) follows trivially from (|9.1.7a|) . 
This concludes the proof of Proposition 19.1.11 

□ 

9.2. Estimates of the error terms, m tJ,l/ , and d^Q, (/i, v = 0, 1, 2, 3). In this section, we state and prove the 
second proposition that will be used to deduce the energy inequalities of Section 111.11 

Proposition 9.2.1. Let N > 3 be an integer, and assume that the bootstrap assumptions (|7.0.1ap - (|7.0.1cp 
hold on the spacetime slab [0, T) x T 3 for some constant c± > 1 and for r] = ?] m in- Then there exists a constant 
e" > and a constant C > 0, where C depends on N,c\, and r) m i n , such that if and S^(t) < e" on [0, T), 
then the following estimates also hold on [0, T) for the quantities Aa aui AcoO) an d ^Cfijy defined in (I5.4.10ap 
- (|5.4.1Ucp and (|5.4.12a|) - (I5.4.12bp .- 



(9.2.1a) 


|Aa,00 


\h n 


< Ce- 2 « Q S 2 g . N , 


(9.2.1b) 


W^Afijl 


\h n 


< Ce^ 2 ^ n S 2 . N , 


(9.2.1c) 


\\^A,jk\ 


\h n 


< Ce( 2 ~ 2 ^ n S 2 . N , 


(9.2.1d) 


\\^c,oo\ 


\h n 


< Ce~ 2 ^S 2 . N , 


(9.2.1e) 


ll A C,0j| 


\h n 


< Ce^ n S 2 . N . 



Additionally, for the quantities A^ u defined in (|5.4.4ip - (|5.4.4kp . we have the following estimates on [0,T) : 

(9.2.2a) IIAooIIh" < Ce~ 2qn S N , 

(9.2.2b) IIAojUz/iv < Ce^-^^Sjy, 

(9.2.2c) \\A jk \ H N < Ce~ 2qn S N . 

For the commutator terms from Corollary \6. 2.,% we have the following estimates on [0,T) : 

(\3\<N), 

(H < N )- 

For the terms from Corollarv \6. 2., SI where Ag.^) [v, dv] is defined in (|6.2.6D . we have the following estimates 
on [0, T) : 

(9.2.4a) e 2 ^||A £ . (700)5oo) [%( 5o o + l),d(d 3 g 00 )]\\ Ll < Ce- qU S goo+1 . N S N , (\a\ < N), 

(9.2.4b) e 2(, - 1)n ||A £;( T t) ., tfo . ) [a S flb J -,S(aafldi)]||Li < Ce~ qn S g0t]N S N , (|a| < N), 

(9.2.4c) e 2qQ \\A £ . im [0,d(d s h jk )]\\ Ll < Ce^ n S htt . N S N , (|a| < N). 



(9.2.3a) \\[6 g ,d s ](g 00 + 1)\\ L 2 < Ce~^ l S N , 

(9.2.3b) \\[n g ,d s ]g 0j \\ L 2 < Ce^ 2 ^ n S N , 

(9.2.3c) \\[n g ,d s ]h jk \\ L 2 < Ce~ 2qn S N , 
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For the error terms corresponding to the fully raised Christoffel symbols of Lemma \ 5. 4 ■ 5\ we have the following 
estimates on [0, T) : 



(9.2.5a) 
(9.2.5b) 
(9.2.5c) 
(9.2.5d) 
(9.2.5e) 
(9.2.5f) 



a 000 1 1 ^ r< -i n Q 

^(F) \\H N S L-e Jg;N, 

*$\\h» < Ce~^ n S, 



'g;N, 

< Ce-^ n S g , N , 

< Ce~^ n S g . N , 
g;N- 



\A°$\\ HN < Ce-^ n S, 

I A Oj/C II 

l A (r)llff* 
^(r) Wh n 



For the reciprocal acoustical metric error terms A( m ), 
respectively, we have the following estimates on [0, T) : 



an 



d A a $ defined in (pDUf)) . (|5.4.4gD, (|5.4.4hp 



(9.2.6a) 
(9.2.6b) 
(9.2.6c) 
(9.2.6d) 
(9.2.6c) 



\A( m )\\H N < Ce qn SN, 
|A( m )||L°° < Ce~ qQ SN, 

\A{l } \\ L oo<Ce-^ n S N , 



e KQ \\A d <s,\\ H N < Ce~ qa S 



N ■ 



For the reciprocal acoustical metric components m? k and m ^ defined in ()5.4.4cp and (|5.4.4dp . we have the 
following estimates on [0, T) : 



(9.2.7a) 
(9.2.7b) 

(9.2.7c) 

(9.2.7d) 
(9.2.7e) 
(9.2.7f) 
(9.2.7g) 



\m 0j \\ H N < Ce~ n S N , 



\m 



i < Ce~ n S N , 



m 



jk 



1 



2s + 1 



jk 



m jk \\L°° < Ce 



-2Q 



C-^XaX,, < e 2n m ab X a X b < C5 ab X a X b , 
\\dm jk \\ H N-i < Ce~ 2n S 



V(x 1 ,x 2 ,x 3 ) e 



AT, 



\\dm Jk \\ L oo < Ce~ Zil S N . 
For the commutator terms from Corollary \6. 3.21 we have the following estimates on [0, T) : 



(9.2.8) m m ,d s ]n L * < Ce~^ n S 2 N , (\a\ < N). 

For the Square of the enthalpy per particle a = —g a P(d a $>)(dp<&), we have the following estimates on [0,T) 



(9.2.9a) 
(9.2.9b) 

(9.2.9c) 



cr > 0, 

I ~-ll s~i „ — 2xVL 

1 



a 



< Ce 



2xQ. 
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For 3& and the ratio Zj = e 37^— , (j = 1, 2, 3), defined in (|5.4.4e|) . we have the following estimates on [0, T) : 

(9.2.10a) < CS d ^. N , (1 < |a| < JV - 2), 

(9.2.10b) 11%^-llioc, < Ce {K ~ 1)n S d ^ iN , (\a\ < N - 3), 

(9.2.10c) Ikjllffw < CSq^-n. 

For the time derivatives of the fluid-related quantities, we have the following estimates on [0, T) : 



(9.2.11a) 


\\d?* 


|l°° 


< Ce - *", 


(9.2.11b) 


WdtZj 






(9.2.11c) 


\\9t^(m) 




< Ce~ qn S N , 


(9. 2. lid) 


11 1 (m) 




< Ce~^ n S N , 


(9. 2. lie) 






< Ce-^ n S N , 


(9.2.11f) 


+ 2wm J ' fc 


L°° 


< Ce~^ n S N , 


(9.2.11g) 


\\d t m jk 


L°° 


< Ce~ 2n S N . 



In the above estimates, the norms S gm+1 . N , S gQji . N>S.h**;Ni S 900+ i-n, S go *;N, Sh**;N, Sg;N, Sqq.n, and Sn are 
defined in Definition \6.1.1\ 

The proof of Proposition 19.2.11 is located in Section 19.41 Before we prove the proposition, we first prove 
some Lemmas that will be useful for estimating quantities involving d&. We also prove the estimates (I9.2.9ap - 
(|9.2,9cp . since we need them for the proofs of some of the Lemmas. See Remark 19. 1.1 1 for some conventions that 
we use throughout our proofs. 

9.3. Some preliminary lemmas for estimating d^, (/U = 0, 1,2,3), and proofs of (|9,2.9ap - (j9.2.9c|) . 

Lemma 9.3.1. Under the assumptions of Proposition [9. 2. 1\ and Remark \9.1.1[ there exists a constant C > 
such that 

(9.3.1) \\a-a\\ H N < Ce- 2 * n Sjv. 

Proof. Recall that a = -g a P(d a <&){dp$>) and a = -g a ^ i (5 a $)(9^$) = (^e~ xQ ) 2 . We therefore can decompose 
a — a as follows: 

(9.3.2) a - a = (d t $ - * e -* n )(d t $ + te - *") - (g 00 + l)(d t $) 2 - 2g 0a (d t $){d a <S>) - g ab (d a <S>){d b $). 

Inequality A9.3. lf> now follows from Proposition IA-5[ the definition (I6.1.2fj) of Sat, Sobolev embedding, (19. 1 .2cH . 
()9.L2ej) . (l9X2fl) . (|9.1 2hjl . and (19.1. 7a|) . 

□ 

Lemma 9.3.2. Under the assumptions of Proposition [9. 2. II and Remark \9.1.1[ there exists a constant C > 
such that 



(9.3.3a) 
(9.3.3b) 

(9.3.3c) 
(9.3.3d) 



a > 0, 

< Ce~ 2><n , 

< Ce 2 ^, 



a 



L°° 



[a )\\ h n- 



.1 < Ce~ 2Ksn S N . 
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Proof. Using the fact that a = (#e~ ) 2 , (|9.3.ip . and Sobolev embedding, it follows that if Sn is sufficiently 
small, then 

(9.3.4) a > 0, 

(9.3.5) |H| L oo < Ce' 2 ^, 

(9.3.6) 



1 



< Ce 2 * n , 



a 

which implies (|9.3.3ap - (|9.3.3cp (and which are the same as (|9.2.9ap - (|9.2.9cp ). 

To prove (19.3.3dp . we first note that because a is a function of t alone, it follows from (19.3. ip that 

(9.3.7) \\da\\ HN -i = \\B(a - o)\\ h n-i < \\a - a\\ H N < Ce~ 2>(n S N . 
By Corollary IA-31 with v = a and F(a) = a s in the corollary, we have that 

N 

(9.3.8) \\d{o- s )\\ H N-i <C\\da\\ H N-iJ2h S ~ l h™M l L~- 

i=i 

Inequality (I9.3.3dp now follows from ()9.2.9bp . (|9.2.9cD . (19X71) . and (HEMP - 

□ 

Corollary 9.3.3. Under the assumptions of Proposition [9.2. 1\ and Remark \9.1.1[ there exists a constant C > 
such that 

(9.3.9) \\a s - a s \\ H N < Ce~ 2xsn S N . 

Proof. By Taylor's theorem, we have that a s — a s = s(a*) s ~ 1 (a — a) for some a* lying between a and a. By 
LemmaEXSl it follows that (a*) 8 " 1 < Ce" 2 *^ 1 ^. By LemmaEXH it therefore follows that 

(9.3.10) \\a s - a s \\ L2 < Ce" 2 ^- 1 ^^ - a\\ L 2 < Ce~ 2>esa S N . 

To bound the spatial derivatives of o~ s — a s in L 2 , we note that a is a function of t alone. Therefore, by 
(|9.3.3dp . it follows that 

(9.3.11) \\B(a s - a^WvN-! = \\d(a s )\\ H N-x < Ce~ 2 * sn S N . 
Combining (I9.3.10p and (19. 3. lip , we deduce (l9~3T9jh 

□ 

Corollary 9.3.4. Let f(d<$>) be defined as in ()5.4.41[) : 

(9.3.12) f(d$) = f 2a s (d t <S>) 2 —a s+1 . 

s + l 

Then under the assumptions of Proposition \9.2.1\ and Remark \9.1.1\ there exists a constant C > such that 

(9.3.13) - f(d$)\\ HN < Ce- 2 ^ S+1 ^S N . 

Proof. Using Proposition EH the definition (16.1. 2fp of Sat, Sobolev embedding, (I9.1.7aj) . (I9.2.9b|) . (I9.3.3dp . and 
(|9.3.9p . we estimate 



(9.3.14) ||/(d$) - f(d$)\\ HN < C{\\(d t $ - *e-^)(0 t $ + Ve- M )a s \\ HN 

+ (^e^) 2 ))^ - a s \\ H N + \\a s+1 - a s+1 \\ HN } 
< Ce- 2 ^ S+1 ^S N . 
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Proof. We now provide a proof of Proposition 19.2.11 We prove statements in logical order, rather than in the 
order in which they are listed in the conclusions of the proposition. See Remark 19. 1.11 for some conventions that 
we use throughout the proof. 

Proofs of (19.2. lOaj) - (j9.2. lOc^ : By Sobolev embedding and the definition (I6,l.ldp of Sg$.\a\+2i the following 
inequality holds for 1 < \a\ < N — 2 : 

(9.4.1) \\dtd 3 $\\ L ao < Ce-" n S dHd \ +2 . 

Integrating dtd^Q in time, using (I9.4.ip . and using that Sg^.^^it) 1S increasing, we conclude that 



1+2 



(*), 



(9.4.2) \\d S mh~ < c(^.| S | +2 (0) + f\^ n ^S dnd \ +2 {r)dT) < CS dHS 

which implies inequality (j9.2.10aj) . 

To prove (|9.2.10b[ ) and (j9.2.10c|) . we use the definition (16.1.2ep of Sq§-n, Corollary IA-31 with v = dt$ and 



~ 37¥ * n coroi l ar y) together with (|9,l,7ap and (|9.1.7bp to conclude that 



N 



(9.4.3) 



l) 



HN . X ^ Ce-^ n S d ^ N J2e K{l+1)Q e-^ 1 -^ = Ce" n S d ^ N . 



i=i 



By (|9.1.7bp . Sobolev embedding, and (l9X3|) . it follows that 

1 



(9.4.4) 



do 



< Ce 



(\a\ < N- 2). 



Now by the definition Zj = 6 , the product rule, Proposition [A^l (|9.1.7bp . (19.2. lUal) . (19X31) . and (f9X4"D . 
we have that 



(9.4.5) WdsZjUo* <Ce^ n S d ^ N , 

1 



(9.4.6) 



\Zj\\ H N 



< Ce 



\d$\\ H N + ||5$||lc 







(\a\<N-3), 
f-1 \ < CS d * ;N , (j = 1,2,3), 



which proves (19.2.10bft and (I9.2.10cl) . 

Proofs of (|9.2.1ap - (|9.2.1ep : To prove (|9.2.1ap . we first recall equation (|5.4.10ap : 

(9.4.7) A Ai00 = (g 00 ) 2 {(dtg 00 ) 2 - (r 00 o) 2 } + g 00 g 0a {2(d t goo)(d t g 0a + d a g 00 ) - 4r 00 or oa} 
g 00 g ab {(dtgoa)(d t g ob ) + (d a g o)(d b g o) - 2r 00a r 00b } 
g 0a g ob {2(d t g o)(d a g ob ) + 2{d t g 0b ){d a g m ) - 2T ma T a()b - 2r ofer oa} 
g ab g ol {2(d t g 0a )(d igob ) + 2(d b g 00 )(d a g 0l ) - 4r 00 ar /0b } 



+ . 



+ g ab g lm (d a g i)(d b g 0m ) + -g lrn ( g ab d t g al - 2 W tf )(d b90m + d m9ob ) 

2 - „ ' 



e 2n g ab d t h al -2ujg° b g i 
1 

4 ( 



7 g ab g lm (dagoi+digoa)(d b go m + d m g ob )-^( g ab d t9al - 2u5\ )( g lm d t g bm - 2u5 l b 



e ™ g ab dthal _2ug0bg 0l e^g^d t h brn ~2Lug0'g 0b 
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We now use Proposition the definition (|6.1.2f|) of Sat, Sobolev embedding, (j9.1.2c|) . (j9.1.2d|) . (I9.1.2ep . 
(HXHD, (|9.1.2hl) . dEH2i]), (|9.1.3aj> . (I9.1.3cl) . f|9.1.4aj) . (I9.1.4bp . and the relation T^ au = ±{8^ + d v g a ^-d a g,, v ) 
to conclude that 

(9.4.8) I|Aa,oo||^ < Ce^ n S 2 . N . 

This proves (19.2. la[) . The proofs of (19.2. \h\ - (19.2. left are similar, and we omit the details. We also remark that 
slight variations of these inequalities were proved in Lemma 12 of |28| . 

Proofs of (I9.2.2ap - (I9.2.2cj) : To prove (I9.2.2ap , we first use equation (15.4.4ij) and Proposition IA-51 to arrive at 
the following estimate: 

(9.4.9) 

IIAooIIh" < c{||A Ai00 |biv + ||A C)0 o||^ + ||/(9$) - f(d$)\\ HN + \f(d$)\\\g 00 + 1\\ h n 

+ 1 1 Poo + 1 II^II°' s+1 ||l o ° + 1 1 Poo + 1 \\L°°\\d(o- s+1 )\\ H N-i + \u- H\\\d t goo\\ H N + \uj + H\\uj - H\\\g 00 + 1\\ h n 

We now use (I9.1.2cl) . (|9.1.2d|) . (|9.2.1aj) . ()9.2.1dll . (I9.3.3dp . ()9.2.9bD . (19.3.131) . the definition ()6.1.2fl) of S N , the 
fact that f(d$) = (|±2)^2( s +i) e -2( s +i)^ ) and the triyial estimate e -2( s +i)xQ < e -2 g n to conclude that 

(9.4.10) IIAooIIh" < Ce- 2qn S N , 
which proves ()9.2.2ap . 

Inequalities ()9.2.2bp and (I9.2.2cj) can be proved using similar reasoning; we omit the details. 
Proofs of (|9.2,4ap - (j9.2.4cH : To begin, we first recall equation (16. 2. 6ft : 

(9.4.11) A g . iltS) [v,dv] = -~/H(d a g ab )vd b v - 2 1 H(d a g Qa )vd t v - 2 1 Hg Ga (d a v){d t v) 

- (d a g 0a )(d t v) 2 - (d a g ab )(d b v)(d t v) - l -{d t g™){d t v) 2 

+ (^d t g ab + ug ab ^ (d a v)(d b v) + (H - u;)g ab (d a v)(d b v) 

- 7 H(d t g 00 )vd t v - -fH(g 00 + l)(d t v) 2 . 

We now claim that the following inequality holds for any function v for which the right-hand side is finite: 

(9.4.12) \\A £ . {7>5) [v,dv]\\ L i < C{e-^||^||| 2 + e ~^ n \\dv\\ 2 L , + C h) e-« n \\v\\ 2 L2 }, 
where C( 7 ) is defined in (|6.2.4p To obtain (|9.4.12p . we use the Cauchy-Schwarz inequality for integrals, (|9. 1.2dp . 



(|9.1.2ep . (|9.1.2gp, (|9.1.2ip . ()9.1.5bp . and (|9.1.5cp . Inequalities (|9.2.4ap - (|9.2.4c|) now easily follow from definitions 



flgXHD - (l6~l~2fT) and (|9.4.12l) . 
Proofs of (|9.2.5ap - ()9.2.5f[) : To estimate A°£°, we first recall equation (|5.4.14ap : 



Ag? = ^ 00 ) 3 ^oo + \g m g 0a g Qh (d t g ab + 2d a g 0b ) 



+ \{g m ?g Qa {id tgQa + «9 a50 o) + lg 0a g ob g ol d a g bl . 

By Proposition lA^U the definition (I6.1.2dp of S g;A r, Sobolev embedding, (I9.1.2cp . (j9.1.2hjl . (I9.1.2ijl . and (19.1 .3d) , 
it follows that 

(9.4.14) ||AgJ|| H iv < Ce-« n 5 g;J v, 
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which proves (|9.2.5ap . The estimates ()9.2.5bp - (|9.2.5f|) can be proved similarly using Propositions 19. 1 ,T1 and |A^B"1 
Proofs of (j9.'2.6a|) - (j9.2.6djl : To prove (I9.2.fctp . we first recall equation (I5.4.4fl) : 

(9.4.15) A (m) = (1 + 2s)(g 00 + l)(g 00 - 1) + 2(1 + 2s)e n g 00 g 0a z a + e m (g 00 g ab + 2 S g 0a g 0b )z a z b . 

By Proposition EH the definition (|6.1.2fp of S N , Sobolev embedding, (I9.1.2cj) . (|9.1.2dp . (19.1. 2el) . (19.1.2flh 
(I9.1.2hp . (|9.1.2il) . (|9.2.10bp . and (19.2. lOcll . it follows that 

(9.4.16) l|A (m) || H iv <Ce-i u S N , 

which proves (|9.2.6ap . Inequality ()9.2.6bp now follows from (|9.2.6ap and Sobolev embedding. 

Inequality (|9.2.6cp follows from equation (5.4.4g) using similar reasoning. Inequality ()9.2.6dp then follows 
from ()9.2.6cp and Sobolev embedding. 



Proof of (|9.2.6ep : To prove ()9.2.6ep . we first multiply each side of equation ()5.4.4hp by e to deduce that 
(9.4.17) 

e^Aa* = we^{3F(A M )-x}d t $-e^F(A (m) )^ 

+ 2 s (a°°J$$ + [Ag» + A Oao + A«*]0 a $ + e 2 ^[Ag + A$ 6 + A^jz^S + e 2n A^z a z fe d c <&) j, 

where 



(9.4.18) F(A {m) )^'[(l + 2 S ) + A {m) ] \ 

and A( m ) is defined in (I5.4.4fp . 

We begin by estimating 3-F(A( m )) — x where x = 1 _^ 2g = 3-F(0) (as in (I5.4.4ap ). Using Corollary IA-4[ with 
v = A( m ) in the corollary, and (|9.2.6aj) . it follows that 



(9.4.19) 



3F(A 



(m) j 



3F(0) 



< C\\A 



(m) II ff^ 



< Ce _ « n Sv. 



As a simple consequence of (|9.4.19p . we note that 



(9.4.20) 
(9.4.21) 



F(0) 
' 1 ' 

1 + 2s 



< Ce- qn S N , 



||a(F(A M ))||^-i <Ce~^S N . 
We now estimate the terms on the right-hand side of (19.4.17p . The bound 



3F(A 



(m), 



3F(0)} 



< Ce~ qQ S N follows from Proposition EIJ the definition (I6.1.2fl) of S 



Sobolev embedding, (19.1.7ap . and (19.4.19p . 
The remaining term to be estimated is e 



f(a m ){g} 



N, 



where G is an abbreviation for the terms in 



the large braces on the right-hand side of (|9.4.17p . By Proposition IA-51 the definition ()6. 1 .2f|) of Sat, Sobolev 
embedding, (I9.1.2cl) . (I9.1.2el) . (l9~l~2fl) . (I9.1.2hl) . (I9.1.7al) . (|9.2.5ap - (l9~2~5f|) . (19.2. lOap . (I9.2.10bl) . and (19.2. lUcI) . we 
have that e^HGH^jv < Ce~ qU SN- We may therefore again apply Proposition IA-51 using (|9.4.20p and (|9.4.2ip 

to deduce that F(A (m) )|G| ^ < Ce~ qQ S N . This concludes our proof of (I9.2.6el) . 
Proofs of (|9.2.7ap - (|9.2.7g[): To prove (|9.2.7c|) - (|9.2.7gP, we use (|5.4.4cp to decompose 
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(9.4.22) mP k = F(A {m) ) 
where .F(A( m )) is defined in (|9.4.18p . and 

(9.4.23) M^ k = g jk -A^ k 

By Proposition EU (|9.1.2ej) . (|9.1.2f) . (|9.2.6ap . and (|9.2.6b|) . we have that 

(9.4.24) \\M jk \\ LO o < Ce~ 2n , 

(9.4.25) \\dM jk \\ H N-i < Ce~ 2n S N . 

It therefore follows from Proposition EH (|9.1.2ep . and (|9.4.20p - (|9.4.25p that 



(9.4.26) 



m jk —g jk 

2 S + r 



< Ce-^ n S N , 



(9.4.27) \\m jk \\L°° < Ce~ 2Q , 

(9.4.28) \\dm jk \\ H N-i < Ce~ 2n S N , 



which proves (|9.2.7cp . ()9.2.7dp . and ()9.2.7f[) . Inequality ( |9.2.7g ) now follows from Sobolev embedding and 



(|9.2.7fp . Inequality (|9.2.7el) follows from (17.0. 4dp and (|9.2.7cp . 

Inequalities (|9.2.7ap and ()9.2.7bp can be proved using the ideas similar to the ones we used to prove (|9.2.7c|) 
- ( |9.2.7gP ; we omit the details. 

Proofs of (|9.2.3ap - (|9.2.3cl) : To prove (|9.2.3ap - (I9.2.3cp . we first use equations (I5.4.2al) - (|5.4.2cl) to conclude 
that 

(9.4.29a) d 2 g 00 = (g 00 )" 1 ( - g ab d a d b g 00 - 2g 0a d a d t g o + 5Hd t g 00 + e^ 2 ^ + 1) + A 00 



(9.4.29b) d 2 g 0] = (5 00 )" 1 ( - 9 ab d a d b g 0j - 2g 0a d a d t g 0j + 3Hd t g 0j + 2H 2 g 0j - 2Hg a0 T a]b + A 0j ) , 

(9.4.29c) d 2 t h jk = (g 00 )- 1 ( - g ab d a d b h jk - 2g 0a d a d t h jk + 3H8 t h jk + A jk 



We will now estimate \\d 2 goo\\ H N-i,\\d 2 goj\\ H N~i, and \\d 2 hj k \\ H N-i , (j,k = 1,2,3). Using (I9.4.29ap . Corollary 
IA-41 with v = g 00 + 1 and F(v) = (v - l)" 1 = {g 00 )~ l in the corollary, Proposition |AJ5J the definition f)6.1.2f[) 
of Sat, Sobolev embedding, (I9.1.2ep . (I9.1.2fll . (I9.1.2hl) . and (I9.2.2ap . it follows that 

(9.4.30a) ||^ 2 5oo||^-i < Ce~ qQ S N . 

Using similar reasoning, we also deduce that 

(9.4.30b) \\d 2 g 0j \\ H N-x<Ce^ n S N , (j = 1,2,3), 

(9.4.30c) \\d 2 h jk \\ H N-i < Ce~ qn S N , (j, k = 1, 2, 3). 

Note that because of inequalities (j9.1.6aj) and (19.2. 2bj) . the inequality (I9.4.30bp has an additional power of e 
on its right-hand side compared to (|9.4.30ap and (|9.4.30cp . 

We now use Proposition EH the definition (I6.1.2fp of Sat, Sobolev embedding, (I9.1.2ap . (I9.1.2fp . f!9.1.2hp . 
(|9.4.30ap to obtain the following inequalities, valid for \a\ < N : 
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(9.4.31) 

|p s ,%]G?oo + 1)||^ < \\g m d a {d 2 g m ) - d d (g 00 d?g 00 )\\ L 2 + \\g ab d d (d a d b g 00 ) - ddg ab d a d b g 00 )\\ L 2 

+ 2\\g 0a d (1 (d t d a goo) - d s (g 0a d t d a g 00 )\\ L 2 

< c(ll^ 00 ||^-i|l^oo||^-i + II^H^-xll^^ooll^-! + \\dg 0a \\ H N-4dadtgoo\\ H ^ 

< Ce~ 2qn S N . 

This completes the proof of (|9.2.3ap . Inequalities ()9.2.3bp and (|9.2.3c|) can be proved similarly; we omit the 
remaining details. 

Proof of (|9.2.8p : To prove (|9.2.8[) . we use Proposition IA-61 the definition (|6. 1 .2f|) of Sat, Sobolev embedding, 
(|9.2.7fj> . and (|9.2.7ap to obtain 

(9.4.32) \\[nm,dsm L 2 < \\m ab d s (d a d b $) - d s {m ab d a d h $)\\ L * + 2\\nP a d <i {d t d a <S>) - d s (m 0a d t d a $)\\ L2 

< c(\\dm ab \\ H N-4d a d b <s>\\ H N-i + \\dm 0a \\ H N-t\\d a dt$\\ H K-i 

This demonstrates (|9.2.8p . 

Proof of (|9.2.11ap : To prove (|9.2.11ap . we first solve for <9 t 2 $ using f|5.4.3j) : 

(9.4.33) d 2 ^ = m jk djd k <S> + 2m 0j d t d j <S> - xwd t $ - A 9 $. 

By the definition ()6.1.2fl) of Stv, Sobolev embedding, (I9.1.7al) . (I9.2.6el) . (I9.2.7dl) . (I9.2.7bp . we have that 

(9.4.34) 11$ < Ce"* . 
This proves (|9.2.11aD . 

Proof of (I9.2.11bl) : By the definition (j5.4.4e|) of Zj, we have that 



e 



(9.4.35) d+Zj = —ujZj + e 

Using the definition (16.1.2fft of Sat, Sobolev embedding, (I9.1.7bp . (I9.2.10ap . (19.2.10bl) . (19.2. llah . (I9.4.35p . and 
Sobolev embedding, we have that 

(9.4.36) WdtZjU- < Ce^ n S d ^ N . 
This demonstrates ()9.2.11bp . 

Proofs of (|9.2.11cp - (I9.2.11d[) : We first recall equation (15.4.411) : 

(9.4.37) A (m) = (1 + 2s)(g 00 + l)(g 00 - 1) + 2(1 + 2s)e n g 00 g 0a z a + e 2n (g 00 g ab + 2sg 0a g 0b )z a z b . 

Differentiating (I9.4.37p with respect to t and using (|9.1.2dp . (|9.1.2ep . (I9.1.2il) . (I9.1.5cp . (I9.1.5dp . (I9.1.5el) . 
(|9.2.10bft . (I9.2.11bp . it follows that 

(9.4.38) II^A^Hi- <Ce-« n S N . 
This concludes the proof of (j9.2. llc|) . 
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The proof of (19.2.11d[ ) follows similarly, and we omit the details. 
Proof of (|9.2.11ej) : To prove (|9.2.11cp , we first recall equation (|5,4.4dD : 

(9.4.39) m 



g 00 g 0j (l + 2s) + 2e n [(a + l)g 0j g 0a + sg 00 g a i] z a + e 2n (g°ig ab + 2sg a i g ob )z a z b 



(l + 2a) + A (m) 

Differentiating (|9.4.39l) with d t and using (I9.1.2dp . (I9.1.2ep (|9.1.2iD . (I9.1.5cl) . (I9.1.5dl) . (|9.1.5ej> . (|9.2.6bl) . (|9.2.10bp . 
(I9.2.11bl) . and (|9.2.11cp . it follows that H^m ^' || L oo < Ce _ ( 1+9 ) n Sjv. This proves (19. 2. Hep . 



Proo/s of (|9.2.11fl) - (|9.2.11gP: First, using (I9.2.6bft and (19.2. llcp . it follows that 



(9.4.40) ||^(F(A (m) ))|| L oo < Ce-^8 N) 

where F(Ar m \) is defined in (|9.4.18p . Next, using the decompositions (|9.4.22[) and (|9.4.23|) . together with simple 
triangle inequality estimates, we have that 

(9.4.41) ||cW fe + 2aW fc || L oc < ||^(F(A (m) ))|| i o ||^ fc || Loo + ||a 4 (F(A (m) ))|| L c ||A^ ) || LO o 

+ ||F(A( m ))||x,<x. ||9 t A^||x<x> + \\F(A {m) )\\ L ^\\d t g jk + 2cog^ k \\ L ^ 
+ 2u\\F(A {m) )\\ L o \\A^ ) \\ L ~. 

Finally, using f|9.1.2ej) . ()9.1.5bl) . ()9.2.6dll . ()9.2.11dl) . (19.4.201) . and (I9.4.40p to bound the right-hand side of (|9.4.41l) . 
we conclude that 



(9.4.42) \\d t m jk + 2wm J ' fc || L o < Ce-^ +q)n S N . 

This proves (19.2. 1 lf[) . 

Inequality ( |9.2.11g ) follows from (|9.2.7dp and ()9.2. 1 lf[) . This concludes our proof of Proposition 19 . 2 . H 

□ 



10. The Equivalence of Sobolev and Energy Norms 

As is typical in the theory of hyperbolic PDEs, our global existence proof is based on showing that the 
energies of Section [6] remain finite (they happen to be uniformly bounded for t > in the problem studied 
here). However, the boundedness of the energies does not in itself preclude the possibility of blow-up; to show 
that the blow-up scenario (4) from the conclusions of Theorem 15.21 does not occur, we will control appropriate 
Sobolev norms of dg and d&. In this short section, we supply the bridge between the energies and the norms. 
More specifically, in the following proposition, we prove that under appropriate bootstrap assumptions, the 
Sobolev- type norms and energies defined in Section [6] are equivalent. 

Proposition 10.0.1. (Equivalence of Sobolev norms and energy norms) Let N > 3 be an integer, and 
assume that the bootstrap assumptions (|7.0.1ap - ([7.0. lc|) hold on the spacetime slab [0, T) x T 3 for some constant 
c\ > 1 and for ij = r\ m i n . Let ((5,7) be any of the pairs of constants given in Definition \6. 2. 1\ and let Cr^\ be 
the corresponding constant from Lemma \6. 2.1[ There exist constants e'" > and C > depending on N, c\ , 
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rjmin, 7, and 5, such that if Sn < e'", then the following inequalities hold on the interval [0, T) for the norms 
and energies defined in (16.1.2a!) - (16.1.2fp . fllT2~2]) . (16.2.9a!) - (I6.2.9dj) . (I6.3.1cl) . and (I6XTD : 

(10.0.1a) C- 1 !!!^!!^ +C {l) \\v\\ L 2 +e- n \\Bv\\ L 2} <£ {7>s) [v,dv] < C{\\d t v\\ L 2 + C M \\v\\ L 2 + e- n \\Bv\\ L 2}, 

(10.0.1b) C~ 1 Eg 0Q + l;N < Sg m + l-N < CEg 00 + \-N, 

(10.0. 1C) C l Eg Qt -M < Sg Qt] N < CEg Qf -N, 

(lO.O.ld) C^E^s < S h ^ ;N < CE h ^. N , 
(lO.O.le) C~ l E g , N < S g . <N < CE g , N , 

(10.0. If) C~ Eq$. n < Sg^ ]N < CE 9 ^. N , 
(lO.O.lg) C _1 Ejv < Sn < CB N . 

Analogous inequalities hold if we make the replacements (E goo+ \^N , E gQt ^N , E^^-^n , E g -N , Eq^n) 
-> (E.g 00 +l;NiE got]N ,E h ^. N ,E g]N ,E 9 $. N ) and (S goo+ i-N, Sg *;N, Sfi**;N> Sg-N, Sq$-n) 

— > (S_g 00 + 1;N > S_g 0t -N, S_ htt . N , S_ g .N, S_Q$. N ) . 

Proof. The proof of the proposition will follow quite easily from the estimates already derived. The inequalities 
in (llO.O.lap follow from the definition (^2]) of S^s) [ v , 9v] , the definition (j6.1.2f) of S N , jS23]), and ()7.0.4dl) . 
The inequalities in (jlO.O.lbp - ()10.0.1d|) then follow from definitions (|6.1.2ap - (|6.1.2cp . definitions (|6.2.9ap - 
(|6.2.9cp . and (|10.0.1ap . The inequalities in (110.0. lf'P follow from definitions (|6.1.2e|) and (|6.3. lc|) . and from 
(j7.0.4dp plus (I9.2.7cp . Finally (110.0. lei) and QlO.O.lgP follow trivially from definitions (16.1.2dj) . (I6.1.2fl) . (I6.2.9dl) . 



and (|6.4.ip . and from the previous inequalities. 

□ 



11. Global Existence 

In this section, we use the estimates derived in Sections [9] and [10] to prove two main theorems. In the first 
theorem, we show that the modified system (|5.4.2ap - ([5.4.3P has global solutions for initial data near that of the 
background solution (g, d&) on [0, oo) x T 3 , which is described in Section |4j As described in Section [5.6| if the 
Einstein constraint equations and the wave coordinate condition = 0, (fj, = 0, 1, 2, 3), are both satisfied along 
the Cauchy hypersurface S = {x € M \ t = 0}, then the solution to the modified equations is also a solution to 
the irrotational Euler-Einstein system. The main idea of the proof is to show that the energies satisfy a system 
of integral inequalities that forces them (via Gronwall-type estimates) to remain uniformly small on the time 
interval of existence. Since Proposition 110.0. II shows that the norms of the solution must also remain small, 
Theorem 15.21 can be applied to conclude that the solution exists globally in time. In the second theorem, we 
provide for convenience a proof of Propositions 3 and 4 of [28], which provide criteria for the initial data that 
are sufficient to ensure that the spacetime they launch is a future geodesically complete solution to irrotational 
Euler-Einstein system. 



11.1. Integral inequalities for the energies. In this section, we derive the system of integral inequalities 
that was mentioned in the previous paragraph. 

Proposition 11.1.1. (Integral inequalities) Let N > 3 be an integer. Assume that on the spacetime slab 
[0, T) x T 3 , (g^y^d^), {n,v = 0,1,2,3), is a classical solution to the modified system (I5.4.2ap - (|5.4.3p . and 
that the bootstrap assumptions (j7.0.1aj) - (I7.0.1cj) hold for some constant c\ ^ 1 and for rj — f]min- 

Then there 
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exist constants e"" > and C > 0, where C depends on N,c\, and rj m i n , such that j/Sjv(i) < e"" on [0, T) and 
ti £ [0, T), i/ien the following system of integral inequalities is also satisfied for t G [ii,^) : 



(11.1.1a) ^ ;JV (t) < ^ ;JV (*i) + C / e-^E^(r) dr, 
(11.1.1b) 

^oo+i;iv(*) < ^oo+i ; jv(<i) + e^TS^r) dr, 
(11.1.1c) ^L;iv(*) < ^L;iv(*0 + / -^^.^(r) + CE h „, N {r)E g ^ N {r) + Ce^E^r)^.^) dr, 
(ll.l.ld) EL ;J v(i)<^L;iv(ii) + C f e~ qHT V 2 N {T)dT. 

Proof. We apply Corollary EL2l using (I9.2.6el) . fl9.2.7gD , (19.2.7blh ([9X8]) . and (19.2.11fp to estimate the terms 
on the right-hand side of (|6.3.8p . using Proposition 110, OTl to replace the norms with corresponding energies, and 
dropping the term (x — l)u Yl\d\<N Jt 3 e2>c ^ mab {dadd^)(dbda^) d?x on the right-hand side of (|6.3.8p . which by 
(|9.2.7ep is non-positive for x < 1, thereby arriving at the following inequality: 



(11.1.2) -(^ N (t))<Ce-« m B 2 N . 

Integrating from t\ to t gives (jll.l.lap . 

To prove (jll.l.lcj) . we apply Corollary [Oj2 using tj9.1.6aj) . (|9.2.2b|> . (|9.2.3bp . and (|9.2.4bp to estimate the 
terms on the right-hand side of (|6.2.12bp . using Proposition IIO.OTTI to replace the norms with corresponding 
energies, and using definition (|8.2.2p to deduce that 2(q — 1) — 770* < — 4g, thereby arriving at the following 
inequality: 



(11.1.3) -(tf g ^ N (t)) < -AqHE*^. N {t) + CE h „, N E got;N (t) + Ce^ m E N (t)E g ^. N (t). 

Inequality flll.l.lc|> now follows by integrating from t\ to t. Inequalities (jll.l.lbj) and (jll.l.ldj) can be proved 
similarly; we omit the details. □ 

Remark 11.1.1. The term CE h ^. N E_ g0f . N in inequality (jll.l.lcp arises from the S_ 9Qt -N S?=i e ^ q ~ 1 ^\\9 ab '^ ajb\\H N 
term on the right-hand side of (|6.2.12bp . This term is dangerous in the sense that it does not contain an expo- 
nentially decaying factor, and looks like it could lead to the growth of E gQif . N . However, as we shall see in the 
proof of Theorem lll.l| there is a partial decoupling in the integral inequalities in the sense that the CE^^-n 
factor in the dangerous term can be controlled from inequality (jll.l.ldp alone. We will then insert this infor- 
mation into inequality (jll.l.lcp . and also make use of the negative term —4qHE? 0f . N to obtain a bound for 

For completeness, we state the following version of Gronwall's inequality; we omit the simple proof. We will 
use it in Section 111.21 

Lemma 11.1.2. Let b(t) > be a continuous function on the interval [ti,T], and let B(t) be an anti- derivative 
ofb(t). Suppose that A > and that y(t) > is a continuous function satisfying the inequality 

(11.1.4) y{t)<A+[ b{T)y(r)dT 

Jti 

for t G [ti, T]. Then for t G [ii, T], we have 



(11.1.5) 



y(t) < Aexp B(t) - B(t 
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In addition, in Section 111.21 we will apply the following integral inequality to (jll.l.lcj) in order to estimate 
the energy E got - N (t). 

Lemma 11.1.3. Let b(t) > be a continuous non- decreasing function on the interval [0,T], and let e > 0. 
Suppose that for each t\ G [0,T], y(t) > is a continuous function satisfying the inequality 

(11.1.6) y\t) < y 2 (h) + / -b(r)y 2 (T) + ey(r) dr 

Jr=ti 

for t G [ti, T]. T/ien /or any ti,t £ [0, T] wi/i ii < t, we have that 

(H-1-7) y(t)<y(h) + - s ^r. 

Proof. Let C be the "highest" curve in the (t, y) plane on which the integrand in . 1 .6j) vanishes; i.e. C = 
{(t,y)\y = jm}- Then by (|11.1.6p . above C (i.e. for larger y values), y(t) is strictly decreasing. Let y(t) have 
achieve its maximum at t max G [ii,T]. We separate the proof of fjl 1 . 1 . Tj) into two cases. Case i) assume that 
tmax = ti- Then y(t) < y{t max ) = y{t\) for t G [t%, T], which implies (|11.1.7|) . Case ii) assume that t max G (t\,T\. 
We claim that y(t max ) < — r. For otherwise, the point {t max ,y{t max j) lies above C. Since y(t) is then strictly 
decreasing in a neighborhood of t max , it follows that there are times t* < t max , with t* G (t\ ,T), at which 
y(i*) < yitmax)- This contradicts the definition of t max . Using also the fact that my is non-increasing, it follows 
that y{t) < y(t max ) < — r < this concludes the proof of (|11.1.7|) . □ 

11.2. The global existence theorem. In this section, we state and prove our main theorem, which provides 
global existence criteria for the modified system ()5.4.2ap - (|5.4.3p . 

Theorem 11.1. (Global existence) Let N > 3 be an integer, and assume that < c s < y/l/3, where c s 
denotes the speed of sound. Let {g^ v ,K^ v ,B^,^), (fJ,,v = 0,1,2,3), be initial data (not necessarily satisfying 
the wave coordinate condition or the Einstein constraints) on the manifold T 3 for the modified Euler-Einstein 

system (|5.4.2a[) - (I5.4.3p . and let Sn = 5 900+ i ; tv + Sg ,;N + Sh**;N + Sq^-^n be the norm defined in (I6.1.2fp . 
Assume that there is a positive constant c\>2 such that 

(11.2.1) -5 ab X a X b < g ab X a X b < ^5 ab X a X b , V(X\X 2 ,X 3 ) G M 3 . 

ci 2 

Then there exists a small constant eo, with < eo < 1, and a large constant C*, both depending on c\ and N, 
such that if e < eo and Sat(0) < C~ 1 e, then the classical solution (g^jd^Q) provided by Theorem \5.1\ exists on 
[0,oo) x T 3 , and 



(11.2.2) S N (t) < e 

holds for all t > 0. Furthermore, the time T max from the hypotheses of Theorem 15.^1 is infinite. 

Proof. See Remark 19.1.11 for some conventions that we use throughout this proof. Our proof relies upon a stan- 
dard bootstrap-style argument that ultimately relies on Theorem 15. 2\ i.e., we will make assumptions concerning 
the size of the energies and concerning g^ on a spacetime slab of the form [0,T) x T 3 , and we will use these 
assumptions, together with assumptions on the data, to deduce an improvement on the same slab. In effect, we 
will avoid the four breakdown possibilities of Theorem 15. 2| which will allow us to conclude global existence. 

Let e = Sjv(0), and let r} m i n be the constant defined in (|8.2.ip . We are assuming that I < C~ 1 e, where C* 
and e will be adjusted throughout the proof. To begin the detailed analysis, we first note that if e is sufficiently 
small and C* is sufficiently large, then by (|8.2.3ap and (|9.2,7ep . it follows that all of the hypotheses of Theorem 
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15.11 are satisfied. Therefore, there is a local solution (g, d&) existing on a (non-trivial) maximal interval [0, T) 
on which the following bootstrap assumptions hold: 

(11.2.3a) Sjv(i) < e, 

(11.2.3b) q 1 5 ab X a X b < e~ 2n g ab X a X b < Cl 5 ab X a X b , V(X\X 2 ,X 3 ) £ M. 3 . 

We remark that if e is sufficiently small, then (jll.2.3ap - (jll.2.3b| ) imply that the rough bootstrap assumptions 
(|7.0.1ap - (|7.0.1cp . with 77 = 7] m i n in (|7.0.1ap . hold on [0, T) x T 3 ; these assumptions were used the proofs of the 
two propositions in Section [9j 

By maximal interval, we mean that 

(11.2.4) T = sup {t > | The solution exists classically on [0, t) x T 3 , and (|11.2.3ap - (lll.2.3bp hold}. 

The remainder of this proof is dedicated to showing that T = oo if e is small enough and C* is large enough. 
Throughout the proof, we assume that e is small enough so that Propositions 19. 1.1} 19.2. 1\ 110.0. 1( and 111. 1.11 
are valid on [0, T). We will make repeated use of Proposition IIO.OTTI throughout this proof without explicitly 
mentioning it each time. 

We now address the bootstrap assumption (|11.2.3bp . with the intent of showing an improvement. First, we 
note that the assumption Sjv < e implies that 

(11.2.5) \\dt(e- 2n g jk )\\L°° = \\d t h jk \\ L oo < Cee~ qHt . 
We then use fjl 1.2.5f> to integrate in time from t = 0, concluding that 

(11.2.6) \\e~ 2n g, k (t, •) - g jk (-)\\ L «> < Ce. 

By (lll.2.ip and (lll.2.6p . it follows that if e is small enough, then on [0, T) x T 3 , we have that 

(11.2.7) ^-S ab X a X b < e- m g ab X a X b < ^±5 ab X a X b , V(X\X 2 ,X 3 ) G R 3 . 

Inequality (jll.2.7p shows that if e is small enough, then the bootstrap assumption (|ll,2.3bp can be strictly 
improved on the interval [0,T). 

To complete our proof of the theorem, we will show that if e is small enough and C* is large enough, then 
the bootstrap assumption (|11.2.3ap can be improved by replacing e with e/2; the primary tool for deducing an 
improvement is of course Proposition 111. 1.11 To begin our proof of an improvement of (jll.2.3ap , we use a very 
non-optimal application of Proposition Ill.lTTI with t\ = 0, deducing that on [0,T), we have that 

(11.2.8) Ejv(t) < Ejv(0) + f cE N (r)dT. 

Jt=0 

Applying Lemma 111.1.21 (Gronwall's inequality) to (jll.2.8p . using Sjv(0) = e, and using Proposition I10.0TT1 we 
conclude that the following preliminary inequalities hold on [0,T) : 

(11.2.9) S N (t) < Cee ct , 

(11.2.10) Ejv(i) < Cee ct . 

Remark 11.2.1. By modifying the argument in the last paragraph of this proof, Theorem 15.21 and inequality 
(jll.2.9p can be used to deduce that the time of existence is at least of order c _1 |ln(Ce)|, if I is sufficiently small. 

We now fix a time t\ G [0,T); t\ will be adjusted at the end of the proof. Roughly speaking, it will play 
the role of a time that is large enough so that the exponentially damped terms on the right-hand sides of the 
inequalities of Proposition Ill.lTTI are of size <C e. To estimate Egq>-^(t) on [t\,T), we simply use (jll.2.3aj) and 
(|11.2.10p to estimate the two terms on the right-hand side of (jll.l.la[) : 
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(11.2.11) E 2 ^. N {t) < E^. N (tx) + Ce 2 f e~ qHr dr < C{ie ct ' + ee~ qHh / 2 } 2 . 

Jr=ti 

Using (I11.2.10P again to estimate Eg^-^(t) on [0, ti], we thus conclude that the following inequality is valid on 
[0,T): 

(11.2.12) E d *. N (t) < C{ie ct > + ee~ qHt ^ 2 }. 

Applying similar reasoning to inequalities (lll.l.lbj) and (jll.l.ldj) . we also have the following inequalities on 
[0,T): 

(11.2.13) E goo+1 . N (t) < C{ee ctl + ee~ qHt ^ 2 }, 

(11.2.14) E htt , N (t) < C{ie ctl + ee~ qHtl/2 }. 

To estimate E got: -N(t), we use ([11.1. left , the bootstrap assumption (|11.2.3ap . and (jll.2.14|) to arrive at the 
following inequality valid for t £ [ti,T): 

(11.2.15) E 2 go ^ N (t) < E 2 0tt . N (h) + f -AqHE 2 got . iN (r) + C{ee ct > + ee~ qHt ^ 2 }E g ^ N {r) dr. 

Jr=ti 

Applying Lemma 111. 1.31 to ([11 .2. 15[) . with y(r) = E_ got . N (T) and 6(r) = 4qH in the lemma, and also using 
(|11.2.1Up . we conclude that the following inequality holds on [0, T) : 

(11.2.16) E gom]N (t) < C{ie ctl + ee~ qm ^ 2 }. 

Adding (|11.2.12p . (jll.2.13|) . ([11.2.14p . and (jll.2.16)) . referring to definition ([6.1.21) . and using Proposition 
110.0. 1[ it follows that on [0, T), we have 

(11.2.17) S N (t) < C{ee ctl + ee" 9 ^/ 2 }. 

We now choose t\ such that Ce~ qHtl ^ 2 < and I such that Cee ctl < ^e, which implies that the following 
inequality holds on [0, T) : 

(11.2.18) S N (t) < X -e. 

We remark that in order to guarantee that the solution exists long enough (i.e. that T is large enough) so that 
t\ £ [0, T), we may have to further shrink eo; see Remark lll.2.U We also remark that the constant C* from the 
conclusions of the theorem can be chosen to be 4Ce ctl , where C is from the right-hand side of (jll.2.17p . 

We now claim that T = oo. We argue by contradiction, assuming that T < oo. Then by combining (I11.2.7P 
and (jll.2.18p , it follows that none of the four existence-breakdown scenarios stated in the conclusions of Theorem 
!5.2l occur: (1) is ruled out by the Sobolev embedding result H^oo + < Ce~ qQ S]\f, (2) is ruled out by (11 1 .2.T|) : 

(3) is ruled out by ([9.2.7eP ; and (4) is ruled out by the Sobolev embedding results ||goo + 1 1 1 c 1 + II^^OoIIl 00 < 

Ce~ qn S N , £f =1 {hojWcl + Ijftgti j M < C e^y N g|fc =1 \\dg jk \ \ L ~ < Ce 2n S N , \\8$ \\ L ~ < Ce^ n S N , 
together with inequalities (19.1. left . (j9.1.3c|> . (|9.1.7aj> . (I9.2.7M . (|9.2.7dj> . ( |9.2.7gP , fl9.2.11ej> . and ( |9.2.11gp . Also 



using the continuity of Sjv(i), it thus follows from Theorem 15.21 that there exists a 5 > such that the solution 
can be extended to the interval [0, T + 5) on which the bootstrap assumptions (jll.2.3ap - (jll.2.3bp hold. This 
contradicts the definition of T, which shows that T is not finite. 

□ 



Stability of Irrotational Euler-Einstein 

54 

11.3. On the breakdown of the proof for c s > 1/3 (i.e. s < l,x > 1). In this short section, we give a 
brief example of how our proof breaks down when c s > 1/3. If x > 1, we cannot use our previous reasoning to 
bound the following term, which is the last term on the right-hand side of (|6.3.5p : 

(11.3.1) - V ( e 2Kn (d t m ab + 2xiom ab )(d a d^)(d b d d ^)d 3 x. 

2 I ^ AT 

\a\<N 

Previously, we had split this term into two pieces, one of which is (x— l)co X^|<3|<tv Jr 3 mab (dadd&)(dbdd&) d 3 x, 
which, as is explained in our proof of Proposition 1 1 1 . ITT] could be discarded from the energy inequality ([ll.l.lap 
because it is non-positive when x < 1. Obviously, we can no longer discard this term when x > 1. Furthermore, 
even in the case x = 1, inequality (|9.2.11fp must be replaced with \\dtm^ k + 2uim^ k \\L' xi < C{e~ 2 ^S 2 N + 
Positive terms}. Ultimately, this fact can be traced to the fact that the L°° estimate for Zj from (|9.2.10bp must 
be replaced with ||^j||i<» < CSat. 

Using Proposition I1U.U.11 it can be shown that the net result in both the case x = 1 and the case x > 1 is 
that the term (|11.3.1|) leads to (as in our proof of Proposition lll.l.ip an integral inequality of the form 

(11.3.2) E^. N (t) < E%Q N (h) + positive terms + Ce [ E^, N (r) dr. 

Jr=t x 

Inequality (jll.3.2p allows for the possibility of growth of Eq^-^()\, i.e., unlike in the cases < c 2 s < 1/3, there 
is no e~ qHr factor with q > in the integrand. Therefore, this inequality provides no means to improve the 
bootstrap assumption Sj\r(i) < fi- 
ll. 4. Future geodesic completeness. In this section, we prove our second main theorem, which provides 
criteria for the initial data under which the global solutions provided by Theorem 111.11 are future geodesically 
complete. The theorem and its proof are based on Propositions 3 and 4 of [28]. 

Theorem 11.2. (Future geodesic completeness) Let N > 3 be an integer, and assume that < c s < a/1/3, 
wherec s denotes the speed of sound. Let ([0, oo) xT 3 , g^ u , 9^), (fi, v = 0, 1, 2, 3), be one of the FLRW background 
solutions described in Section^ Note that the background is a solution to both the un-modified system (13.3. 10ap 
- (I3.3.10bp and the modified Euler-Einstein system (I5.4.2ap - (15 .4.3j) . Let (g^ u , LC^ u ,d^,^>), ({J.,v = 0,1,2,3), 
be initial data for the modified Euler-Einstein system on the manifold T 3 that are constructed from initial data 
(T 3 , g~jk, Kjk, ^), (j,k = 1,2,3), for the un-modified system (13.3. 10ap - (I3.3.10bp (and which are viewed as 
a perturbation of the data for the FLRW solution) as described in Section \3. 2. H Let Sn be the norm defined 
in (I6.1.2fp . and let S g - t N, a norm for the metric components, be the quantity defined in (j6.1.2dj) . Assume that 
the data for the modified system are near that of the FLRW solution in the sense that Stv(0) < C~ 1 eo, where 6q 
and C* are the constants from the conclusion of Theorem \ll.l\ (note that Sn = for the FLRW solution). Also 
assume that the perturbed data satisfy the inequality (|11.2.ip . so that all of the hypotheses of Theorem \ll.l\ are 

satisfied. Let (A4 =^[0, oo) x T 3 , g^ u , d^&) be the global (to the future) solution (which by Proposition \5.6.1\ is 
a solution to both the un-modified system and the modified system) guaranteed by Theorem \ll.l\ and let 7 be a 
future- directed causal curve in M with domain [sQ,s max ) such that 7(so) = 0. Let 7^ denote the coordinates of 
this curve in the universal covering space of the spacetime (i.e. [0, 00) x M 3 ). Then there exist constants C > 
and e±, where < e\ < eo, such that if Sg-pf(t) < e\ for all t > 0, then 7 > 0, and furthermore, the length of 
the spatial part of the curve as measured by the metric = g^ satisfies 

I* Smax j 

(11.4.1) / \Jgab(no-/)j"jbds<C, 

J SO 

def 

where ir denotes projection onto spatial indices; i.e., ir 3 07 = 7^. We remark that the constants C and e\ can be 
chosen to be independent of~f. Additionally, if 7 is future-inextendible, then 7°(s) f 00 as s f s max . Finally, 
the spacetime {M^g^u) is future causally geodesically complete. 

Remark 11.4.1. Our smallness assumption concerns S g ^ because this is the quantity that controls g^ u ; i.e., 
the influence of the fluid plays no direct role in the conclusions of this theorem. 
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Proof. See Remark 19.1.11 for some conventions that we use throughout this proof, but as noted above in the 
statement of the theorem, the smallness assumption on Sat mentioned in the remark can be replaced with a 
smallness assumption on S g -N. We also note that 7°(s) can be identified with the spacetime coordinate t. 

We begin by showing that 7 > 0, where • = f Since 7 is causal and future-directed, and since dt is 
future-directed and timelike, we have that 

(11-4.2) ff^fV < 0, 

(11.4.3) 50 a7 Q < 0. 

Our first goal is to prove that if ei is small enough, then the following estimates hold: 

(11-4.4) g ab ^ b < C,(7 ) 2 , 

(11.4.5) S ab n b < Ce- 2n (f) 2 , 



where rj = Ce\ is from the right-hand side of (j7.0. lap . To this end, we use the Cauchy-Schwarz inequality, 
(|7.0.1cp . the Cauchy-Schwarz inequality again, and (|7.0.1b|) to conclude that 

(11.4.6) |2 90 a7°7 a | < V 1/2 (i°) 2 + rr 1/2 baff < + 

<r, 1/2 (i°) 2 +V 1/2 e-W gai j°j>. 
Combining (|11.4.2p . (|11.4.6|) . and the estimate —goo < 1 + rj, which follows from (|7.U.lap . we have that 

(11.4.7) g ab j a i b < -good ) 2 + \2goai°i a \ < (1 + v)(i°) 2 + V 1/2 (i°) 2 + ^V^^T* ■ 

It now easily follows from (jll.4.7p that if ei is small enough, then there exists a constant C v > such that 
(|11.4.4p holds. Inequality (|11.4.5p then follows from fjll.2.3b|) and (|11.4.4p . 

We now claim that | -y° | > 0. For if |7°| = 0, then inequality (jll.4.5p shows that X^a=o l^"! = ^' wn i cn 
contradicts All .4.3|> . We may therefore divide each side of (jll.4.6p by |7°| > and use (jll.4.4p to arrive at the 
following inequality: 

(11.4.8) |<7oa7 a | <^ /2 (l + ^ e - 2 ^)| 7 °|. 

Using (I11.4.8P and the estimate \goo\ > 1 — r], we conclude that if ei is sufficiently small, then 
(11-4.9) sgn(# a7 a ) = sgn(£oo7°), 

where sgn(y) = 1 if y > 0, and sgn(y) = — 1 if y < 0. Since (jl 1 .4.3[) implies that the left-hand side of f 1 1 1 . 4 . 9 [) is 
negative, and since goo < 0, we conclude that 7 > 0. 

We now show the estimate ([11.4. ip . First, by the hypothesis (|11.2.ip . the fact that g^ = gj^, and the estimate 
()7.0.1bp . it follows that 

(H.4.10) g ab n b < ^S ab j a i b < §e~ 2n g ab j a i b . 

Integrating the Square root of each side of ([11 .4. lOj) from sq to s max , recalling that e~^^°^ < Ce~ H ~ 1 °( s \ and 
using (jl 1.4.4j) . we have that 

(11.4.11) / V& 6 (7ro T )ffd s < / Ce- H ^°fds = -j ! -e~ H "' ds < C, 

J Sq " Sq J Sq 

which proves (jll.4.ip . 

We now show that 7 converges to infinity as s f s ma x if 7 is a causal future-inextendible curve. Since 7 > 0, 
it follows that either 7 converges to infinity as s f s max , which is the desired result, or that 7°(s) converges to a 
finite number. In the latter case, by flll.4.4[) . we also conclude the 7 J (s) converge to finite numbers as s "f" 
which contradicts the definition of future-inextendibility. 
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To show that the spacetime is future causally geodesically complete, we consider a future-directed causal 
geodesic 7. We recall that the geodesic equations are 7^ + B^i = ^> which in the case of \i = reads 

(11.4.12) 7° + r Q %7 a 7 /? = 0. 

To analyze equation (|1 1.4. 12j) . we will use the following estimates for the Christoffel symbols: 



(11.4.13) |r ° | < Ce ie -" m , 

(11.4.14) ir^-l < Ce ie ^ m , 
(H.4.15) |r/ fe -^ fc |<Ce ie ( 2 -^. 

We will prove the estimate f 1 11 . 4 . 1 5 f> : the estimates ()1 1 .4. 13j) - (jl 1 .4. 14j) can be shown similarly. To begin, we 
use the definition (|3,0.2dp of T A, and the triangle inequality to obtain the following estimate: 



(11.4.16) 

1 111 

\Y 3 \ - iog jk \ < -Ig^Wdjgok + d k g 0j \ + -\g 00 + l||^5jjt| + -\d t g jk - 2ojg jk \ + -\g 0a \\djg ak + d k g ja - d a g jk \. 

Inequality (|11.4.15|) now follows easily from ([11 .4. 16|) . the assumption S 9 -n < ei, Sobolev embedding, (j9. 1.2cip . 

dSXSH , (I9.1.3bl), and (|9.1.3cp. 

We now use (|11.4.4j) . (|11.4.5|) . and (|11.4.13[) - (|11.4.15|) to arrive at the following inequality: 



(11.4.17) |r |(7°) 2 + 2|r ° a7 Vl + |r a ° fe - cog ab \n b < Ce ie -^°( s )(7 ) 2 , 

where it is understood that both sides of (jll.4.17p are evaluated along the curve 7(5). Using (|11.4.12p . (jll.4.17p . 
and the positive-definiteness of the 3x3 matrix gj k , it follows that 



(11.4.18) 7° = -r °o(7°) 2 " 2r a7 °7 a - r o 6 7 a 7 6 < K \(i°) 2 + 2|r ° a7 °7 a | + \Y a ° b - ug ab \n b 

Since we have already shown that 7 is positive if €\ is small enough, we may divide inequality (111.4. 18|> by it 
and integrate: 

(11.4.19) = < C, t\-^ns'W = f'^-^-'Us' 
\T{s )J J so 7 u (s') J so qH J SQ ds' 

~ qH 

We therefore conclude that 7 is bounded from above: 

(11.4.20) 7 < C. 
Integrating (jll.4.20p from sq to s, we have that 

(11.4.21) 7°( s )_ 7 0( So)= [ S f( s ')ds'<C\ S -s \. 

Js 

Finally, since we have already shown that 7 f 00 as s j s max , it follows from (jll.4.2ip that s max = 00. 

□ 
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In this section, we strengthen the conclusions of Theorem 111.11 by showing that gnv,g^ v , dn&, v = 0, 1, 2, 3), 
and various coordinate derivatives of these quantities converge as t — > oo. Because our strategy is to integrate 
bounds for time derivatives, we will lose at least one order of differentiability in our convergence estimates. 
Furthermore, we note that although our bootstrap assumptions were sufficient to close the global existence 
argument, they are far from optimal from the point of view of decay rates. Thus, at the cost of a few more 
derivatives, we will also revisit the modified equations and derive improved rates of decay compared to what 
can be directly concluded from the estimate Sat < e. In particular, the Counting Principle of Section [9] is not 
precise enough to detect these refinements. These results should be viewed as an initial investigation of the 
asymptotics; it is clear that more information could be extracted at the expense of more work. This theorem is 
analogous to |28[ Proposition 2]. 

Theorem 12.1. (Asymptotics) Assume that initial data (g^, K^ u , 3<&, x &), (fi, v = 0,1,2,3), for the mod- 
ified system (|5.4.2ap - f|5.4.3j) satisfy the assumptions of Theorem \11.1\ including the smallness assumption 
Sjv < C~ 1 e, where < e < eo- Let g^ v denote the inverse of g^ u . Assume in addition that N > 5, and let 
{g^y^d^) be the global solution launched by the data. Then there exists a constant ei satisfying < £2 < eo 

such that if e < €2, then there exists a Riemann metric g^\ (j,k = 1,2,3), with corresponding Christoffel 
symbols r\°$ , {i,j,k 



1,2,3), and inverse <7^\ on T 3 , and (time independent) functions d&^j, ^(00) 071 1" a 

_1 , and ^(oo) — ^ G H N ~ l where ^> is defined in 



Soo) 



9jk e H , g J (oo) 



such that g - k 

(14.1. 8ft . and such that the following estimates hold for all t > : 



gjk g R N^ QQ {oQ) g H l 



(12.0.1a) 
(12.0.1b) 



1 (°°) 1 1 ^ ri 

\9jk ~ 9jk\\H» ^ Ce > 



I jk 

15(00) 



g jk \\ H » < Ce, 



(12.0.2a) 
(12.0.2b) 
(12.0.2c) 
(12.0. 2d) 
(12.0.2e) 
(12.0.2f) 
(12.0.2g) 
(12.0.2h) 

(12.0.3a) 
(12.0.3b) 



-2Q 

2 gjk 

\\e 2n g jk 
\\e 2n g jk 



-2H 



9jk - 9jk \\h n S Cee 
g\f\\n^<Cee-^ 

Jk II ^ r<-2Ht 

'(00) lltf^- 2 ^ Cee 

'd t g jk -2ug\f\\ HN <Cee-* m 
) t g jk -2Log^ ) \\ HN -2 <Cee~ 2Ht 
\\e 2Q d t gi k + 2^^ <Cee-« m 
\\e 2Q d t gi k + 2^11^-. <Cee- 2Ht 



\„ u~ l r, ab r(°°)ii n^ a -i Ht 
ISoj - ti g (oo) L ajb \\ h n-z < Cee 



-2Q. 



9 (^) L ajt 

\\dtgoj\\H N -z 



< Cee- qHt 



(12.0.4a) 
(12.0.4b) 
(12.0.4c) 
(12.0.4d) 



||flbo + l||if" <Cee- qHt , 
||fldo + l||^-a <Ce(l + t)e- 2Ht , 
ll^oollffiv < Cee~ 9m , 
\\d t g m + 2uj{g m + 1)\\ h n-2 < Cee~ 2Ht , 
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[12.0.5a) \\e- 2Q K jk -ujg^ ] \\ H N-i <Cce~ qHt , 

[12.0.5b) \\e~ 2Q K jk -ug^WjjN-2 < C(l + t)ee 



In the above inequalities, Kj k is the second fundamental form of the hypersurface t = const. 
Furthermore, we have that 

(12.0.6a) ||e*°$$ - f (oo) \\ h n-i < Cee~ qH \ 

(12.0.6b) \\e Kn d t <5> - * (oo) \\ h n-2 < Cee' 2 ^-^ 111 , 

(12.0.6c) - d<$> {oo) \\ H N-i < Cee~ xHt , 

(12.0.6d) ll^oo) ||ffjv-i < Ce. 

Remark 12.0.2. We assume that N > 5 so that we can use standard Sobolev-Moser estimates during our 
proofs the improved rates of decay. 

Proof See Remark 19.1.11 for some conventions that we use throughout this proof. Furthermore, in our proofs 
below, we will introduce new energies, and the differential inequalities that we will derive are valid only under 
the assumption that the energies are sufficiently small; we don't explicitly mention the smallness assumption 
each time we make it. In the interest of brevity, we will only sketch the proofs of the estimates involving the 
improved decay rates. We also remind the reader of the conclusion (111. 2. 2ft of Theorem lll.l[ which is that 

Sat = S goo+ i- N + S g0t - N + S htt . N + S d 0. N satisfies S N (t) < e for t > 0. 



Proofs of (|12.0.1a| . (I12.0.1bl) . (|12.0.2ap . (I12.0.2cl) . (|12.0.2el) . and ( |12.0.2 g p : It follows from the definition (I6.1.2fp 
of Sjv that 

(12.0.7) \\dth jk \\ H N <Cee~ qHt . 

Integrating dthj k and using (112.0.7j) . it follows that for t\ < t2, we have that 



(12.0.8) ||M*2) - M*i)llff" < C 



ce 



-qHt 



Using (|12.0.8p and the fact that hj k = e 2n gj k , it easily follows that there exist functions g)^ (x 1 , x 2 , x 3 ) such 
that 

(12.0.9) ||e" 2 V - g^y* < Cee~ qH \ 

and such that 

(12.0.10) Wffi ~ 9jk\\ H » <Ce. 

We have thus shown (112.0. lap and (I12.0.2al) . Inequality (I12.0.2el) follows from (I12.0.2ap and (112.0.71) . 
To obtain the asymptotics for g^ k , we use (I9.1.5ap . which implies that 

(12.0.11) \\d t (e 2n g jk )\\ H N < Cee~ qHt . 

From ([12.0. lip and the estimates (|9.1.2ep . ()9.1.2f|) at t = 0, it follows as in the previous argument that there 
exist functions g^^Jx 1 , x 2 , x 3 ) such that 

(12.0.12) || e 2 V fe -^)||^<Cee-^, 
and such that 
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;i2.0.13) ll^-^II^^Ce, 
;i2-0.14) ll^ll^-i <Ce, 

;i2-0.15) \\fy\L«><C, 



where gi k = g jk \ t=0 . This proves (I12.0.1bp and (Il2.0.2cl) . ( |12.0.2g| then follows from (|12.U.2cj) and (112.0.111) . 
Furthermore, since 

(12.0.16) g aj 9ak + g 0] gok = 5i 

and since Proposition EUl S N < e, (I9.1.2hj) . and (|9.1.2ij) imply that 
(12.0.17) \\9 03 90k\\H^<Ce 2 e-^ m , 

it follows that are the components of a Riemannian metric g^°°^ and that are the components of its 
inverse p (oo) . 

Proofs of (|12.0.4ap and (|12.0.4cp : The estimates (|12.0.4ap and (|12.0.4cp follow trivially from definition (|6.1.2ap . 

Proofs of (|12.0.6ap . ()12.0.6c|) . and ()12.0.6dp : To prove (|12.0.6ap . we first recall equation (|5.4.3p . which can be 
re-expressed as follows: 

(12.0.18) dtfe^dt® - tf) = e^A'g,, 

where ^> is defined in (|4.1.8p and 

(12.0.19) A' d<s> = -m ab d a d b <S> - 2m 0a d a d t $ + A 9$ . 

Using Proposition EH1 the definition (I6.1.2fp of Sat, Sobolev embedding, (I9.2.6ep . (j9.2.7dp . fl9.2.7fl) . and (I9.2.7al) . 
it follows that 

(12.0.20) \\e" n A' d(S> \\ H N-i < Cee- qHt . 

As in our proof of (I12.0.9|) . it easily follows from (|12. 0.181) . (112.0.201) . and the initial condition \\e KQ( -^d t ^{0, •) - 
^llifjv < Ce that there exists a function ^(^(x 1 , x 2 , x 3 ) with ^(oo) — ^ £ H N ~ l such that 

(12.0.21) He^W, •) - * {oo) ll^-x < Cee"^, 

(12.0.22) ||* (oo) - ^ll^-x < Ce, 

which proves (I12.0.6ap . 

Furthermore, the bound Sjy(i) < e implies that 

(12.0.23) \\d t d$\\ H N-i < Cee-* Ht . 

It follows easily from (|12.0.23[) and the initial condition ||9$(0, -)||^jv < Ce that there exists a function 
9<l ) ( 00 )(x 1 , x 2 , x 3 ) (see Remark fl. 0.1 I) satisfying d^^) G H 1 ^" 1 such that 



(12.0.24) •) - a* (oo) ||^-i < Cee- xm , 

(12.0.25) l|^(oo)ll^-i <Ce, 
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which proves (I12.0.6cl) and fj!2.0.6dj) . 

Proof of (|12.0.5a|) : We first observe that N, the future-directed normal to the surface {t = const}, can be 
expressed in components as 

(12.0.26) A> = -(-5 00 r 1/2 <?° M - 

By the definition of K, it thus follows that 

(12.0.27) K jk = g ak DjN a = -(d 3 [{-g™)- l / 2 g^])g ak - (-g 00 )^ 2 g 0a T jka . 

The dominant term on the right-hand side of (|12.0.27p is the one involving Tj k Q : it follows from Corollary 
IA-41 with v = <7°° + 1 and F(v) = (1 — v)~ 1 / 2 = (— g 00 )" 1 / 2 in the Corollary, Proposition IA-51 the definition 
(I6.1.2fp of S N , Sobolev embedding, (19.1. left . (I9.1.2cl) . (I9.1.2hl) . (I9.1.3ap . and definition (|12.0.53p that 

(12.0.28) 

„-2Q||7 r , / f)0\-l/2 no-p II _ „-2Q\\ /a f/ J30\-l/2 rial \ / J0\-i/2 > r II ^n co -i m 

e \\K jk + (-g ) g 1 j k o\\ HN -i - e || _ l°j|A _ 3 ) 9 \)9ak-{-g ) 9 *- jka\\ H N-i < Lee 

(12.0.29) 

i-9 00 )- 1/2 9°%ko + % ik nN _ t < \e-^\\{-9 Q0 r l/2 9 m (d mk - d k g Qj )\ 



e -2U 



, 1 2f2|| r/ JOl-l/2 JO , ila . || 

+ 2 e ) 9 + l \Ot9jk\\ H N-i 



< Cee~ qHt . 

Inequality (I12.0.5ah now follows from combining (I12.0.2el) . (I12.0.28p . and (112.0.2911 . 



□ 



Proofs of (j!2.0.3aj) - (j!2.0.3bj) : Our proofs of (I12.0.3ap - (112.0.3bp are based on two refined versions of the 
energy inequality (jll.l.lcj) . The main point is that the even though the energy E_ gQt . N defined in (|6.2.8bp allows 
us to efficiently close the bootstrap argument of Theorem 111.11 there is room for improvement. In particular, 
Theorem 111.11 only allows us to conclude that E_ gQt . N < Sn < Ce, which implies that \\dtgoj \\jjn < ee^ l ~ q ^ Ht 
and Hfi'oill//^ < ee ( - 1 ~ q ^ m . As we will see, it is possible to improve these estimates by a factor of e^ q ~ l ^ Ht . This 
is a preliminary step that we will need in our remaining proofs. This improvement is based on the following 
simple identity: 

(12.0.30) g Q i = -—g aj g 0a . 

goo 

We begin our proof of the improvement by defining a new energy (£ 90>> .jv-i f° r the 9oji (j = 1> 2, 3), by 



3 

( 12 -0-31) <£L;iV-i= £ E^,A.)[ 9 * 1 9(3 59 o J )], 

\S\<N-lj=l 

where £ 2 ^ q Sq ^[dd90j,d(da90j)] is defined in (|6.2.2|) . and the constants 70*, ^0* are defined in Definition 16.2.11 
Note that the scaling in f| 1 2 . . 3 1 [) differs from the scaling used in definition ()6.2.8bp by a factor of e( 1 ~~ q ^. 
Furthermore, we are using an A?" — 1 st order energy, rather than an N th order energy, because we will make 
use of the improved rates of decay for lower derivatives that are already discernible from the fact that S^v < e 
(consider e.g. inequality (|9,1.6bp ). Now using (|10.0.1ap . we have the following comparison estimate: 
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3 

(12.0.32) C-^jv-i < H^Oill^-i + e^WdgojWHN-! + C {j0t) \\g 0j \\ H N-i < C^.^. 

3=1 

Using the definition of <£_ g0t -N_i and the comparison estimate, it follows that the energy inequality (|6.2. 12b|) 
can be replaced with 

(12.0.33) 

d 3 3 

^QdL^-i) - - r lo*Hg 2 got . N _ 1 + C£ got . N ^ 1 ||ff ab r aife || H jv-i + C^ got . N _ 1 || AojU^jv-i 

3=1 3=1 
3 3 

+ C£g *;N-l Y\^9^ d s)90j\\L^ + X] X] ll A £;(7o,,5o,)[ 5 a5'Oj,9(0 5 5oi)]||Li- 

|cf|<JV-lj'=l |a|<JV-lj'=l 

We now claim that the following improvements of (|9.2.1bp . (|9.2.1e|) . (|9.2.2b|) . ()9.2.3b|) . and (|9.2.4bp hold: 



(12.0.34) IIAa.OjIIh"-! < Cte-* m g^. N _ x + Cee"***, 

(12.0.35) HAc^ll^-i <Cee-' lHt g^. N _ 1 , 

(12.0.36) ||A 0j || H iv-i < Ce-' 111 *^.^ + Cee~ qH \ 

(12.0.37) H^^Hl* ^Cee-^^.jv.i + Cee"^*, (|a| <N-1), 

(12.0.38) ||A £ . (70s! ^ ) [a 5ff0 i,5(^)]|| i i <Ce e -^ 0if . Ar _ 1 + Ce-^ 0)> . iV _ 1) (|a| < N - 1). 



These improved estimates can be derived using the same methods as in our proofs of the original estimates, 
together with the identity (|12.U.3U|) . More specifically, in our proofs of (|9.2Tb|) . (|9.2.1ej) . (|9.2.2b|> . (|9.2.3bp . 
and (|9.2.4bp . we used the estimates ||c^<7qj||.h-jv < e ^ 1_9 ^i5g 0j ,-iVj llffojll/f^ < e^ q ' n S_ go ^. N , and H^ojllffw < 
e( 2_9 ) n 5 .jy, which follow directly from the definition of S_ g0t .j\f, together with (|9.1.2hp . which reads 
Ce~( 1+q ^ n S g] N- However, whenever it is convenient, these estimates can be replaced with 



(12.0.39) \\9t90jWHN-i < 
(12-0.40) bojWHN-i^Cl^.^, 

(12.0.41) \\dg 0j \\ HN -,<Ce n l g0t . N _ 1 , 
(12.0.42) ||0 Oj ||tf"-i <Ce~ 2n l qot[N _ 



respectively, where (112.0.391) - (|12.U.4ip follow from (112.0.321) . while (|12.U.42p follows from applying (112.0. 321) . 
Proposition IA-51 Sobolev embedding, (|9.1.2e|) . and (|9.1.2fp to the identity (|12.0.30[) . We remark that the 
Ce~ qHt S_ g0f . N _i term on the right-hand side of (|12.0.36[) arises from from e.g. the {u — H)dtgoj term on the 
right-hand side of Q5.4,4jp . Similarly, the Ce~ qHt £^ /Qt . N _ 1 term on the right-hand side of (|12.0.38[) arises from 
the e.g. the (H — ui)g ab (d a v)(dbv) term on the right-hand side of (|6.2.6p . See |28| Section 14] for additional 
details on these improved estimates. 

Now using (I9.1.6bp . (|12.U.33I) and (|12. 0.361) - (112.0.381) . we argue as in our proof of (|11.1.3|) to deduce the 
following inequality: 

(12-0-43) < -TfrHg^^_ 1 + Ce-« m g^^_ 1 + Ceg^. N _ 1 , 

where the Ce(o_ go ^. N _ 1 term on the right-hand side of (|12.0.43p arises from applying inequality (|9. 1.6b[) to the 
second term on the right-hand side of (|12.0.33|) . Integrating (|12.0.43|) from to t, using the smallness condition 
— go*-N-i(®) — C £ l° r sma ll t, and applying Lemma 111. 1.31 for large t (the Ce~ qHt £^ /Qt . N _ 1 term is dominated 
by the —r]o*H^] j0t . N _ 1 term for large t), we conclude that the following bound holds for all t > : 
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(12.0.44) l 90 *;N-i < Ce. 

This completes the proof of our preliminary improved estimate. 

We are now ready for the proofs of (|12.0.3ap and (|12.0.3bp . Defining 

(12-0.45) ^o.-tf-yVS?' 
and recalling that goj is a solution to 

(12.0.46) n g g 0j = 3Hd t g 0j + 2H 2 g 0j - 2Hg ah T a]h + A 0j 

it follows that Vj is a solution to 



(12.0.47) U g Vj = mdtVj + 2H 2 Vj + Aj, 

where 

(12-0.48) A, = A 0i + 2H(gf oo) T ( ^ - g ah T a]b ) - g lm d^g^T^) . 

To estimate Vj, (j = 1, 2, 3), we will use the energy 

(12-0.49) H ;N -3 = ^ Ht jZ £ lM^ v ^ d ^ v ^ 

|5|<iV-3i=l 

where S? s ^ ^[d^Vj,d(d^Vj)] is defined in (|6.2.2p . and the constants 70*, ^0* are defined in Definition 16.2.11 
Note that in response to the last term on the right-hand side of (|12.0.48p . we have further reduced the number 
of derivatives in the definition of our energy by two. Using (110.0. lap , it follows that 

3 

(12.0.50) C-^.jv.g < + C ( ^ m \\ Vj \\ HN - 3 + e^^WdvjW^-a < C^. N _ 3 . 

Arguing as in our proof of (|12.0.33|) . we have that 

(12.0.51) 

d 3 
^GdU-a) < (29 - Vo*)Hll ;N _ 3 + Ce^ Ht S Vt . tN ^ l|Ai||^-a 

3 3 

|cf|<JV-3i=l \a\<N-3j=l 

Note that the inequality (I12.0.5ip does not have a term corresponding to the ||(? a6 r a jb||jyjv-i term in (|12.0.33|) ; 
the remnants of this term are present in || A,.||^jv-3. 

We will now estimate ||Aj||^iv-3, our goal being to show that it decays like ee~ qHt . We begin by estimating 
the term 2Hg ab T a ji, from the right-hand side of (|12.0.48p . Let us first recall the definitions of the lowered 
Christoffel symbols and T^j, corresponding to g^°°^ and g respectively: 

(12-0.52) rg? d = f \{d i9 ^ + d k9 \f - dj g™), 

(12.0.53) T^ av = f ^{d^gav + d u g m - d a g^ u ). 
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Using Proposition^!! the definition (I6.1.2fl) of S N , Sobolev embedding, (|12.0.1al) . ()12.0.1bl) . (I12.0.2al) . (I12.0.2cl) . 
(I12.U.15P . (I12.0.52p . and (|12.0.53p . we conclude that 

(12.0.54) \\g ab (t,.)T ajb (t,.) -g^T^Wn^ < \\e 2n 9 ab (t, •) - 9^1^-4^^, .)|| HJV -! 

_i_ iu ab II n -2^r /'f \ T^t°°) ll 

+ ||a< ) ||^- 2 ||e- 2n r aj6 (t,.)-r^ ) ||^ 

(12-0.55) ll^r^ll^-i^Ge. 

We now estimate the term H~ 1 g lm did m (g®^r^° l j) from the right-hand side of (112.0.48p . By Proposition 
bV5l (|9.1.2e.p . (j9.1.2fp . and (112.0.551) . it follows that 

(12-0.56) Wg^Oidmig^T^W^-s < Cee~ 2Ht . 

Applying the estimates (112.0.36j) . (I12.0.44p . (j!2.0.54j) . and (I12.0.56|) to the terms in (112.0.48j) . we conclude the 
desired estimate for ||Aj||^jv-3 : 

(12.0.57) \\Aj\\ H N-a < Cee~ qHt . 

In addition, arguing as in our proof of (19.2,3bp and ()12.0.38p . and in particular making use of the improved 
estimates (|12.0.39[) - (j!2.0.42j) and (|12.0.44D . it can be shown that 

(12.0.58) H[n 5 ,daKllL2 < CSquareone- qHt £ Vt . N _ 3 + Cee~ qH \ (|a| <N-S), 

(12.0.59) WA^s^v^did^)}]^ < Cee~ 2qHt ^ N _ 3 + Ce~^ m ^ N _ 3 , (|a| < JV - 3). 

Using (|12.U.5ip . (|12.U.57I) . (|12.U.58I) . and (I12.0.59D . it follows that 
(12-0.60) ^(^U-s) ^ ( 2 9 " m*)Hll ]N _ 3 + Ce-^<^_ 3 + Ce^. N _ 3 . 

Using the fact that 2q — 770* < 0, we argue as in our proof of (|12.0.44|) to conclude that the following inequality 
holds for all t > : 

(12-0.61) l Vt . N _ 3 < Ce. 

Inequalities (|12.0.3a|) and ()12.0.3bp now follow from the comparison estimate ([12.0. 50[) and from (j!2.0.6ip . 

Proofs of (I12.0.2bp . (I12.0.2dl) . ()12.0.2fp . (I12.0.2hp . (I12.0.4bl) . and (I12.0.4dl) : 

We begin by using equations (|5,4,2ap and (|5.4.2cp . together with the identity g 00 + 1 = g 00 (goo + 1) + g 0a go a 
to derive the following equations: 

(12.0.62) d?[e 2n (g o + l)] = -udt [e 2n {goo + 1)] + 2ue 2n (g 00 + 1) + e 2n A' 00 , 

(12.0.63) d t [e 2n h jk ] = -ue 2n d t h jk + e 2n A' jk , 

where 
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(12.0.64) 

A'oo = (g 00 )-^ ~ 2g 0a d a d t g 00 - g ab d a d b g m + 5a; [g 00 (g 00 + 1) + g 0a g a\ d t9oo + 6a; 2 [g 00 (g o + l) 2 + <? 0a <70a(<7oo + 1)] 

+ 2(< 7 00 )- 1 {a Ai oo + A Ci00 + [/(£>£) - /(£>$)] _ ( goo + _ ^(ffoo + 1)^ +1 }, 

(12.0.65) 

= (5 00 )" 1 ! - 2g 0a d a d t h jk - g ab d a d b h jk + 3u;[<? 00 (<?oo + 1) + 9° a 9oa]d t h jk } 

+ 2( 5 00 )- 1 {e- 2n A yijfc - [ 5 °°( 5oo + 1) + g 0a g 0a ]a s+1 h jk - 2ug 0a d a h jk - 2e~ 2n a s (d J <l>)(d k ^) 

Using the fact that Sat < e, the improved estimates (whenever they are convenient) (|12.0.39[) - (|12.0.42|) . 
(|12.0.44p . and the Sobolev-Moser type inequalities in the Appendix, we derive the following inequalities: 

(12.0.66) \\A' 00 \\ h n-2 < Cee~ m + Cee~ qHt { \\g m + 1\\ h n-i + \\d t g o\\ H N-2 + \\dth jk \\ H N- a }, 

(12.0.67) IIA^H^-2 < Cee~ 2n + Cee- qHt \\d t h jk \\ HN -2 . 

Since the derivation of the above inequalities is similar to many other estimates proved in this article, we have left 
the tedious details up to the reader. However, we remark that the ee~ qHt \\goo + lll/f^- 2 term on the right-hand 
side of (|12.0.66p arises from e.g. the first term on the right-hand side of (|5.4.12ap . that the ee~ f/ ^* 1 1 c^i^oo 1 1 i 1 /^— 2 
term on the right-hand side of (|12.0.66[) arises from e.g. the first term on the right-hand side of (|5.4.10ap . that 
the ee~ qHt \\dthj k \\ H N-2 term on the right-hand side of (|12.0.66[) arises from e.g. the last term on the right-hand 
side of (|5.4.10ap . and that the ee~ qHt \\dthj k \\jjN-2 term on the right-hand side of (|12.0.67[) arises from e.g. terms 
on the fourth and fifth lines of (|5,4.10cp . Furthermore, we remark that we are using H N ~ 2 norms because of the 
presence of the terms on the right-hand sides of (|12.0.64[) - (|12.0.65|) that contain second derivatives, e.g. the 
term <9 a <9f>gooi by examining e.g. the definition (|6.1.1ap . we see that the H N ~ 2 norm of such a term has a more 
favorable rate of decay than its H N ~ l norm. We need this additional decay to deduce (|12.0.66[) - (j!2.0.67|) . 

Now in order to derive our desired inequalities, it is convenient to introduce the following non-negative energies 
^ e 2fi (900+1);7V „2, £&[e«»(floo+l)]iJV-2> and which are defined by 



(12.0.68a) l 2 e 2n igoo+1)]N _ 2 = E / {d s {e 2n (g 00 + l)]) 2 d'x = \\e™(g 00 + 

(12.0.68b) £ft[e*Wi)];tf-2 = E /. (^[e 2n (,?oo + l)]) 2 d 3 x = \\d t [e 2n (g 00 + 1 



\a\<N-2' 



3 



(12.0.68c) l 2 £ 2n 9thtt . N _ 2 = EE/ e 4n (d d d t h jk ) 2 d 3 x=^2\\e 2n d t h jk \\ 2 HN _ 2 . 

\a\<N-2j,k=l J 1 j,k=l 

Observe the following simple consequence of the definitions (j!2.0.68ap and (|12.0.68bp : 

(12.0.69) e 2n \\d t g o\\ H N-2 < C{£^m (soo+l y )N _ 2 + £s t [e in {goo+l)];N-2} ■ 

Now from (|12.U.66I) and (|12.U.67j) . and (|12.0.69|) . it follows that 

(12.0.70) e 2n ||A' 00 ||^iv-2 < Ce + Cee~ qHt {£ e 2tl( goo+1 y iN ^ 2 + ^d t [e^(g 00 +l)};N~2 + ^e™d t h„;N-2}: 

(12.0.71) e 2n \\ A' jk \\ HN - 2 <Ce + Cee~ qHt l_ e 2n dth ^. N _ 2 . 
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We may therefore use equations (|12.0.62j) - (|12.0.63|) . together with (|12.0.70j) - (|12.0.71|) to derive the following 
system of differential inequalities, which is valid if £L e 2n {g 00 +i)-,N-2i£Ld t [e 2n {g w +i)}-,N-2-> an( i ^.e 2n d t h tt ;N-2 are 
sufficiently small: 



(12.0.72a) 

-j t {^ e 2n {gm+l) . N _ 2 ) < 2^ e 2n ( 3 00+1);A r_ 2 ^a t [ e 2n(goo+l)];A r -2' 
(12.0.72b) 

^{£ 2 dt [ e 2n {goo+1)];N _ 2 ) < -2w^ t [ e an( S0o+ i)] i j V _ 2 + 2e 2 ^^ af [e 2n ( 500 +i)] ; 2 1 1 A Oo\\h"- 2 

+ 4|w|^ e 2n( S00+1 ).jv_2^[ e 2n( S00+1 )].jv-2 

< -2w^ t[e2 n( ff00+1 )] ; jv_2 + C e ~ 9m ^ e 2n (9oo+l);A r -2^i[e 2n (3oo+l)];A r -2 + ^ e ^9 t [e 2n (goo+l)];JV-2 
+ Cee 9m ^ 9t [ e 2n( 900+1) ] ;A r_2{^ e 2n( 300+1 ) ;A r_2 + &d t [ e ™(g 00 +l)];N-2 + ^U 2n dth**;N-2\ i 

(12.0.72c) 

^ {^2n 9th ^. N _ 2 ) < - 2uj £ 2 e 2n dth ^. N _ 2 + 2e 2n g e 2n dthtt . N _ 2 \\A' jk \\ H N-2 

< -2ug 2 e 2n dth ^. N _ 2 + Ce£^2n dthtt . N _ 2 + Cee~ qHt ^_ 2 e2 n dth ^. N _ 2 . 

We remark that (|12.0.72ap follows from a simple application of the Cauchy-Schwarz inequality, after bringing the 
time derivative under the integral over T 3 . Using the initial conditions <C e 2«( 900+ i) ; 7v-2(0) — C £ > <£dt[e 2n (go +i)];7v-2(C)) < 
Ce, and <o_e2ng thtt . N _ 2 (0) < Ce, and assuming that e is sufficiently small, we apply a Gronwall-type inequality 
to the system (112.0.72ap - (I12.0.72cj) . concluding that the following inequalities hold for all t > : 

(12.0.73a) ^( SO o+i);JV- 2 < Ce(l + t), 

(12.0.73b) ^a t [ e 2n( 90 o+i)];A r -2 ^ ^ 

(12.0.73c) i^dth^iN-z < Ce. 

Inequalities (|12.0.4bp and (|12.0.4dp now follow from fll2.0.73ap and (I12.0.73bp . Similarly, (112.0. 73cl) implies that 

(12.0.74) \\d t h jk \\ HN -2 < Cee~ 2n , 

from which (I12.0.2bl) and (I12.0.2fp easily follow. 

To obtain ()12.0.2dp . we first use the improved estimates (I12.U.39P - (|12.0.42p . (H2.U.44p . and (|12. 0.741) to 
upgrade the inequality (|9.1.9p as follows: 

(12.0.75) \\d t g jk + 2u;g jk \\ H N-2 < Cee~ m . 

(112.0.2dl) then follows from (I12.0.75p as in our proof of (I12.0.2cp . (112.0.2hl) then follows from (112.0. 2dl) and 
(I12.0.75p . 

Proof of (|12.0.5b| ) : To prove ()12.0.5bp . we use the identity g 00 + 1 = ^~[(<?00 + 1) ~~ 9° a 90a] an d the improved 
estimates (I12.0.4bp . (I12.0.4dp . (112.0.391) - (I12.0.42p . (I12.0.44p . and (I12.0.74p to upgrade inequalities (I12.0.28p - 
(I12.0.29P as follows: 



(12.0.76) e- 2n \\K jk + ( VT V V°I>|U- 2 < Cee~ 2H \ 

(12.0.77) e- ffl ||(Y)-n °r jfe o + ^4,_ 3 <^(l + ^- 2fft 

Combining (|12.0.2f) . (|12. 0.761) . and (112.0. 771) . we deduce (I12.0.5bp (as in our proof of (I12.0.5al) ). 
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Proof of (|12.0.6bD : To prove (I12.0.6bl) . we use the improved estimates (I12.0.4bl) . (I12.0.4dp . (112.0.391) - (112.0.421) . 
(I12.U.44|) . and (I12.0.74p to upgrade inequality (112.0.201) to 

(12.0.78) \\e Kn A' d ^\\ HN -2 <Cee~ 2{x -^ m . 

(I12.0.6bl) then follows from (I12.0.78p . as in our proof of (|12.0.6al) . 
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Appendices 

A. SOBOLEV-MOSER INEQUALITIES 

In this Appendix, we provide some standard Sobolev-Moser type estimates that play a fundamental role in 
our analysis of the nonlinear terms in our equations. The propositions and corollaries stated below can be 
proved using methods similar to those used in |19[ Chapter 6] and in |21] , The proofs given in the literature 
are commonly based on a version of the Gagliardo-Nirenberg inequality [26j , which we state as Lemma IA-11 
together with repeated use of Holder's inequality and/or Sobolev embedding. Throughout this appendix, we 
abbreviate LP = LP(T 3 ), and H M = H M (T 3 ). 

Lemma A-l. // M,N are integers such that < M < N, and v is a function on T 3 such that v G 
L°°,\\d( N ) v \\ L 2 < oo, then 

(A.l) M \\\ L2N/M < C(M,iV)||«||i~^||aWt;||J. 

Proposition A-2. Let M > be an integer. If {f a }i< a <z are functions such that v a G L°°, ll^^Vallz, 2 < oo 
for 1 < a < /, and c?i, ■ • ■ c?j are spatial derivative multi-indices with \a±\ + • • • |d?/| = M, then 

l 

(A.2) ||(^i)(^ 2 ) ■ ■ ■ (d dlVl )\\ L2 < C(l, M) (\\d {M) v a \\ L2 J] ||« 6 || L cc) . 

a=l b^a 

Corollary A-3. Let M > 1 be an integer, let & be a compact set, and let F £ C^ 1 {&) be a function. Assume 
that v is a function such that v(T 3 ) C ^ and Bv 6 H M ~ X . Then B(F o v) G H M ~ 1 , and 

M 

(A.3) \\B(F o v)\\ h m-i < C(M)\\dv\\ HM -i Y \F { %\Ml~- 

i=i 

Corollary A-4. Let M > 1 be an integer, let & be a compact, convex set, and let F G C^' 1 (&) be a function. 
Assume that v is a function such that v (T 3 ) C ^ andv — v G H M , where v G A is a constant. Then Fov — Fov G 
H M , and 

M 

(A.4) \\Fov -Fo V \\ h m < C(M){ \FW\si\\v - v\\l> + \\8v\\ H u-i IF^UM] 1 ^}. 

i=i 

Proposition A-5. Let M > 1, 1 > 2 be integers. Suppose that {i> a }i<a<z are functions such that v a G L°° for 
1 < a < I, that vi G H M , and that Bv a G R M - X for 1 < a < I - 1. Then~ 

l-i l-l 

(A.5) \\vtV2 ■ ■ ■ Vi\\ h m < C(1,M){\\vi\\ h m JJ ||fa|U°° + Y \\d v a\\HM-l Y{ \\ v b\\L°°y 

a=l a=l b^=a 

Remark A.l. The significance of this proposition is that only one of the functions, namely vi, is estimated in 
L? . 
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Proposition A-6. Let M > 1 be an integer, let 8. be a compact, convex set, and let F G Cff (&) be a function. 
Assume that V\ is a function such that ui(T 3 ) C that Bv\ G L°°, and that d^ M ^v\ G L 2 . Assume that v 2 G L°°, 
t/iai 5^ M_1 )f2 G L 2 , and Zei a be a spatial derivative multi-index with with \a\ = M. Then d$ ((F o vi)v 2 ) — (F o 
vi)d$V2 G L? , and 

\\d d ((F o vi)v 2 ) - (F o «i)5 5 « 2 ||i,2 

A/ 

1=1 

Remark A. 2. The significance of this proposition is that the M th order derivatives of v 2 do not play a role in 
the conclusions. 
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